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Abstract 



During the last decade a lot of efforts have been undertaken to ex- 
plain the unusual normal state properties of high temperature super- 
conductors (HTS) in the framework of unconventional theories based 
on strongly interacting electrons, pre-formcd Cooper pairs, polaron 
mechanism of superconductivity etc. A different approach to this 
problem would be to develop the perturbation theory for interact- 
ing electrons in the normal phase of strongly anisotropic supercon- 
ductors without specifying the origin of this interaction. The Cooper 
channel of interelectron interaction is equivalent to the superconduct- 
ing fluctuations which are unusually strong in HTS. We show that 
the peculiarities of such systems not only lead to the increase of the 
magnitude but are also frequently responsible for the change of the 
hierarchy of different fluctuation effects and even of the sign of the to- 
tal corrections. As a result the fluctuation contributions can manifest 
themselves in very unusual forms. 

The flrst and well known result is that that now one has the "pre- 
formed Cooper pairs" automatically, from ab initio calculations: tak- 
ing into account thermal fluctuations leads to the appearance of a non- 
zero density of fluctuating Cooper pairs (with finite lifetime) within 
layers without the establishment of long range order in the system. 
The fluctuation Cooper pair density decreases with temperature very 
slowly (~ In j7}rp in 2D case) . The formation of these pairs of normal 
electrons leads to the decrease of the density of one-electron states 
(DOS renormalization) at the Fermi level, and this turns out to be 
the key effect in our discussion. 

The DOS contribution to the most of characteristics is negligible 
for traditional superconducting materials. However it becomes domi- 
nant when highly anisotropic materials are discussed, and therefore is 
very important in HTS, especially when transport along the c-axis is 
considered. We analyze the role of the DOS fluctuations in the prop- 
erties of HTS and show how, taking into account this effect, many 
puzzling and long debated properties of HTS materials (such as the 
steep increase of the electrical resistivity along the c-axis just above 
Tc, the anomalous magnetoresistance, effects of the magnetic held on 
the resistive transition along the c-axis, the c-axis far infrared absorp- 
tion spectrum, NMR characteristics around the critical temperature 
etc.) can be understood leading to a simple, consistent description in 
terms of the fluctuation theory. 
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1 Introduction 



There are no doubts that the puzzhng anomahes of the normal state proper- 
ties of high temperature superconductors (HTS) are tightly connected with 
the physical origin of superconductivity in these materials. Among them are: 

- a peak in the c-axis resistivity above Tc followed by a decrease to zero 
as temperature is decreased |1|, ; 

- the giant growth of this peak in the presence of an external magnetic 
field applied along the c-axis and its shift towards low temperatures 0; 

- the giant magnetoresistance observed in a wide temperature range above 
the transition ^, ^; 

- the deviation from the Korringa law in the temperature dependence of 
the NMR relaxation rate above Tc 0; 

-the opening of a large pseudo-gap in the c-axis optical conductivity at 
temperatures well above Tc [§, § ; 

- the anisotropic gap observed in the electron spectrum by angular re- 



solved photo-emission experiments |T0 



the gap-like tunneling anomalies observed already above |jTT|, O, 



14, 15, 16 



the anomalies in the thermoelectric power above Tc |T^, |T^ . 



- the anomalies in the Hall effect above Tc p9| , P0| , pT| p^ . 

- the anomalies in the heat transport above Tc [pSj p4| . 

These effects have been attributed by many authors to the opening of a 
"pseudo-gap". Naturally this has led to numerous speculations about the 
physical origin of such a gap. 

During the last decade a lot of efforts have been undertaken to explain 
the unusual normal state properties of HTS materials using unconventional 
theories of superconductivity based on ideas of spin-charge separation, pre- 
formed Cooper pairs, polaron mechanism of superconductivity, etc. (see for 



instance |25|). They have been widely discussed and we will not overview 
them here. In the case of HTS with a well developed Fermi surface (i.e. in 
the optimally doped or overdoped part of the phase diagram) one can ap- 
proach this problem from another side, namely to develop the perturbation 
theory for interacting electrons in the normal phase of a strongly anisotropic 
superconductor. We will not specify the origin of this interaction: for our 
purposes it is enough to assume that this interaction is attractive and leads 
to the appearance of superconductivity with Cooper pairs of charge 2e at 
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temperatures below Tc- Of course, the smallness of the effects magnitude 
(necessary for the apphcabihty of the perturbative approach) is a serious 
hmitation of the proposed theory. Nevertheless the current state of HTS in- 
vestigations in some respects reminds one of the situation which occurs in the 
study of metal-insulator transitions in the 1970s. The weak localization the- 
ory did not describe consistently the Anderson transition, but was successful 
in the explanation of a set of anomalous properties of the disordered metal 
systems; besides this it gave a hint to the development of the renormalization 
group approach to the description of the metal-insulator transition. 

The Cooper channel of interelectron interaction is equivalent to taking 
into account superconducting fluctuations which are unusually strong in HTS 
[^. The reasons for this strength are the effective low dimensionally of the 
electron spectrum, the low density of charge carriers and the high values of 
critical temperature of HTS. We will show that these peculiarities lead not 
only to the increase of the magnitude of the fluctuation effects, but frequently 
change the hierarchy of the different fluctuation contributions, leading to 
the appearance of competition among them and even to the change of the 
habitual (in conventional superconductivity) sign of the overall correction. 
As a result, the fluctuation effects can manifest themselves in very unusual 
form, so that their origin cannot be identified at first glance. 

The first well-known result is that in the metallic phase one automati- 
cally has the non-equilibrium analogue of pre-formed Cooper pairs above Tc. 
Indeed, taking into account thermal fluctuations (or interelectron interaction 
in the Cooper channel) leads to the appearance of some non-zero density of 
fluctuation Cooper pairs (in contrast to pre-formed pairs with finite lifetime) 
in the superconducting layers without the establishment of the long range 
order in the system. It is important that in the 2D case, typical for HTS 
materials, the density of Cooper pairs decreases with temperature extremely 
slowly: ~ In ^^^^ . One should therefore not be surprised that precursor ef- 
fects can often be detected in the normal phase well above T^. (especially in 
underdoped samples, see section 5). 

The formation of fluctuation Cooper pairs of normal electrons above 
has an important though usually ignored consequence: the decrease of the 



density of one-electron states (DOS) at the Fermi level This circum 



stance turns out to be crucial for the understanding of the aforementioned 
effects and it will constitute the quintessence of this review. In this way 
the following phenomena can be at least qualitatively (and in many cases 
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quantitatively as well) explained: 

The behavior of the c-axis resistance which has been explained in 
terms of the suppression of the one-electron DOS at the Fermi level and the 
competition of this effect with the positive Aslamazov-Larkin (AL) paracon- 



ductivity l2|, |g, P |2|, |3 



The giant growth of the c-axis resistance peak in the presence of an ex- 
ternal magnetic field applied along the c-axis is explained using the same 
approach which was shown to fit well the experiments ^ ^ . 



The anomalous negative magnetoresistance observed above Tc in BSCCO 
samples ^ ^ was again explained by the same DOS fluctuation contri- 
bution Moreover, its competition with the positive Aslamazov-Larkin 
magnetoresitance gave good grounds for the prediction of a sign change in 
the magnetoresistance as temperature decreases towards Tc [|^. The latter 
effect was very recently confirmed experimentally on YBCO samples [0. 

The decrease of the thermoelectric power at the edge of transition turns 
out to be the result of the DOS fluctuation contribution which dominates over 



the AL term |^ previously assumed to play the crucial role 

The temperature dependence of the NMR rate can be explained as 
the result of the competition between the positive Maki-Thompson (MT) 
correction to the Korringa law and the negative DOS contribution at the 
edge of the transition |H6|, E7 . 



The observed pseudo-gap-like structure in the far infra-red optical con- 
ductivity along c-axis can also be attributed to the suppression of the one- 
electron DOS at the Fermi level. This leads to the appearance of a sizable 
negative contribution in optical conductivity, which shows up in a wide range 
of frequencies (up to ujdos ~ 'T""^); exceeding the positive AL and MT con- 
tributions in magnitude and range of manifestation ||48|| . 

We believe that the fact that even the simple approach proposed here 
was able to explain most of the anomalies of the normal state properties 
of HTS mentioned above is not accidental. It shows the importance of the 
interelectron interaction in the problem discussed and demonstrates that 
even the way of "up-grading" the traditional BCS theory to include the HTS 
peculiarities is creative in the explanation of the HTS properties. Further, the 
approach considered provides clear results which can be compared with those 
obtained from the alternative viewpoints and an attempt to match them in 
the region of the intermediate strengths of interaction may be undertaken. 

The review is organized in the following way. The sections 2 and 3 in- 
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troduce the reader to the short story and simplest notions of the fluctuation 
theory. The existence of a finite non-equihbrium concentration of Cooper 
pairs at temperatures above the critical one is shown and the effect of fluctu- 
ations on the order parameter and critical temperature is discussed. In the 
section 4, the renormalization of the one-electron density of states is con- 
sidered, and in section 5, its consequences on the tunneling properties arc 
discussed. This preliminary introduction prepares the reader for the central 
section 6 where the different fluctuation effects are discussed in their variety, 
first at a qualitative level and, then, within the microscopic approach. The 
results obtained in section 6 are then applied in section 7 to the analysis of 
the experimental data for the c-axis electrical transport. Section 8 is devoted 
to the effect of an external magnetic field on the c-axis transport. In section 
9 we demonstrate that the effect of DOS fluctuation renormalization causes 
the opening of a pseudo-gap type structure in the c-axis optical conductivity. 
The importance of the DOS fluctuation contribution for the thermoelectricity 
is shown in section 10. The last two sections are devoted to the discussion of 
the effect of fluctuations on the NMR characteristics at the edge of transition 
and the possibility of apphcation of the effect discussed as a tool to study 
the order parameter symmetry of HTS. 
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2 Excursus to superconducting fluctuation 
theory 

During the first half of the century, after the discovery of superconductivity 
by Kammerhngh-Onnes, the problem of fluctuations smearing the supercon- 
ducting transition was not even considered. In bulk samples of traditional 
superconductors the critical temperature Tc sharply divides the supercon- 
ducting and the normal phases. It is worth mentioning that such behavior of 
the physical characteristics of superconductors is in perfect agreement both 
with the Ginzburg-Landau phenomenological theory (1950) and the BCS mi- 
croscopic theory of superconductivity (1957). However, at the same time, it 
was well known that thermodynamic fluctuations can cause strong smearing 
of other second -order phase transitions, such as the A-point in liquid helium. 

As already mentioned, the characteristics of high temperature and organic 
superconductors, low dimensional and amorphous superconducting systems 
studied today, differ strongly from those of the traditional superconductors 
discussed in textbooks. The transitions turn out to be much more smeared 
here. The appearance of thermodynamically nonequilibrium Cooper pairs 
(superconducting fluctuations) above critical temperature leads to precursor 
effects of the superconducting phase occurring while the system is still in 
the normal phase, often far enough from T^. The conductivity, the heat 
capacity, the diamagnetic susceptibility, the sound attenuation etc. may 
increase considerably in the vicinity of the transition temperature. 

So, what is the principal difference between conventional and unconven- 
tional superconductors with respect to fluctuation phenomena, and in gen- 
eral, what determines the role and the strength of fluctuations in the vicinity 
of the superconducting transition? How smeared out is the transition point in 
existing superconducting devices, and how can one separate the fluctuation 
contributions from the normal state ones? Which microscopic information 
can be extracted from the analysis of the fluctuation corrections in different 
physical characteristics of superconductors? 

These questions, along with many others, find their answers in the theory 
of fluctuation phenomena in superconductors. This chapter of the supercon- 
ductivity has been developed in the last 30 years. 

The first numerical estimation of the fiuctuation contribution to the heat 
capacity of superconductors in the vicinity of Tc was done by Ginzburg in 
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1960 In that paper he showed that superconducting fluctuations in- 

crease the heat capacity even above Tc. In this way the fluctuations smear 
the jump in the heat capacity which, in accordance with the phenomeno- 
logical Ginzburg-Landau theory of second order phase transitions (see for 
instance [^), takes place at the transition point itself. The range of tem- 
peratures where the fluctuation correction to the heat capacity of a bulk 
clean conventional superconductor is relevant was estimated by Ginzburg as 




where a is the interatomic distance, Ep is the Fermi energy and ^ is the super- 
conductor coherence length at zero temperature []. It is easy to see that this 
is many orders of magnitude smaller than the temperature range accessible 
in real experiments. This is why fluctuation phenomena in superconductors 
were considered experimentally inaccessible for a long time. 

In the 1950s and 60s the formulation of the microscopic theory of su- 
perconductivity, the theory of type-II superconductors and the search for 
high-Tc superconductivity attracted the attention of researchers to dirty sys- 
tems, and the properties of superconducting films and filaments began to 
be studied. In 1968, in the well known paper of L. G. Aslamazov and A. 
I. Larkin the consistent microscopic theory of fluctuations in the nor- 



mal phase of a superconductor in the vicinity of the critical temperature was 



formulated. This microscopic approach confirmed Ginzburg's evaluation [49 



for the width of the fluctuation region in a bulk clean superconductor, but 



much more interesting results were found in for low dimensional or dirty 
superconducting systems. The exponent u of the ratio (a/^o), which enters 
in (|l]), drastically decreases as the effective dimensionality of the electron 
motion diminishes: z/ = 4 for 3D, but u = 1 for 2D electron spectrum (in the 
clean case) which is the most appropriate for HTS materials. 

Another source of the effective increase of the strength of fluctuation 
effects is the decrease of the coherence length, which occurs in dirty super- 
conductors because of the diffusive character of the electronic motion. This 



"'^The same expression for the width of the strong fluctuations region was obtained by 
Levanyuk in [5l[| . So in the modern theory of phase transitions the relative width of 
fluctuation region is called the Ginzburg-Levanyuk parameter G«(£)) and its value strongly 
depends on the space dimensionality D and on the impurity concentration |5Q] . 
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means that fluctuation phenomena are mainly observable in amorphous ma- 
terials with removed dimensionality, such as films and whiskers, where both 
factors mentioned above come into play. HTS is of special interest in this 
sense, because their electronic spectrum is extremely anisotropic and their 
coherence length is very small. As a result the temperature range in which 
the fluctuations are important in HTS may reach tens of degrees. 

The manifestation of superconducting fluctuations above critical temper- 
ature may be conveniently demonstrated considering the case of electrical 
conductivity. In the first approximation there are three different effects. The 
first one, a direct contribution, consists in the appearance of nonequilibrium 
Cooper pairs with the characteristic lifetime tgl = 7^h/8kB{T — T^) in the 
vicinity of the transition. In spite of their finite lifetime, a non-zero number 
of such pairs (which depends on the proximity to T^) is always present in the 
normal phase (below they are in excess in comparison with the equilib- 
rium value (see the next section)). Their presence gives rise, for instance, to 
the appearance of the precursor of the Meissner-Ochsenfeld anomalous dia- 



magnetism in the normal phase: that is ||53| the anomalous increase of the 
diamagnetic susceptibility at the edge of the transition. As far as the con- 
ductivity is concerned, one can say that above Tc, because of the presence of 
nonequilibrium Cooper pairs, a new, non-dissipative, channel of charge trans- 
fer has been opened. This direct fluctuation contribution to the conductivity 



is called paraconductivity or the Aslamazov-Larkin (AL) contribution |]52[ 



Another consequence of the appearance of fluctuating Cooper pairs above 
Tc is the decrease of the one-electron density of states at the Fermi level. In- 
deed, if some electrons are involved in the pairing they can not simultaneously 
participate in charge transfer and heat capacity as one-particle excitations. 
Nevertheless, the total number of the electronic states can not be changed 
by the Cooper interaction, and only a redistribution of the levels along the 
energy axis is possible [27, (see Fig. 2 in section 4). In this sense one can 
speak about the opening of a fluctuation pseudo-gap at the Fermi level. 

The decrease of the one-electron density of states at the Fermi level leads 
to the reduction of the normal state conductivity. This, indirect, fluctuation 
correction to the conductivity is called the density of states (DOS) contribu- 
tion and it appears side by side with the paraconductivity. It has an opposite 
(negative) sign and turns out to be much less singular in (T — Tc)~^ in com- 
parison with the AL contribution, so that in the vicinity of it was usually 



omitted. However, in many cases p9|, |33, p3, E5|, p6[, when for some special 
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reasons the main, most singular, corrections are suppressed, the DOS cor- 
rection becomes of major importance. Such a situation takes place in many 
cases of actual interest (quasiparticle current in tunnel structures, c-axis 
transport in strongly anisotropic high temperature superconductors, NMR 
relaxation rate, thermoelectric power). In this context the study of the DOS 
contribution will be our main goal in this review. 

The third, purely quantum, fluctuation contribution is generated by the 
coherent scattering of the electrons forming a Cooper pair on the same elastic 
impurities. This is the so called anomalous Maki- Thompson (MT) contribu- 
tion |5^ which often turns out to be important in conductivity and other 



transport phenomena at the edge of the transition. Its temperature singular- 
ity near is similar to that of paraconductivity, but this contribution turns 
out to be extremely sensitive to electron phase-breaking processes. In HTS 
materials there are several sources of strong pair-breaking above T^. (such as 
localized magnetic moments, thermal phonons etc.). So the MT contribu- 
tion turns out to be depressed by these phase-breaking processes and can 
usually be omitted in HTS fluctuation analysis. Nevertheless in some special 
cases (like NMR relaxation rate) it has to be taken into account even in its 
overdamped form. 

Finally, for completeness, we have to mention the regular part of the Maki- 
Thompson diagram which is much less singular and has an origin similar to 
the DOS renormalization contribution. 

We will see below that the strong anisotropy of the electron motion in HTS 
makes the DOS contribution particularly important for c-axis transport and 
related phenomena, changing significantly the hierarchy of the fluctuation 
corrections in comparison with the conventional case. 
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3 The effect of fluctuations on the order pa- 
rameter and critical temperature 

3.1 Introduction 

The order parameter temperature dependence and the critical temperature 
are among the main characteristics of superconductors. The traditional BCS 
approach gives simple expressions for both of them, but it turns out to be 
valid for 3D systems only. The effect of fluctuations in the 2D case becomes 
crucial and leads to the break down of the fundamental idea of BCS the- 



ory: the association of superconductivity with long-range order in the system. 



The works of Berezinski , Thouless and Kosterlitz |61[ demonstrated that 
the requirement (\E'(0)\E'(r))|f|^oo = const for an homogeneous 2D system is 
too rigid, actually the flow of supercurrent takes place even if the system pos- 
sesses some "stiffness" only. More precisely the aforementioned correlator of 
the order parameter at different points only has to decrease with distance 
|r| as some power law (in contrast to its exponential decrease in the normal 
metal). 

We will not discuss here the well known properties of the Berezinski- 
Thouless-Kosterlitz state (see, for instance, the review |^ |6^ ) but will 
concentrate on the crucial, for our purposes, fact of the nonequilibrium fluc- 
tuation Cooper pair formation above in quasi-two-dimensional systems. 
As we will show below, in spite of the exponential decrease of spatial super- 
conducting correlations above Tc, the density of nonequilibrium pairs in this 
case decreases only logarithmically with temperature. 

3.2 Fluctuation Cooper pairs above Tc 

Let us start from the calculation of the density of fluctuation Cooper pairs 
in the normal phase of a superconductor. We restrict ourselves to the re- 
gion of temperatures near the critical temperature, so we can operate in the 
framework of the Landau theory of phase transitions |Q . 



When we consider the system above the transition temperature, the order 
parameter \E'(r) has a fluctuating origin (its mean value is equal to zero) and 
it depends on the space variables even in the absence of magnetic field. This is 
why we have to take into account the gradient term in the Ginzburg-Landau 



15 



functional for the fluctuation part of the thermodynamical potential ^(fi)'- 



(2) 



where e = In (T/T^ 



T-Tc 



<^ 1 for the temperature region discussed and 
The positive constant rjo of the phenomenological Ginzburg- 
Landau may be expressed in terms of microscopic characteristics of the metal: 



a 



1 



D 



1 1 
^ ' 2 ^ A-ktT 



1 



1 



8DT I 7C(3)/(27r3rT) 



^ \2J AtttT 
for tT « 1 
for tT » 1 



(3) 



where vp is the Fermi velocity r is the quasiparticle scattering time, D is 
the space dimensionality, ipi^) ci-nd ip'{x) are the digamma function and its 
derivative respectively, and ({x) is the Riemann zeta function 0. Dealing 
mostly with 2D case we will often use this definition omitting the subscript 
"2": r]2 = v- 

Dealing with the region of temperatures £ > 0, in the first approximation 
we can neglect the fourth order term in (H). Then, carrying out the Fourier 
transformation of the order parameter 



(4) 



one can easily write the fluctuation part of the thermodynamic potential as 
a sum over Fourier components of the order parameter: 



^1 



(5) 



Here 



- 27r - 27r - 
k = —n^t + —riyj 



2n 



where L^^y^z are the sample dimensions in appropriate directions; are 
unit vectors along the integer numbers; V is the volume of 

the sample. 

^We will mostly use the system h = c = ks = ^ everywhere, excluding the situations 
where the direct comparison with experiments is necessary. 
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In the vicinity of the transition the order parameter \1/ undergoes equi- 
hbrium fluctuations. The probabihty of the fluctuation reahzation of a given 
configuration \l/(r) is proportional to 



V oc exp 



a 



(6) 



Hence the average equilibrium fluctuation of the square of the order param- 

,(/0 2 



eter Fourier component 



may be calculated as 



{/O 



/ 




2 

exp 






2" 


d 


2 


/exp 






2" 


d 




2 



(7) 



T 



The concentration of Cooper pairs Mc.p. is determined by the average value 



of the square of the order parameter modulus |5^. For the two-dimensional 
case, which is of most interest to us, one finds: 



= - / = 2AAf ^^In - (8) 

where A/'J^^ = -^Ep is the one-electron concentration in 2D case, and r]2 is 
defined by the expression 0. 

■^One can notice that the number of Cooper pairs in (||) surprisingly does not depend 
on the concentration of electrons even when ad ahsurdum the number of electrons tends 
to zero. In this relation it is worth to mention two circumstances. The first one consists 
in the fact that we used the degenerate Fermi gas model from the very beginning, so the 
electron density cannot tend to zero. The second comment concerns the very special role 
of the dimensionality "2" for electron systems ( disordered, or superconducting). The well 
known universality of the paraconductivity expression in this case is directly related with 
the discussed property of ( 0) (see section 6). 
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We see that in the 2D case the density of fluctuation Cooper pairs de- 
creases very slowly as the temperature increases : logarithmically only. Of 
course these are nonequilibrium pairs, their lifetime being determined by the 
Ginzburg-Landau time tcl = s{t-t ) there is no long range order in the 
system. Nevertheless, one can see that even in the normal phase of a super- 
conductor at each moment there is a non-zero density of such pairs which 
may participate in charge transfer, anomalous diamagnetism, heat capacity 
increase near transition. In this sense we can speak about the existence of 
the average modulus of the order parameter (which is defined as the square 
root of the average square of modulus (H)). 

The participation of normal electrons in nonequilibrium Cooper pairing 
above is an inelastic process leading to some decay of the phase coherence 
between initial and final quasiparticle states. This means that fluctuations 
themselves act as a source of some phase-breaking time T^{e) side by side 
with paramagnetic impurities and thermal phonons. The consequence of 
this fact is the shift of the transition temperature toward lower temperatures 
with respect to its mean field value Tcq. 

This shift is easy to estimate by taking into account the next order cor- 
rection (~ in the Ginzburg-Landau functional. We make the Hartree 
approximation by replacing the |^'(r)|^ term in (H) by (|\E'(ji)p)|\E'(r)p. This 
leads to the renormalization of the reduced critical temperature value 

^* = ^ + ^(l%0l') (9) 

and to appropriate reduction of the critical temperature T* with respect to 
its BCS value T^q. Using the microscopic values of a and b and also the 
results (0)-® for ^'^^ easily find within logarithmic accuracy: 



' c 



J-c J-cO ^i{2) 

where ^^[2) ~ ]Fid "^i^ty) ^^2? ~ p]^^ clean film of thickness d (/ is 

the electron mean free path). In the case of HTS single crystal Gi^^^^ ~ 

The consistent description of fluctuations in superconductors above Tc is 
possible only in the framework of the microscopic approach based on BCS 
theory and it will be presented later (section 6). 
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3.3 Fluctuations below T, 



c 



The description of fluctuations is much more comphcated at temperatures 
below Tc- In contrast to the simple fluctuation picture of the normal phase, 
where nonequilibrium pairs appear and decay, below the consequent con- 
struction of the fluctuation theory requires one to go beyond the simple BCS 
picture of superconductivity and to take into account correctly the different 
channels of interelectron interaction. 

The electron-electron interaction in a normal metal may be considered as 
the sum of the dynamically screened Coulomb interaction and virtual phonon 
exchange. Generally speaking both of them involve momentum transfers in 
the range of \Ap\ < 2pp. However, in normal phase, it is usually possible to 
consider the interaction as taking place in just two channels: low-momentum 
and high-momentum transfer. 

The first one is appropriate for the momentum transfers Ap and is 
called the diffusion channel or the dynamically screened Coulomb interac- 
tion. Virtual phonon exchange is neglected in it. The role of this interaction 
becomes really pronounced in amorphous systems where the delay of the 
screening, due to the diffusive character of the electron motion, is consid- 
erable. It leads to the renormalization of the density of states and other 



thermodynamical and transport properties of a metal (see [^] and section 
5). 

The second, so-called Cooper, channel is appropriate for the values of the 
momentum transfers Ap 2pp, while the total momentum of two particles 
Pi + P2 ^ 0. Here both Coulomb and virtual phonon exchange interactions 
are important and it is this part of interaction which leads to superconduc- 
tivity. In the BCS model it is described by the effective constant of the 
electron-electron interaction only. 

Above Tc these two channels of interaction do not mix, so the interaction 
corrections from both channels can be considered separately |6^. Account- 



ing for interaction in the diffusion channel leads to corrections important 
in the description of the disordered systems properties. An effective attrac- 
tion in the Cooper channel, on the other hand, leads to the reconstruction 
of the electronic ground state below to form the superconducting state. 
Taking into account the Cooper channel above is equivalent to treating 
superconducting fluctuations. 

Below Tc the interelectron interaction, in the framework of the BCS the- 
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ory, shows up itself via Tc only. Fluctuations are not taken into account here 
and the equilibrium coexistence of Cooper pairs condensate with one-particle 
excitations is supposed. Fluctuations may be taken into account by devel- 
oping the perturbation series in the interaction beyond the BCS picture. It 
is important to stress that the interaction in a superconductor cannot be 
reduced to separate diffusion (screened Coulomb interaction) and Cooper 
channels only, as was done above T^. the presence of the condensate mixes 
up the channels leading to a variety of collective processes. These include: 
the appearance and decay of non-equilibrium Cooper pairs, different types 
of quasiparticle scattering processes involving the condensate, and scattering 



of Cooper pairs between themselves |6^, pi 

Another way of considering fluctuation phenomena in the superconduct- 
ing phase is to speak in terms of fluctuations of the modulus and phase of the 
order parameter, to calculate their correlators in different points of the super- 
conductor and to study the physically measurable values |]69|, ITOl, I7l|. In this 



scheme, the inseparability of scalar potential (dynamically screened Coulomb 
interaction) and order parameter phase fluctuations below Tc becomes evi- 
dent. Indeed, above Tc the dynamical screening of the charge fluctuation 
originates by the space and time redistribution of electrons only. Below the 
critical temperature such quasiparticle currents cause charge redistribution 
by means of supercurrent flows too, i.e. the appearance of gradients of the 
phase of order parameter. One can see that this means the linking of scalar 
potential and phase fluctuations, i.e. the appearance of off-diagonal elements 
in the matrix correlation function. 

We can discuss the effect of fluctuations on Tc and the modulus of the 
order parameter \E'(T) starting from the superconducting phase (in contrast 
to the previous section). Evidently fluctuations suppress both of these with 
respect to their BCS values and, as we will see below, this effect is very 
similar to the effect of paramagnetic impurities. However, while in the case 
of paramagnetic impurities, it is possible, at least in principle, to clean up 
the sample and to determine the value of Tco, it is impossible to do this 
with fluctuations. They reduce the values of the critical temperature and 
the order parameter modulus with respect to their BCS values, but there is 
no way to "switch them off" at finite temperatures and especially near the 
critical temperature Q. 

^The effect of fluctuations in principle can be suppressed by external magnetic field but 
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We confine ourselves here to the discussion of the most interesting case of 
2D fluctuations (films with thickness d <^ (T) or strongly anisotropic layered 
superconductors). The renormalization of the order parameter modulus may 
be calculated directly from the order parameter self-consistent equation by 
including the fluctuation corrections to Gorkov F-function. In the vicinity of 
Tc one finds for a thin film |^ : 



|vI/(^,)|^(T,T,o)) = vl/^(T,T,o)- 

(11) 

^^^VJln^ + 21nA + J^l,3^Tl 



ApFdyEpJ \ Tco-T Lt 63C(3) d 

where Lt = \J^: is the diffusion length and D is diffusion coefficient. The 
inclusion of Tco amongst the arguments of and \l/o is done on purpose 
to underline that this parameter can vary due to the effect of fluctuations 
too. For a layered superconductor the analogous expression can be found by 
replacing the film thickness d by the interlayer spacing s in the coefficient, 
and replacing d — > max{s,^x ~ "5?} in the argument of the last logarithm 
(here J is a hopping integral describing the Josephson interaction between 
layers, see section 6). 

The first correction to {'^\fi)) is due to fluctuations of the modulus of 
the order parameter and is primarily responsible for the temperature depen- 
dence. One can see that this term is analogous to (\l/(j;-)(T)) from above the 
transition. 

The second correction comes from the fluctuations of the phase of order 
parameter and depends on some longitudinal cut-off parameter L which has 
to be chosen in accordance with the problem considered. Its origin is con- 
nected with the destruction of long range order by phase fluctuations in low 
dimensional systems [^] 0. Thus in the framework of such an approach the 

this will change Tc with respect to the BCS value itself. 

'"^Let us mention that whilst above Tc both degrees of freedom (modulus and phase 
of the order parameter) fluctuate in the same way, below the critical temperature their 
behavior is quite different. Namely, as one can see from the result (|l^), the modulus 
fluctuates qualitatively in the same way as above Tc, while the phase fluctuations lead 
to the destruction of the long range order. If the electron Fermi liquid were uncharged, 
one could associate phase fluctuations with the appearance of the Goldstone mode. In 
superconductors, due to the electro-neutrality condition, the massless Goldstone boson 
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longitudinal dimension of the specimen L has to be taken as a cut-off param- 
eter in (p!TD and the divergence which appears there indicates the absence 
of long-range order for infinite 2D superconductors. However, as mentioned 
above, a supercurrent can flow in the system even when the average value 
of the order parameter is not well-defined. The only requirement is that the 
correlation function of the order parameter behaves as a power of r in the 
long-range limit r — > cxd. In the calculation of thermodynamical functions 
in such a state, the presence of vortices in 2D superconductors gives rise 
to another cut-off parameter for the contribution of phase fluctuations: the 
characteristic distance between vortices. 

The last term on the right-hand side of ( pT]) represents the contribution 
both of the scalar potential fluctuations and off-diagonal phase - scalar po- 



tential interference terms |6^, |70[. One can see that it becomes important 
for very thin films with d <^ Lt, and for strongly anisotropic layered super- 
conductors. 

Such considerations lead to the renormalization of the critical tempera- 
ture which can be obtained from (|Tl]) |^9|, |70|, 0. For not too thin (but with 
d <^ ^{T)) superconducting films: 

(12) 

in accordance with (p!oD . 

In the important for HTS case of a clean strongly anisotropic layered 
superconductor the consideration of the first two terms in ([Tl|) is analogous 
but the last term In^ — ~ lii^ — tWtj generally speaking, cannot be 
omitted. As a result the shift of the critical temperature due to the inter- 
electron interaction (which includes fluctuations of the modulus and phase 
of the order parameter side by side with fluctuations of the scalar potential) 
in this case is determined by the formulae 



6T, n 



Ef 4:71^ , , Ef 



In — H — - In-^ 



Tc 63C(3) max{ J, Tc} 



(13) 



and the last term can even dominate in the case of extreme anisotropy. 



cannot propagate, but some traces of this phenomenon can be observed at finite frequencies 
l73,|7l. 
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3.4 Discussion 



We discuss here the full picture of the temperature dependence of the order 



parameter modulus y (|\E'(/i)p(T)) renormalized by fluctuations for quasi- 
two-dimensional systems (see Fig. 1): 

We start from the unperturbed BCS curve '^o{T,Tcq) (dashed-dot line 
at the Fig. 1). The effect of fluctuations is described by equation ([Tl|) and 



results in the deviation of -y/(|\E'(/i)P(T)) below '^q{T,Tco) (solid line) due 
to the growth of order parameter modulus fluctuations with the increase of 
temperature. 

The second effect of fluctuations is the decrease of the critical temperature 
with respect to Tco (see (^2])), so we have to terminate our consideration based 



on ( pTD at the renormalized transition temperature T*. One can see that the 



value of (|\E'(/i)p(T*, Tco)) is of the order of NeGi(^2)^^ and it matches 
perfectly with the logarithmic tail calculated for the temperatures above Tc 



in section 3.1. The full curve y (|^(j/)P(T)) is presented in Fig. 1 by the 
solid line. 

One further comment should be made at this point. In practice the 
temperature T^o is a formal value only, T* is measured by experiments. So 
instead of "^oiTjTco) the curve \E'o(T, T*) more naturally has to be plotted 
(dashed line at Fig. 1). It starts at T* and, in accordance with [^], finishes 
at zero temperature a little bit below the BCS value \E'o(0,Tco). The shift 
5\I'o(0) = \E'o(0, Tco) — ^^0(0, T*) due to quantum fluctuations turns out to be 
have the same relative magnitude as ST^. 

^ ' ^ — ^ Gdn (14) 
^(0) T,o ^ ^ 



One can see that the renormalized by fluctuations curve y (|\E'(//)p(T)) (pJ]) 
passes above to this, more natural parametrization of the BCS \I'o (T,T*) 
temperature dependence. 

The renormalization of the critical temperature by fluctuations (p!2D was 
first obtained by Yu.N. Ovchinnikov in 1973 by means of the expan- 



sion of the Eilenberger equations in the fluctuations of the order parameter 
and scalar potential degrees of freedom. Then it was reproduced in different 
approaches in following studies |6^, 0. Nevertheless, the studies of the ef- 
fect of fluctuations on the superconducting transition properties are still in 
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progress. For instance, the authors of |76| have revised the BCS theory in- 
cluding Goldstone bosons fluctuations. They have stated that this procedure 
led them to find the increase of the critical temperature due to fluctuations. 
Other authors |7^, ^ (including the authors of this review), are much 
less optimistic: it seems today evident that fluctuations can only reduce the 
critical temperature, and even suppress the zero-temperature value of the or- 



der parameter [^. The growth of the critical temperature discovered in [76 



is very likely related to a lack of the correct account of the "phase - scalar 
potential" correlations which leads to the exact cancelation of the logarithmic 
contribution from boson-like Goldstone phase fluctuations (see |]70|, |69|). 

The decrease of the critical temperature is specially pronounced in the 
case of bad metals such as synthetic metals and optimally doped or shghtly 
underdoped HTS compounds |^ |7^. We will see in the following that as 
oxygen concentration decreases the Gi number grows and the role of fluctua- 
tion effects (Cooper channel interaction) becomes more and more important, 
making inapplicable the perturbation theory methods far enough from the 
optimal oxygen concentration. The first example of this process is the strong 
suppression of Tc seen from the formula (|12D. 
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4 The effect of superconducting fluctuations 
on the one-electron density of states. 

As was already mentioned, the appearance of non-equilibrium Cooper pair- 
ing above leads to the redistribution of the one-electron states around the 
Fermi level. A semi-phenomenological study of the fluctuation effects on the 
density of states of a dirty superconducting material was first carried out 
while analyzing the tunneling experiments of granular Al in the fluctuation 



regime just above Tc [|^. The second metallic electrode was in the supercon- 
ducting regime and its well developed gap gave a bias voltage around which a 
structure, associated with the superconducting fluctuations of Al, appeared. 
The measured density of states has a dip at the Fermi level ^ reaches its nor- 
mal value at some energy Eq{T), show a maximum at an energy value equal 
to several times £"0, finally decreasing towards its normal value at higher 
energies (Fig. 2). The characteristic energy Eq was found to be of the order 
of the inverse of the Ginzburg-Landau relaxation time tgl ^ t- = Tce~^ 

1 I c 

introduced above. 

The presence of the depression at = and of the peak at ~ (1/ tgl) 
in the density of states above Tc are the precursor effects of the appearance of 
the superconducting gap in the quasiparticle spectrum at temperatures below 
Tc- The microscopic calculation of the fluctuation contribution to the one- 
electron density of states is a nontrivial problem and can not be carried out 
in the framework of the phenomenological Ginzburg-Landau theory. It can 
be solved within the diagrammatic technique by calculating the fluctuation 
correction to the one-electron temperature Green function with its subse- 
quent analytical continuation to the real energies ^ . We omit here the 
details of the cumbersome calculations and present only the results obtained 
from the first order perturbation theory for fluctuations. They are valid 
near the transition temperature, in the so-called Ginzburg-Landau region, 
where the deviations from the classical behavior are small. The theoretical 
results reproduce the main features of the experimental behavior cited above. 
The strength of the depression at the Fermi level is proportional to different 
powers of tgl, depending on the effective dimensionality of the electronic 
spectrum and the character of the electron motion (diffusive or ballistic). In 
a dirty superconductor for the most important cases of dimensions D=3,2 

^Here we refer the energy E to the Fermi level, where we assume E = 0. 
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one can find the following values of the relative corrections to the density of 
states at the Fermi level ||27[| : 



TV2 



6Njf{0) 



{TcTgl 



>3/2 



^3/2 



D 



D 



(15) 



where D = ^ is the diffusion coefficient. At large energies E ^ TqI the 
density of states recovers its normal value, according to the same laws (P^) 
but with the substitution tgl — >■ E~^. 

It is interesting that in the case of the density of states fluctuation cor- 
rection the critical exponents change when moving from a dirty to a clean 
superconductor p8| : 



(0) 



(TcTgl 
(TcTgl) 



1 1/2 



D 



D 



(16) 



(the subscripts (cl) and (d) stand here for clean and for dirty cases respec- 
tively). Nevertheless, as it will be seen below, due to some specific properties 
of the corrections obtained, this difference between clean and dirty systems 
does not manifest itself in the physically observable quantities (tunneling 
current, NMR relaxation rate etc.) which are associated with the density of 
states by means of some convolutions. 

Another important respect in which the character of the density of states 
renormalization in the clean and dirty cases differs strongly is the energy 
scale at which this renormalization occurs. In the dirty case (^o ^ t^i^ 



energy turns out to be [27]: 



while in the clean one (^o ^ 



T-Tr 



~ r, 



(17) 



E 



(d) 



TriT-Tr 



(18) 



To understand this important difference one has to study the character of 



the electron motion in both cases discussed 281. We recall that the size of 
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the fluctuating Cooper pair is determined by the coherence length 



e(T)=eo(^^T:f^j (19) 

of the Ginzburg-Landau theory. The zero-temperature coherence length 
differs considerably for the clean and dirty cases: 

so,d = (21) 
To pass from the dirty to the clean case one has to make the substitution 

Br^Pfir^EEl^^. (22) 

m m mTc 

The relevant energy scale in the dirty case is the inverse of the time necessary 
for the electron to diffuse over a distance equal to the coherence length ^(T). 
This energy scale coincides with the inverse relaxation time tgl'- 

t-i = Dr' ~ rjl ~ T - n. (23) 

In the clean case, the ballistic motion of the electrons gives rise to a different 
characteristic energy scale 



t^' - VfC' - {T^r^lY" - V^c(r - TJ. (24) 

The fluctuation corrections to the density of states may be presented as a 
function of the energy and the temperature in a general form, for any dimen- 
sionality of the isotropic electron spectrum and any impurity concentration 
but the relevant expressions are very cumbersome and we restrict our- 



selves to report the 2D clean case only 



-N, 



8aTr 



T} 



(2) 



1 - 



2E 



X 



In 



(25) 



2E + {4:E^ + aTcTal 



1x1/2- 
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where A'"(2) = is the 2D density of electron states in the normal metal and 
a is some number of the order of unity. One can check that the integration 
of this expression over all positive energies gives zero. This is a consequence 
of conservation of the number of particles: the number of quasiparticles is 
determined by the number of cells in the crystal and cannot be changed by 
the interaction. So the only effect which can be produced by the interelectron 
interaction is the redistribution of energy levels near the Fermi energy. This 
statement can be written as the "sum rule" for the fluctuation correction to 
the density of states: 



This sum rule plays an important role in the understanding of the mani- 
festation of the fluctuation density of states renormalization in the observable 
phenomena. As we will see in the next section the singularity in tunneling 
current (at zero voltage), due to the density of states renormalization, turns 
out to be much weaker than that in the density of states itself (Ine instead 
of or see (|T3p-(|TBp). The same features occur in the opening of the 
pseudo-gap in the c-axis optical conductivity, in the NMR relaxation rate etc. 
These features are due to the fact that we must always form the convolution 
of the density of states with some slowly varying function: for example, a 
difference of Fermi functions in the case of the tunnel current. The sum rule 
then leads to an almost perfect cancellation of the main singularity at low 
energies. The main non-zero contribution then comes from the high energy 
region where the DOS correction has its 'tail'. 

Another important consequence of the conservation law (p6|) is the consid- 
erable increase of the characteristic energy scale of the fluctuation pseudo-gap 
opening with respect to Eq: this is cVq = tcT for tunneling and a; ~ r^^ for 
c-axis optical conductivity. 




(26) 
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5 The effect of fluctuations on the tunnel 
current 



5.1 Introduction 

One of the most currently discussed problems of the physics of high temper- 
ature superconductivity is the observation of pseudo-gap type phenomena 
in the normal state of these materials []. In this Chapter we deal with the 
"pseudo-gap" type behaviour of the tunneling characteristics of HTS ma- 
terials in the metallic phase (slightly under-, optimally or over-doped com- 
pounds, where the Fermi surface is supposed to be well developed). By this 
we mean the observation at temperatures above critical one of non-linear 
I — V characteristics usual for the superconducting phase where a real gap 
in the quasiparticle spectrum occurs. 

First of all we would like to attract the readers attention to a possible, 
relatively old, scenario of a "pseudo-gap" opening in the tunnehng resistance 
. This is caused by the fluctuation renormalization of the one-electron 
density of states in a very narrow interval of energies near the Fermi surface as 
discussed in the previous section. It turns out that, due to the conservation of 
particle number, this sharp renormalization manifests itself in the tunneling 
conductance by means of the appearance of an unexpectedly wide "pseudo- 
gap" type structure with a central minimum and two lateral maxima at a 
characteristic bias eV = ztvrT and a strong temperature dependence of its 
magnitude in the vicinity of the transition. The appearance of a characteristic 
"kink" (of the same origin as the fluctuation growth of the c-axis resistance 
^ |3^) in the temperature dependence of the zero-bias conductance 
G(y = 0, T) is also predicted at temperatures close to Tc. This theory 
was checked and confirmed long ago on the conventional superconducting 
junction Al — I — Sn and we now apply it to analysis of the very recent 
tunneling experiments on HTS materials. 

^It is worth mentioning that some confusion with the concept of "pseudo-gap" takes 
place in hterature. It is used often as synonym for the spin gap (even being far from the 
anti- ferromagnetic phase), the same definition is used in the description of a variety of HTS 
normal state anomalies observed by means of transport and photo-emission measurements, 
NMR relaxation, optical conductivity and tunneling at temperatures above Tc in all range 
of oxygen concentrations without any proof that it has the same origin in the different 
experiments. 
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We start with a discussion of the effect of the inter-electron interaction 
on the tunnehng properties of superconductors. In the following subsection 
these results are then used to explain the anomaly in the zero-bias tunneling 
conductance G{V = 0,T) of a YBaCuO/Pb junction above the YBaCuO 
transition temperature |jll], |8^, Finally, we discuss the results of tra- 
ditional electron tunneling spectroscopy |ll3| , 
interlayer tunneling spectroscopy 



STM measurements 14 and 



15 1, where pseudo-gap structures in the 



tunnel conductance of BiSrCaCuO were found in a wide range of tempera- 
tures above Tr. 



5.2 Preliminaries 

It is quite evident that the renormalization of the density of states near the 
Fermi level, even of only one of the electrodes, will lead to the appearance 
of anomalies in the voltage-current characteristics of a tunnel junction. So- 
called zero-bias anomalies, which are the increase of the differential resistance 
of a junction with amorphous electrodes at zero voltage and low tempera- 
tures, have been observed for a long time. They have been explained in 
terms of a density of states depression in an energy range of the order of 
Earn ~ arouud the Fermi level due to the electron-electron interaction in 
the diffusion channel. 

The quasiparticle current flowing through a tunnel junction may be pre- 
sented as a convolution of the densities of states with the difference of the 
electron Fermi distributions in each electrode (L and R): 



roo / E + eV E \ 

I oo V''''^ ~ ^'"''^ 2t) ^^(^)^^(^ + ^^)^^' 

where Rn is the Ohmic resistance per unit area and Nl{Q), Nji{0) are the 
densities of states at the Fermi levels in each of electrodes in the absence 
of interaction. One can see that for low temperatures and voltages the ex- 
pression in parenthesis is a sharp function of energy near the Fermi level. 
The characteristic width of it is i^kcr ~ max{T, l^} <^ Ep. Nevertheless, 
depending on the properties of densities of states functions, the convolution 
( p7D may exhibit different properties. 
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If the energy scale of the density of states correction is much larger than 
ii^ker , the expression in parenthesis in (pTj) acts as a delta-function and the 
zero-bias anomaly in the tunnel conductivity strictly reproduces the anomaly 
of the density of states around the Fermi level: 

G„(0) iV(0) ' ^ > 

where GiV) is the differential tunnel conductance and Gn(0) is the back- 
ground value of the Ohmic conductance supposed to be bias independent, 
=G(l^)-G„(0). 



This situation occurs in a junction with amorphous electrodes [^. In 
the amorphous metal, the electron-electron interaction with small momentum 
transfer (diffusion channel) is retarded and this fact leads to a considerable 
suppression of the density of states in the vicinity of the Fermi level, within 
an energy range Eam^ t^^ T -Ekcr- At zero temperature for the 2D case 
one has: 



5N2iE) = ^\niET), (29) 

where the constant A is related to the Fourier transform of the interaction 
potential. In the 3D case the correction to the density of states turns out to 
be proportional to 

In the framework of this approach Altshuler and Aronov |86, 8^ analyzed 



the experimental data obtained in studies of the tunneling resistance of Al — 
/(O2) — Au junctions and showed it to be proportional to | V]^^"^ at eV <^ T. 
The identification of the "wings" in the I — V characteristics of such junctions 
with (|29| ) was a key success of the theory of the electron-electron interaction 
in disordered metals ||67|| . 

It is worth stressing that the proportionality between the tunnel current 
and the electron DOS of the electrodes is widely accepted as an axiom, but 
generally speaking this is not always so. As one can see from the previous 
section, the opposite situation occurs in the case of the DOS renormalization 
due to the electron-electron interaction in the Cooper channel: in this case 
the DOS correction varies strongly already in the scale of -Eq <^ T ~ i^^er 
and the convolution in (|27|) with the density of states (^) has to be carried 
out without the simplifying approximations assumed to obtain (pBf). This 
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statement represents the central point of this section. We will see in the 
following that this fluctuation induced pseudo-gap like structure in the tunnel 
conductance differs drastically from the anomaly of the density of states ( 
|25|) , both in its temperature singularity near Tc and in the energy range of 
manifestation. 



5.3 The effect of fluctuations on the tunnel current 

Let us first discuss the effect of the fluctuation suppression of the density of 
states on the properties of a tunnel junction between a normal metal and a 
superconductor above Tc. The effect under discussion turns out to be most 
pronounced in the case of thin superconducting films {d <C ^{T)) and layered 
superconductors like HTS cuprates. 

We now derive the explicit expression for the fluctuation contribution to 
the differential conductance of a tunnel junction with one thin film electrode 



close to its Tc. To do this we differentiate (pTj) with respect to voltage, and 



insert the density of states correction given in P3|). This gives (see 
6Gfi{V,e) _ I dE 



G„(0) 2T J- 



J— CO -12'^^ 

cosh 



2T 



(30) 



G.(,)(4vrTr)ln(i)ReV."g-^ 



where '4){x) is the digamma function, and r is the electron's elastic scattering 
time. We have introduced the Ginzburg-Levanyuk parameter Gi(2)(47rTr), 
which characterizes the strength of fluctuations in the 2D case with arbitrary 
impurity concentration, and is given by : 



Gi[2){,x%=i^TT = (31) 
Tc\ _ for47rTr«l 



^^[V'(|)-V^(| + a;) + i^'(i)],^,^^^UW Inci^ for47rTr»l- 

It is important to emphasize several nontrivial features of the result ob- 
tained. First, the sharp decrease (e"^*-^-*) of the density of the electron states 
generated by the inter-electron interaction in the immediate vicinity of the 
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Fermi level surprisingly results in a much more moderate growth of the tunnel 
resistance at zero voltage (In i). Second, in spite of the manifestation of the 

density of states renormalization at the characteristic scales Eq^^ ^ T — 
or Eq^''^ ~ \JTc{T — Tc), the energy scale of the anomaly development in the 
I — V characteristic is much larger: eV = ttT ^ Eq (see Fig. 3). 

This departure from the habitual idea of the proportionality between 
the tunnel conductance and the so-called tunneling density of states (^) 
is a straightforward result of the convolution calculated in ( with the 
difference of Fermi-functions as a kernel. As already explained in the previous 
section the physical reason is that the presence of inter-electron interaction 
cannot create new electron states: it can only redistribute the existing states. 

In the inset of Fig. 3 the measurements of the differential resistance of the 
tunnel junction Al — I — Sn at temperatures slightly above the critical tem- 
perature of Sn electrode are presented. This experiment was accomplished 
by Belogolovski, Khachaturov and Chernyak in 1986 ||82| with the purpose of 
checking the theory proposed by Varlamov and Dorin |8l| which led to the 
result (^). The non-linear differential resistance was precisely measured at 
low voltages which permitted the observation of the fine structure of the zero- 
bias anomaly. The reader can compare the shape of the measured fluctuation 
part of the differential resistance (the inset in Fig. 3 with the theoretical pre- 
diction. It is worth mentioning that the experimentally measured positions 
of the minima are eV ~ ±3Tc, while the theoretical prediction following from 
( |30D is eV = ±7rTc. Recently similar results on an aluminium film with two 
regions of different superconducting transition temperatures were reported 



We will now consider the case of a symmetric junction between two su- 
perconducting electrodes at temperatures above Tc. In this case, evidently, 
the correction (|30|) has to be multiplied by a factor "two" because of the 
possibility of fluctuation pairing in both electrodes. Furthermore, in view of 
the extraordinarily weak (~ In {-)) temperature dependence of the first order 
correction, different types of high order corrections may manifest themselves 



'Cl- 



ou the energy scale eV ~ T — Tc or yTc{T — 

The first type of higher order correction appears in the first order of 
barrier transparency but in the second of fluctuation strength (~ Gi'^) . 
Such corrections are generated by the interaction of fluctuations through the 
barrier and they can be evaluated directly from ( |27|) applied to a symmetric 
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junction 
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T he second order correction in Gi can be written as: 



Because of the evident positive sign of the integrand, the condition ( p6D does 
not lead to the cancellation observed in the first order correction SG^^i . As 
a result this correction turns out to be small (~ GP), but it is strongly- 
singular in temperature and has the opposite sign to SG^^,^ . Apparently it 
leads to the appearance of a sharp maximum at zero voltage in G{V) with 
a characteristic width eV ~ T — in the immediate vicinity of (one can 
call this peak the hyperfine structure). To our knowledge such corrections 
were never observed in tunneling experiments. 

The second type of correction is related to the coherent tunneling of the 
Cooper pairs through barrier. These appear in the second order in trans- 
parency but in the first order of fluctuation strength Gi, so they can only be 
expected to be observed in transparent enough junctions. These corrections 
are similar to the paraconductivity and anomalous Maki- Thompson contribu- 
tions to electric conductivity. Recently, some of them have been found to play 
an important role in the case of transverse resistance of strongly anisotropic 
layered superconductors, where the conducting layers are connected in the 
Josephson way p9|, The following sections will be devoted to the 

detailed discussion of the interplay of all these contributions. 



5.4 Zero-bias anomaly studies in HTS junctions 

The tunneling study of HTS materials is a difficult task for many reasons. 
For instance, the extremely short coherence length requires mono-layer-level 
perfection at surfaces which are subject to long oxygen anneals; in general 
they do not satisfy this stringent requirement. Another problem is the intrin- 
sic, (unrelated to superconductivity), bias dependence of the normal tunnel 
conductance Gn{V) in voltage ranges as wide as ~ 100 mV, in which it is 
necessary to scan for the study of I — V characteristics of HTS junctions. 
Nevertheless in recent years high quality junctions were realized and I — V 
characteristics can be measured today through a wide temperature and volt- 
age range by means of different techniques |13, 14, 15||). 
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Much attention in the scientific hterature is devoted to the temperature 
behaviour of the zero-bias conductance. The appearance of a kink around 
Tc has been often observed [Q) and in this section we demonstrate that 
quantitative fitting of G(y = 0,e) vs temperature can be obtained if the 
fluctuation contribution to the DOS in (|27|) is taken into account. 

For the junction with one thin film (2D) superconducting electrode in the 
vicinity of we have obtained the result (^); in the case of HTS tunnel 
junction it has to be modified to take into account the quasi-two-dimensional 
spectrum. In section 6, where the microscopic theory of fluctuation conduc- 
tivity will be presented, it is demonstrated that the generalization of the 
logarithmic dependence on reduced temperature e (see (pOD) from 2D to the 
quasi-2D spectrum of the corrugated cylinder type (see (0)) is trivial: it is 
enough to replace 

'■■(D-^'K^fw?)' 

Here r is the anisotropy parameter (Lawrence-Doniach crossover temper- 
ature) which is defined precisely in the following section in terms of the 
microscopic parameters of the HTS material. Physically this dimensionless 
parameter determines the reduced temperature at which the c-axis Ginzburg- 
Landau coherence length reaches the interlayer distance: C,c{£ = r) ^ s . 

Experiments have been performed on YBaCuO/Pb planar junctions ob- 
tained by chemically etching the degraded surfaces of YBaCuO single crys- 



tals by a l%Br solution in methanol. Careful quality controls have been 
carried out on these junctions in order to assure that a pure tunneling process 
without any interaction in the barrier takes place. 

In Fig. 4 we show the theoretical fitting (full line) for the normalized fluc- 
tuation part of the tunneling conductance at zero bias, 6Gfi{0,e)/Gn{0,T = 
1407^) to formula (|30D and ( ^3|) with = 0. In the inset the complete set of 
experimental data between 30 K and 180 K is shown, as reported in [HH]. The 



junction's critical temperature and the magnitude of the Ginzburg-Levanyuk 



parameter, have been taken as fitting parameters [|90[ . 

It is worth noticing that the tunneling spectroscopy probes regions of the 
superconducting electrodes to a depth of (2 3)^ in contrast with resistive 
measurements which sense the bulk percolation length. The transition tem- 
peratures determined by the two types of measurement can be quite different 
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In our case the resistive critical temperature measured on the YBaCuO single 
crystal was Tc = 91 K, while a "junction's Tc" = 89K was obtained from the 
fitting procedure. This indicates that a slightly oxygen deficient YBaCuO 
layer in the junction area is probed by the tunneling measurements. However 
the sample turns out to be still in the proper (metallic) region of the phase 
diagram. 

The value of Lawrence-Doniach crossover reduced temperature r = 0.07 
was taken from the analysis of the crossover between 2D and 3D regimes 



in the in-plane conductivity measurements Another independent mea- 



surement of r from the analysis of the non-linear fluctuation magnetization 



5^, leads to the very similar value r = 0.057. 

The magnitude of the fluctuation correction \%lj"(^)\Gi(2) = 0.029 ~ 
Tc/Ep for the sample under discussion leads to the value of Ep ^ 0.3 eV 
which is in the lower range of the existing estimates (0.2 1.0 eV). 

The temperature range in which ( |30D satisfactorily reproduces the be- 
haviour of the zero-bias conductance, extends from the "junction Tc" up to 
110 K. 



5.5 Analysis of the pseudo-gap observation experiments 

As we have mentioned above the fluctuation renormalization of the DOS at 
the Fermi level leads to the appearance of a pseudo-gap type structure in the 
tunnel conductance at temperatures above Tc which is very similar to the 
one in superconducting phase, with the maximum position being determined 
by the temperature T instead of the superconducting gap value: eVm = vrT. 
For the HTS samples this means a scale of 20 — 40 meV, considerably larger 
than in the case of conventional superconductors. 

The observations of the pseudo-gap type anomalies of G{V,T) at tem- 
peratures above Tc obtained by a variety of experimental techniques on 
BiSrCaCuO — 2212 samples were reported very recently []T3|, [1^, |T5|. Let us 
discuss their results based on the theoretical discussion presented above. 

1. The tunneling study of BiSrCaCuO — 22121 Ph junctions (as grown 
single crystal samples) where the pseudo-gap type conductance nonlinearities 
were observed in the temperature range from = 87 K — 89 up to llO/T 
(see Fig. 5 |T^. The maximum position moves to higher voltages with the 
growth of temperature and at T = 100 K it turns out to be of the order of 
30 — 35 mV (depending on the sample). Both the magnitude of the effect 
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and the measured maximum position are in qualitative agreement with the 
theoretical predictions (eKn(100 K) = ±7iT = 27 mV). 

2. The STM study of Bi2,iSri QCaCu20s+s vacuum tunneling spectra 



1^, demonstrates the appearance of pseudo-gap type maxima in G{V,T) 
starting from temperatures far above (T* = 260 K for optimum and T* = 
180 K for overdoped samples). The energy scale characterizing the pseudo- 
gap dip was estimated by the authors to be of the order of 100 mV which 
again is consistent with the predicted maximum position eKn(100 K) = 
±nT = ±27mV. 

3. Impressive results on tunneling pseudo-gap observations are presented 
T5[| . They were carried out by interlayer tunneling spectroscopy using 



m 



very thin stacks of intrinsic Josephson junctions fabricated on the surface of 
Bi2SrCaCu20^ single crystal. The opening of the pseudo-gap was found in 
dl/dV characteristics at temperatures below 180K. The data presented can 
be fitted (Fig. 6) by means of formula (^) with Gi = 0.008 and tunneling 
critical temperature Tc = 87 (in view of the strong anisotropy of the 
BSCCO spectrum r = may be assumed) pO[] . 

One can see that this fit reproduces not only the magnitude of the effect 
and the maxima positions, but also the shape of the experimental curves 
(especially at low voltages). At higher voltages the theoretical curves, being 
background voltage independent (Ohm law), overestimate the experimental 
results. 



5.6 Discussion 

We have demonstrated that the idea of relating the pseudo-gap type phe- 
nomena observed in tunneling studies with the fluctuation renormalization 
of one-electron DOS in the immediate vicinity of the Fermi surface permits us 
to fit the experimental data available with values of microscopic parameters 
(Efjt) consistent with those obtained from independent measurements. 

Two important comments are necessary. Both of them concern the limits 
of applicability of the approach proposed. 

The first concerns the magnitude of the effect. It is clear that (pOD, be- 
ing a perturbative result, has to be small, so the parameter Gi Inl/e has to 
be restricted. The temperature range is evidently restricted by e > Gi, so 
the criterion for applicability of the theory is Gi In 1/Gi -C 1 . The analysis 
of experimental data obtained on samples with different oxygen concentra- 



37 



tion gives us grounds to believe that Gi increases as the oxygen concentra- 
tion is decreased from optimal doping. It follows that the importance of 
the electron-electron interaction increases in the underdoped region of the 
phase diagram. The qualitative extension of our perturbative results to this 
poor metal - insulator region permits us to attribute the huge growth of the 
anomalies observed to the effect of the strong e-e interaction. 

Secondly it is very important to discuss the temperature range in which 
pseudo-gap phenomena are expected to be observed following the approach 
proposed. The result (|30D is obtained in the mean field region e <^ 1, ne- 
glecting the contribution of the short wave-length fluctuations. Nevertheless, 
characteristic for 2D very slow (logarithmic) dependence of the fluctuation 
correction on e = In ^ permits us to believe that the result PD| ) can be 
qualitatively extended on wider temperature range up to several T^. The 
study of high temperature asymptotics (In ^ 1) for the e-e interaction in 
the Cooper channel demonstrates the appearance of an extremely slow 
In In jr dependence which fits In 1/e in the intermediate region and shows the 
importance of the interaction effects up to high temperatures. In such a way 
one can understand the reported gap-opening temperature T* ~ 200 — 300 K 
in the underdoped part of the phase diagram as the temperature where a no- 
ticeable concentration of short-lived (r ~ -^^) fluctuation Cooper pairs first 
manifests itself. 
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6 Theory of fluctuation conductivity in a lay- 
ered superconductor 

6.1 Introduction 

Among all unconventional properties of the high temperature superconduc- 
tors, the transport properties are the most puzzling. The transition in the 
"in-plane" resistivity is usually smeared for tens of Kelvins and then a mys- 
terious linear increase with temperature takes place. If the smearing of the 
transition finds a satisfactory explanation in the framework of the fluctua- 
tion theory, the linear increase of the resistance is still the subject of much 
speculation and discussion. 

The temperature dependence of the transverse resistivity turns out to 
be even more complicated. For not too underdoped samples at high tem- 
peratures (considerably higher than Tc) Rc{T) diminishes linearly with the 
temperature in an analogous way to the in-plane resistivity. Nevertheless 
as temperature decreases, for many samples this moderate decrease is fol- 
lowed by a precipitous growth, the resistance passing through a maximum 
and then abruptly decreasing to zero as the sample passes into the supercon- 
ducting phase. Such behavior in some degree has been observed in all high 
Tc compounds [0, |3^, 0] and even in conventional layered superconductors 



As noted by Anderson ||95|, the difference in temperature dependence 
between the transverse and in-plane resistivities is extremely difficult to ex- 
plain in the framework of conventional Fermi liquid theory. Following the 
general ideology of this review we show that such an explanation is possi- 
ble by taking into account the electron-electron interaction in the Cooper 
channel. Discussing in detail the peculiarities of the fluctuation conductivity 
tensor of a layered superconductor we will identify at least one source of the 
difference between its transverse and longitudinal components: the interplay 
of the suppression of positive paraconductivity along the c-direction by the 
square of the interlayer transparency with the growth of the normal resis- 
tance due to the fluctuation depression of the density of states at the Fermi 
level. 

In this section after a short review of the theory of electrical conduc- 
tivity in a layered metal (section 6.2) the main fluctuation contributions to 
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conductivity (AL, MT, DOS) will be examined, first qualitatively (section 
6.3) and then quantitatively within the full microscopic theory (section 6.4). 
Section 6.5 is devoted to the analysis of the role of the MT contribution as 
the precursor of the Josephson effect. The results obtained are summarized 
in section 6.6. 



6.2 Normal conductivity tensor of a layered metal. 
Electron interlayer hopping time 

The goal of this subsection is to discuss the charge transport mechanisms in 
a normal layered metal with an arbitrary impurity concentration and more 
specifically the c-axis transport. Interest in this problem of classical theory of 
metals was recently revived and the problem was discussed in various models 
|9^, ^ |9^, |100|| in connection with the study of the normal properties of 



high temperature layered superconductors. Below we present the traditional 
results of anisotropic diffusion theory which will be necessary for the following 
analysis. In addition we will demonstrate, on the basis of simple ideas of band 
motion, how the crossover between the normal Drude regime and the hopping 
regime can appear 



100 



We assume that the electronic spectrum of a layered metal has the form: 

e(p) = ^(p) + Jcosip^s) - Ef, (34) 

where E{p) = p^/(2m), p = {p,Pz), P = {Px,Py) is a two-dimensional in- 
tralayer wave-vector, and J is the effective nearest-neighbor interlayer hop- 
ping energy for quasiparticles. We note that J characterizes the width of 
the band in the c-axis direction taken in the strong-coupling approximation. 
The Fermi surface defined by ^(pi?) = is a corrugated cylinder (see Fig. 7), 
and Ep is the Fermi energy. 

We assume that electrons are scattered by an elastic random potential 
of arbitrary origin, and that the scatterers are located in the layers and the 
scattering amplitude U is isotropic. 

We start from the Einstein relation for conductivity expressed through 
the diffusion constant Da/s 

a^p = N{Q)e^Ti^p. (35) 
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The dirty case (r ^ J^^) is trivial, since the usual anisotropic diffusion 
takes place. The conductivity tensor is determined by the formula (|35|) with 
a diffusion coefficient 




where < ... >denotes angular averaging over the Fermi surface. This leads 
to a Drude conductivity with the substitution of vp by its c-axis analogue 
for the transverse component: 

J)=5A'(0)eV(;4,) (37) 

Before starting the discussion of the clean case we recall the old story of 



Zener (or Bloch) oscillations ||101|| . In the case of a normal isotropic metal it 
seems that there are no restrictions on the applicability of the Drude formula 
for large r, and at first glance one might say that the conductivity of a perfect 
crystal in the absence of any type of scattering processes is infinite. However, 
more careful consideration shows that when the work produced by the electric 
field on the acceleration of the quasiparticle is accelerated by the electric field 
to the edge of the Brillouin zone, it will be Bragg reflected. Consequently 
the electrons oscillate with an amplitude which is determined by the inverse 
value of the electric field and is huge on a microscopic scale (of the order of 
centimeters for a field ~ 1—). So no charge transfer takes place in this ideal 

cm / <=> ^ 

case and the conductivity is zero. Naturally any scattering event destroys 
such a localized state and leads to a finite conductivity. 

The same situation must occur in an ideal clean layered superconduc- 
tor. The important point here is that the band along the c-axis direction 
is extremely narrow (for HTCS it is estimated that J ~ 10 — 200K) so the 
conditions on the mean free path are much less rigid than for the normal 
Zener oscillations. So one can believe that for an ideal layered crystal the 
c-axis conductivity has to be zero. Now we consider what happens if a small 
number of impurities are introduced (but r ^ Tr^ J~^, where T is the pe- 
riod of oscillations and r is the elastic scattering time introduced above). In 
this case Zener oscillations of electrons along c-axis will be destroyed from 
time to time by the scattering events. In order to estimate the transverse 
resistance in this clean case one can use formula ( |35|) with the diffusion 
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coefficient describing the process of charge transfer discussed above. The 
hopping of an electron from one plane to the next is possible only as the 
result of scattering on an impurity which dephases the Zener oscillations. So 
the diffusion coefficient in this case may be estimated as 



D?) = ^ ^ ^ (38) 

r r 

(the time average here is taken over an interval much larger then r). The 
transverse conductivity in this case is 



a?^ ~ N{0)e^- (39) 



T 



and one can see that it really vanishes in the clean limit in accordance with the 
Zener statement. This value evidently matches with transverse component 



of (|37D in the crossover region Jr ~ 1. 

Let us try to understand the physical meaning of the result obtained. 
The diffusion in the case of a strongly anisotropic superconductor should not 
be thought of as the continuous traveling of an electron in the insulating 
space between metal layers. Even in the dirty case, the c-axis motion has the 
character of wavepacket propagation, but the impurity scattering destroying 
this state occur too often: r ^ J^^(we remind that the J determines the 
period of oscillations). In a weak enough electric field (in practice the only 
reachable in real experiment) the displacement of the wavepacket during the 
time T is of the order of many interlayer distances and all above consideration 
concerning hopping to a neighbor plane remains the same as in the case of 
full oscillation. 

We can estimate the characteristic time of anisotropic diffusion from one 
layer to a neighboring one (Thop) on the basis of the Einstein relation. The 
average displacement s in the direction of the field (c-axis) takes place in the 
time {Thop), so 



and hence 



D±r£i (40) 



'^hop ~ ' 

(41) 
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Similarly we find 



- r. (42) 

The last result seems much more natural that the previous one. In the 
clean case each scattering process leads to a hop along the c-axis, while in the 
dirty case the hopping time turns out to be proportional to r^^. Nevertheless 
replacing r in ( pSD by (^TJ) reproduces the Drude formula. 

In conclusion, the results of this subsection can be summarized as follows: 

1. In the dirty limit (Jr -C 1) both components of the layered metal 
conductivity have the Drude form: 

(T_L oc (T|| oc r. (43) 

2. In the clean case (Jr ^ 1) along the c-axis direction Zener oscilla- 
tions with infrequent dephasing take place, the conductivity (T± oc has 
a hopping character, and the components of the conductivity tensor satisfy 



the relation M, 97, 98, 99 



const. (44) 



It is worth mentioning, that formally in the extremely clean case Zener 
oscillations can take place also in the a6-planes, but this possibility may be 
excluded from consideration for any real experimental situation. 

6.3 Qualitative consideration of different fluctuation 
contributions 

Let us start with the qualitative discussion of the various fluctuation contri- 
butions to conductivity [|102|| . 

The first effect of the appearance of fluctuation Cooper pairs above tran- 
sition is the opening of a new channel for charge transfer. Cooper pairs can 
be treated as carriers with charge 2e and lifetime tgl = ^(^t-t ) • '^^^^ lifetime 
has to play the role of scattering time in the Drude formula because, as is 
evident from its derivation, carrier scattering or annihilation are essentially 
equivalent. 

Finally we should replace the electron concentration A^by the Cooper pair 
concentration N'c.p. in the Drude formula. This gives the paraconductivity 
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c.p. 



(45) 



In the 2D case TVi^) 



'c.p. 2TTd Ef e' Substituting this result in ( ^Sj ) one 
reproduces with logarithmic accuracy the well known result of microscopic 



calculations 52 : 
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(46) 



The 2D result (^) is rather surprisingly identical for clean and dirty 
cases. For electronic spectra of other dimensions this universality is lost, and 
the paraconductivity becomes dependent on the electronic mean free path 
via the r]D parameter given in (^: 



6a 



AL 




16 hde 
Txr]\/2 

16 e3/25 



three — dimensional case, 

film, thickness :d <^ C,, 
wire, cross section :S <^ 



(47) 



Nevertheless one can see that the paraconductivity critical exponent u = 
{D — 4)/2 depends only on the effective dimensionality of the electronic spec- 
trum, and not on the nature of the scattering process. 

This universality of the critical exponent obtained in the framework of 
the simple mean field theory of fluctuations turns out to be very robust. The 
recent revisions of the theory of paraconductivity for the models of supercon- 
ductivity with Eliashberg's strong coupling ||103|| , self-consistent treatment of 
the coexistence of superconductivity and localization [UM, rU3] and a deriva- 
tion of TDGL theory using stochastic differential equations ||106|| resulted in 
the same critical exponents for 6aA_L- Only the overall coefficient depends on 
the characteristics of the model. 

The correction to the normal state conductivity above the transition tem- 
perature related with the one-electron density of states renormalization can 
be reproduced in analogous way. The fact that some electrons participate in 
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fluctuation Cooper pairing means that the effective number of carriers tak- 
ing part in one-electron charge transfer diminishes leading to a decrease of 
conductivity: 



imp 



imp 



(48) 



m 



m 



which in the 2D case yields 




(T,r) In - 



(49) 



The exact diagrammatic consideration of the DOS fluctuation effect on con- 
ductivity agrees with this estimate obtained in its sign and temperature 
dependence. The impurity scattering time dependence of (|4^) is correct in 
the dirty case -C /) but in the clean case the exact calculations show a 
stronger dependence on r (r^instead of r, see the next subsection). 

Finally, we discuss the Maki-Thompson contribution. This anomalous 
term has the same singularity in e as the AL one (within logarithmic accu- 
racy), but has a purely quantum nature and does not appear in the usual 
TDGL approach at all ^ Its physical nature has remained mysterious since 
1968, when it was calculated in the diagrammatic approach by Maki. This 
contribution is related to coherent electron scattering, manifest itself only 
in transport properties, and is strongly phase sensitive. These facts suggest 
that the MT contribution should be treated in the same way as Altshuler 
and Khmelnitskii ( [Q) have treated weak localization and interaction cor- 
rections to conductivity. 

Let us consider possible types of Cooper pairing above in real space. 
The simplest one is the appearance of Cooper correlation between two elec- 
trons with momenta p and —p moving along straight lines in opposite di- 
rections, (see Fig. 8(a)). Such pairing does not have the characteristics 
mentioned above and has to be attributed to the AL process. 

Nevertheless, another, much more sophisticated pairing process can oc- 
cur: one electron with spin up and momentum p can move along some self- 
intersecting trajectory; simultaneously another electron with spin down and 

^Recently it was reported |l07| that the account of the interference between super- 
fluid and normal motions of charge carriers in TDGL scheme permits to derive the MT 
contribution 
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the opposite momentum —p can move in the opposite direction along the 
same trajectory (see Fig. 8(b)). The interaction of such a pair of electrons 
during their motion along the trajectory leads to the appearance of some 
special contribution similar to the localization and Coulomb interaction cor- 
rections to conductivity but evidently singular in the vicinity of T^. One 



can easily see that such pairing is possible only in the case of diffusive motion 
(necessary for the realization of a self- intersecting trajectory). Finally any 
phase-breaking mechanism leads to the loss of coherence and destruction of 
the Cooper correlation. So all properties of the Maki-Thompson contribution 
coincide with the properties of the process proposed. The contribution to the 
conductivity of such process must be proportional to the ratio of the num- 
ber of interfering Cooper pairs ^Msa. to the full concentration of fluctuation 
Cooper pairs. In the 2D case 

^J^s.i. r-* ivpdt ^^^^ 



A/'c.p. Jtgl (Dt) 7,^ 

where is the one-electron phase-breaking time and -^^ = -^f^ is the ap- 
propriate phase-breaking rate. 

The denominator of this integral, as in , describes the volume available 



for the diffusive electronic motion with the coefficient D during the time t, 
(Dt)'^/^. The numerator describes the volume element of the tube in which 
superconducting correlation of two electron states of opposite momenta, can 
occur. The width of the tube is determined by / (mean free path) while the 
element of arc length is vpdt (see Fig. 8). The lower limit of the integral 
is chosen so that at least one Cooper pair occurs along the trajectory. The 
upper limit reflects the fact that for times t > r<^ an electron loses its phase 
and coherent Cooper pairing of two electron states above Tc is impossible. 

The contribution of Cooper pairs generated by coherent electrons moving 
along self-intersecting trajectories is therefore 

..2..A^2,„(±)(^..S..(±). (51, 

One flnds that this result coincides with the result of microscopic calculations 
of the anomalous Maki-Thompson contribution 



-2 

(TLfV = — In — (52) 
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In a similar way the correct temperature dependence of all contributions 
in all dimensions can be obtained. 

We will now try to understand the effect of fluctuations on the transverse 
resistance of a layered superconductor in the same qualitative manner. As 
we have demonstrated above, the in-plane component of paraconductivity is 
determined by the Aslamazov-Larkin formula (^61). To modify this result for 
c-axis paraconductivity one has to take into account the hopping character 
of the electronic motion in this direction. One can easily see that, if the 
probability of one-electron interlayer hopping is Vi, then the probability of 
coherent hopping for two electrons during the virtual Cooper pair lifetime 
tgl is the conditional probability of these two events: 

V2 = Vi-{Pi-tgl). (53) 
The transverse paraconductivity may thus be estimated as 

~ V2 ■ ~ Pi^l. (54) 

We see that the temperature singularity of erf ^ turns out to be stronger than 
that in a^^ because of the hopping character of the electronic motion (the 
critical exponent "2" in the conductivity is characteristic of zero-dimensional 
band motion). However for a strongly anisotropic layered superconductor 
af^, is considerably suppressed by the square of the small probability of 
inter-plane electron hopping which enters in the prefactor. 

It is this suppression which leads to the necessity of taking into account 
the DOS contribution to the transverse conductivity. The latter is less sin- 
gular in temperature but, in contrast to paraconductivity, manifests itself 
at the first, not the second, order in the interlayer transparency. One can 
estimate it in the same way as above by multiplying the in plane result (^) 
by the one-electron hopping probability: 

Aaf^^~ -Piln i. (55) 

It is important that, in contrast to the paraconductivity, the DOS fluc- 
tuation correction to the one-electron transverse conductivity is obviously 
negative and, being proportional to the first order of Vi, can completely 
change the traditional picture of fluctuations just rounding the resistivity 
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temperature dependence around transition. Excluding temporarily from con- 
sideration the anomalous Maki-Tliompson contribution (which is strongly 
suppressed in HTS by strong pair-breaking effects |108|| ), one can say that 



the shape of the temperature dependence of the transverse resistance is de- 
termined by the competition of two contributions of the opposite sign: the 
paraconductivity, which is strongly temperature dependent but is suppressed 
by the square of the barrier transparency (~ J^) and the DOS contribution 
which has a weaker temperature dependence but depends only linearly on 
the barrier transparency (~ J^), 

afi hV^\ - k2Vi\n -. (56) 

Here ki and /c2 are coefficients which will be calculated in the framework of 
the exact microscopic theory presented in the next subsection. It is this com- 
petition which leads to the formation of a maximum in the c-axis resistivity. 

6.4 Microscopic theory of fluctuation conductivity in 
layered superconductor 

6.4.1 The model 

In a layered superconductor the zero-field resistivity is a diagonal tensor 
(pa-j. = Pa, Pyy = pb, Pzz = Pc), so the calculation of its various compo- 
nents is required. Neglecting weak in-plane anisotropy, we further assume 
that the most appropriate model for high-temperature superconductors has 
isotropic in-plane electronic motion {p^x = Pyy)- Therefore, we will evaluate 
the fiuctuation corrections to the two remaining independent components of 
the resistivity tensor, p^x and p^z- 

We begin by discussing the quasiparticle normal state energy spectrum. 
While models with several conducting layers per unit cell and with either 
intralayer or interlayer pairing have been considered [|109|| , it has recently 
been shown ||11CI|| that all of these models give rise to a Josephson pair poten- 



tial that is periodic in kz, the wave- vector component parallel to the c-axis, 
with period s, the c-axis repeat distance. While such models differ in their 
superconducting densities of states, they all give rise to qualitatively sim- 
ilar fluctuation propagators, which differ only in the precise deflnitions of 
the parameters and in the precise form of the Josephson coupling potential 
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11CI|| . Ignoring the rather unimportant differences between such models in 
the Gaussian ffuctuation regime above Tc{H), we therefore consider the sim- 
plest model of a layered superconductor, in which there is one layer per unit 



cell, with intralayer singlet s-wave pairing |112|| . These assumptions lead 
to the simple spectrum (0) and hence to a Fermi surface having the form of 
a corrugated cylinder. Experiments ||113|| confirm the assumption about the 



rough isotropy of the Fermi surface in the ab-plane for metallic part of the 
HTS phase diagram. 

Some remarks regarding the normal-state quasiparticle momentum re- 
laxation time are necessary. In the "old" layered superconductors such as 
TaS'2 (pyridine) 1/2, the materials were generally assumed to be in the dirty 
limit In the high-Tc cuprates, however, both single crystals and epitax- 

ial thin films are nominally in the "clean" limit, with l/C,ab values generally 
exceeding unity, where / and C,ab are the intralayer mean-free path and BCS 
coherence length, respectively. However, as //^af, ~ 2 — 5 for most of the 
cuprates, these materials are not extremely clean. In addition, the situation 
in the cuprates is complicated by the presence of phonons for T ~ Tc ~ lOOK, 
the nearly localized magnetic moments on the Cu^"*" sites, and by other un- 
specified inelastic processes. In the following, we assume simple elastic in- 
tralayer scattering ||112|| , keeping the impurity concentration Ui and the re- 



sulting mean-free path arbitrary with respect to ^ab Q The phase-breaking 
time is taken to be much larger than r. 

We now consider the various diagrams for the electromagnetic response 
operator Qa/si^u), ( = (2i^ + l)7rTare the Matsubara frequencies) which 
contribute to the fluctuation conductivity of layered superconductors. The 
diagrams corresponding to the first order of perturbation theory in the fluc- 
tuation amplitude are shown in Fig. 9. In this notation, the subscripts a, (3 
refer to polarization directions and thus to the conductivity tensor elements 
according to 

ao,^ = - lim —[Q^pf{uj) , (57) 

Intralayer quasiparticle scattering is included in the Born approximation, 
giving rise to a scattering lifetime r and resulting in a renormalization of the 



^As it will be clear from the following assumptions for the impurity vertex calculations 
the necessity to deal with the anisotropic spectrum (^) restricts us in this section by the 
requirement I < ^abiT) = y= 
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single quasiparticle normal state Green's function to 



G{p, UJr. 



lUJr, 



(58) 



where a)„ = + l/(2|co'„|r)]. Such renormalizations are indicated in the 
vertices of Fig. 9 by shadowing. The resulting expression for the triangle 
vertex, valid for impurity concentration rii with the restrictions mentioned 
above, was calculated in |5^ (see Appendix A): 



X{q, UJn, UJn') 



1 



Q{-UnUJn 



TiUJ. 



D„ — Cbn' ) ^ 



{m-aQ-p)f) 



(59) 



where 0(x) is the Heaviside step function, and (■ ■ ■) denotes an average over 
the Fermi surface. Performing the Fermi surface average, we find. 



1 / 



(60) 



where vp 
layers. 



IPiT'l/m is the magnitude of the Fermi velocity parallel to the 



For (|59D to be valid we have to make the assumption rDg^ <^ 1. This 
will prevent us from discussing non-local corrections in this section. We will 
explain what this means at the end of this subsection. 

With each electromagnetic field component eA^ we associate the external 
vertex eVa{p) where e is the quasiparticle electronic charge, and v^ip) = 
For longitudinal conductivity tensor elements (parallel to the layers, for which 
a = x,y), the resulting vertex is simply epa/m. For the c-axis conductivity, 
the vertex evzip) is given by [ p.ll|| 



Vz{p) 



dm 

dpz 



-Js sin(p^s). 



(61) 



Each wavy line in the diagrams represents a fluctuation propagator L{q, uj^), 
which is a chain of superconducting bubble diagrams. This object, which is 
the two-particle Green's function of the fluctuation Cooper pair, was intro- 
duced in 
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and is calculated from the Dyson equation in the ladder 



approximation (see Appendix B) 



L ^{q,uJ^) = g ^ - U{q,ujf,). 



(62) 
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Here g is the effective constant of tlie electron-electron interaction in the 
Cooper channel and the polarization operator Yi{q,uj^) consists of the corre- 
lator of two one-electron impurity Green functions: 

/d^p 
— --^A(g, ujn+t,, uj^n)G{p + g, w„+/,)G(-p, (63) 

In the absence of a magnetic field and in the vicinity of T^, the inverse of 
L(g, oj^) has the form 



L~\q, uj,) = -iV(0) + V; (- + £^ + a,) - ^ ( 
where e = In (T/Tc) and 



(64) 



<? 9 2* 



(65) 



(we remind, that for simplicity we omit the subscript "2" in 7]2 : ri2 = r]). 

We integrate over the internal momenta q and sum over the internal Mat- 
subara frequencies cu^, with momentum and energy conservation at each in- 
ternal vertex (fluctuation propagator endpoint) in the analytical expressions 
for the diagrams presented at Fig. 9. 

It is worthwhile making some comment there As it was already mentioned 
above, the assumption rDg^ << 1 is necessary in order to derive an explicit 
expression for X{q,u!n,uJn') valid for an arbitrary spectrum. This expression 
will be used for the further treatment of fluctuation conductivity of layered 
superconductor. It is therefore important to understand, which are the effec- 
tive momenta involved in the following integrations. It turns out that in the 
vicinity of Tc the convergence of integrals is determined by the fluctuation 
propagator L{q,0), which means 

"9 = 9 9 ~ £ < 1, (66) 

This does not imply any restrictions on q^ but qe// ~ |~- It follows 
that the vertex part (|59| ) does not permit us to treat the non-local limit 
/ S> ^ = ^abi^) appropriate to extremely clean systems or in the immediate 
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vicinity of Tc- Since we want to apply our result to the analysis of real HTS 
compounds with / ~ 2 — 5^afe we sacrifice the non-local limit here for the 
possibility of treating the anisotropic spectrum (133). We stress again that 



the most important arbitrary (in sense of the relation between / and ^ab) case 
is accessible to our consideration. The treatment of the non-local limit will 
be presented in section 11, where its consequences on the MT contribution 
will be discussed in detail for the example of NMR relaxation rate. 

After these necessary introductory remarks and definitions we pass to the 
microscopic calculation of the different fluctuation contributions represented 
by the diagrams of Fig. 9. 

6.4.2 Aslamazov-Larkin contribution 

We first examine the AL paraconductivity (diagram 1 of Fig. 9). The proper 
contribution to the electromagnetic response tensor has the form: 

Jc^B^^q, ujf,, uju)L{q, ^p,)L{q, ^t, + uju)Bp{q, uj^, tu^), (67) 

where / d?q = J ci^q /T^/^ dq^ is the appropriate momentum space integral 
for a layered superconductor [|1 1 1| , |114|| , and the Green functions block is 
given by 

xG{p,UJn+u)G{p,UJn)G{q - p,U^_n), (68) 

(A(g, Un, ujn') is defined by (|59|), and ujn±u = =t tUy, etc). 



In the vicinity of the leading contribution to the response Qa/3 arises 
from the fluctuation propagators in ( |B7D rather than from the frequency de- 
pendences of the vertices Ba, so we can to neglect the tu^- and cj^-dependences 
of Ba- This approximation leads to [R 



i?„(q,0,0) = -2p4A([^(p)_^(q_p)]2) 

Vp oqa 



(69) 



V j sJ'^sin{qzs) ioi a = z 
Vp \ ^pQa for a = x,y 
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Notice, that since Vz{p) is odd in pz from ([6T|), Bz{q,uj^,Uu) is proportional 
to to leading order in J. 

Using the expression in (^7p followed by analytic continuation of the ex- 
ternal Matsubara frequencies to the imaginary axis (to obtain the appropri- 
ate retarded response Q^^uj)) and integration over q^, the zero-frequency AL 
contribution to the in-plane fluctuation conductivity response was found in 
the static limit to be 

^2g2^2 „ ^2q q2 



AL 

^ XX 



(27r)2[(^q2 + ^)(^q2 + ^ + ^)]3/2 

1 ^ ^ / for ^ « r, 

16s [e{e + r)]i/2 ^ iQs \ I/5, for e » r, 



(70) 



Here 



where t'(Tc) = 4^^(0)/s^ is the usual Lawrence-Doniach anisotropy param- 
eter ||115|| characterizing the dimensional crossover from the 2D to the 3D 
regimes in the thermodynamic fluctuation behavior at T^d, and ^±(0) is the 
zero-temperature Ginzburg-Landau coherence length in the c-axis direction. 

In the same way one can evaluate the AL contribution to the c-axis fluc- 
tuation conductivity ||111| , |114| 



Tie sr^ f d q 



32 J (27r)2 [(^q2 + ^)(^q2 + ^ + ^)]3/2 
e^s f e + r/2 _ ^ ^ | (r/e)^^ for e « r, ^^2) 



32ri\ [e{e + r)]y^ J GAr] \ [r / {2e)]'^ , ioi e » r 



Note, that contrary to the case of in -plane conductivity, for a^z the crossover 
occurs from OD to 3D at T/^z). 

In the region ^±(T) << s/2 of two-dimensional fluctuation behavior, 
a^^ is smaller than the static in-plane fluctuation conductivity a^^ by the 
factor [2^j_ (T)/s]^((T^/cr^), so that the other contributions to the transverse 
fluctuation conductivity need to be considered as well. The normal state 
conductivity tensor components in this model are cr^ = N{0)e'^VpT/2 = 
EpTe'^/{2ns), and cr^/cr^ = J^s'^/vp is the square of the ratio of effective 
Fermi velocities in the parallel and perpendicular directions, respectively. 
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6.4.3 Contributions from fluctuations of the density of states 

The specific forms of the AL and, as shown below, MT contributions to 
the fluctuation conductivity, which are suppressed for small interlayer trans- 
parency, suggest that one should compare these terms with those arising 
from other, less divergent, diagrams which are of lower order in the trans- 
mittance Such diagrams are pictured in diagrams 5-10 of Fig. 9. These 
(DOS) diagrams arise from corrections to the normal quasiparticle density 
of states due to fluctuations of the normal quasiparticles into the supercon- 
ducting state. In the dirty limit, the calculation of the contributions to the 
longitudinal fluctuation conductivity a^x from such diagrams was discussed 
previously Diagrams 9 and 10 arise from averaging diagrams 5 and 6 



over impurity positions. It was shown |5§] that, for a^x-, diagrams 9 and 10 



are less temperature dependent than diagrams 5 and 6, and can therefore be 



neglected. In the dirty limit, diagrams 7 and 8 were shown ||5^ to be equal to 
— I times diagrams 5 and 6, which are evidently equal to each other. In the 
clean limit, diagrams 7 and 8 can be neglected relative to diagrams 5 and 6. 
For general impurity scattering, the ratio of these diagrams depends upon r. 
As we are interested in the results for arbitrary impurity concentration, we 
shall evaluate all these diagrams separately. Calculations show that contrary 
to the case of the AL contribution, the in-plane and out-of-plane components 
of DOS contribution differ only in the square of the ratio of effective Fermi 
velocities in the parallel and perpendicular directions. This allows us to cal- 
culate both components simultaneously. The contribution to the fluctuation 
conductivity due to diagram 5 is 

Q\ti^^y) = 2e^T^ J -^^L{q,u^)TY, j -^^v^{j))vi}{p)\^{q,Un,u:^-n) 

XG'^{p,UJn)G{q - p,UJ^^n)G{p,UJn+u), (73) 

and diagram 6 gives an identical contribution. Evaluation of the integrations 
over the internal momenta p and the summation over the internal frequencies 
UJn are straightforward. Treatment of the other internal frequencies c^;^ is less 
obvious, but in order to obtain the leading singular behavior in e << 1 of 
Q^, it suflices to set = 0. After integration over q^, we have 



X 



5+6 _ _e^smj_ 
r 1 A 



{VF/Sjf 
1 



|q|<qn,ax (27r)2 [(^ + r]q2){e + r + r/q2)] 



1/2 
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where 



{vp/sJf 

1 
1 



ri 



2(Jr 



In 



+ (£ + r + r/gLx)'/'~ 



£1/2 + + r)V2 



£V2 + (£ + r)i/2 



i2 r 



^(2 + 4;^) 



47rTr 



(74) 



(75) 



In the clean hmit, azz in ( [T^ ) reduces to that obtained in [^. In ([7^), we 
have introduced a cutoff in the integral at |q| = gmax, where r/^max ~ 1, as 
in 3T|. This cutoff arises from the g-dependence of the vertices and 

of the Green's functions, which had been neglected in comparison with the 
contribution from the propagator, and is appropriate for both the clean and 
dirty limits. 

In a similar manner, the equal contributions from diagrams 7 and 8 sum 

to 



cr 



7+8 



e s'Kr2 



877 



'{vp/sJf 

1 

vp/slf 

1 



|q|<gmax (27r)2 [(^ + r7q2)(£ + r + r/q^)]^/^ 



.ei/2 + (£ + r)i/2 



where 



^2 



27r3T 



(77) 



Comparing ( [7^ ) and (|76D, we see that in the clean limit, the main contribu- 
tions from the DOS fluctuations arise from diagrams 5 and 6. In the dirty 
limit, diagrams 7 and 8 are also important, having -1/3 the value of diagrams 
5 and 6, for both a^x and azz- Diagrams 9 and 10 are not singular in e << 1 
in the 2D regime, and can be neglected. The total DOS contributions to the 
in-plane and c-axis conductivity are therefore 



D05 

^ r XX 1 



e2 

K 

2s 



ln( 



£i/2 + (e + ^)i/2y' 



(78) 
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where 





+ 2 






7r2 











f 56C(3)/7r^ ^ 0.691, 

\ 8vr2(rT)V[7C(3)] ^ 9.384(rT)2, 



for Tt « 1, 
for Tt » 1 



(79) 



is a function of rT only. 



6.4.4 Maki-Thompson contribution 

We now consider the Maki-Thompson (MT) contribution (diagram 2 of Fig. 
9) to fluctuation conductivity. The contributions from the two other dia- 
grams of the MT type (diagrams 3 and 4 of Fig. 9) are negligible, because 
they are less singular in e. Although the MT contribution to in-plane con- 
ductivity is expected to be important in the case of low pair-breaking, experi- 
ments on high-temperature superconductors have shown that excess in-plane 
conductivity can usually be explained in terms of the fluctuation paracon- 
ductivity alone. Two possible explanations can be found for this fact. The 
first one is that pair-breaking in these materials is not weak. The second is 
connected with the possibility of rf— wave pairing which does not permit the 
anomalous Maki-Thompson process at all (see section 12 for details). 

Concerning the out-of-plane MT contribution, as was stated in ||29|, even 
when the pair lifetime is short, this contribution is proportional to for 
small J above Tld, but is less singular above Tld than the AL diagram, 
and was therefore excluded from that treatment. As we shall show in the 
following, neglecting the MT diagram is usually justified for the out-of-plane 
component of fluctuation conductivity in layered materials. Nevertheless, 
there are a variety of situations where the MT diagram may be important, 
depending upon the material parameters. Because of its dependence on t^, 
the MT contribution to the transverse conductivity can have different tem- 
perature dependencies, and its order in the interlayer transmittance can vary. 
For completeness, we consider the scattering lifetime r and the pair-breaking 
lifetime to be arbitrary, but satisfying > r. The MT contribution to 
the electromagnetic response tensor is then 
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where 



Va{p)vpiq -p)X{q, )A(g 



^2 ^ 3 J ^ P \ I:' J ' '{L-^iyi " fi—'iL — iy/ ' ""/t? "^/i— ny 

xG'(p, U)n^v)G{j)-, ^n)G{q - p, UJ^_n-v)G{q - J9, . {l 

In the vicinity of Tc, it is possible to take the static hmit of the MT diagram 
simply by setting cj^ = in (|80D . Although dynamic effects can be important 
for the longitudinal fluctuation conductivity well above T;,£), the static limit 
is correct very close to T^, as shown in ||116| , |117|| . 

In evaluating the sums over the Matsubara frequencies Un in (|81|), it is 
useful to break up the sum into two parts. In the first part, uon is in the 
domains ] — oo, —uji,[ and [0, oo[. This gives rise to the regular part of the 
MT diagram. The second ( anomalous) part of the MT diagram arises from 
the summation over ci;„in the domain ] — (-o'j.jOf. In this domain, analytic 
continuation leads to an additional diffusive pole in the integration over g, 
with a characteristic pair-breaking lifetime r^. 

We start with the MT contribution to in-plane conductivity. Since in this 
case the regular part of the MT diagram is completely similar to the DOS 
contribution |]3^, we list the result only: 



a. 



MT(reg) 



where 



2s 



-Kin 



■l/2 + (5 + r)V2 



(82) 



K 









7r2 










m 



28C(3)/7r4 R 
7rV[14C(3)] 



0.3455 
=i 0.5865 



for Tr << 1, 
for Tr >> 1 



(83) 



is another function only of tT. We note that this regular MT term is negative, 
as is the overall DOS contribution. 

For the anomalous part of in-plane MT contribution we have: 



MT(an) 



1 



(27r)3 [i/r^ + Dg2][£ + ryq2 + i(i _ cosg.s)] 
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4s(£-7^)"^' ■ ' ^^^^ 
where 

2r/ n f 1 forrr«l, 

^'^ " ^;^rr<^ ^ 8Tr^ I 7C(3)/(27r3rr) for Tr » 1. ^^^^ 

Equation (|8^) indicates that in the weak pair-breaking hmit, the MT diagram 
makes an important contribution to the longitudinal fluctuation conductivity: 
it is of the same order as the AL contribution in the 3D regime, but is 
larger than the AL contribution in the 2D regime above T^d- For finite pair- 
breaking, however, the MT contribution is greatly reduced in magnitude. 

We now consider the calculation of the MT contribution to the transverse 
conductivity. From the forms ofVz{p) and ^^(g— p)in (|8TD obtained from (|6T|), 
the bare electromagnetic vertices are proportional to sin(p^s) sm{qz — Pz)s. 
After integration over the momentum p = {p,Pz), the non- vanishing contri- 
bution is proportional to cosg^s. We take the limit Jr << 1 in evaluating 
the remaining integrals, which may then be performed exactly. The normal 
part of the MT contribution to transverse conductivity is 

,MT(reg) _ _eV7rrK r jPq cosg^g 

4 J (27r)3£: + r7q2+ ^l-cosg^s) 
e^STTK f (fq f e + T]q^ + r/2 ^^^^ 



(27r)2 V [{e + r/q2)(e + r]q^ + r)]V2 
e'^srk / (e + rV^'^ — e-^/^x 2 



e^srk J 1 for £ << r, 

~* 16T[\r/{4:e) for e » r, 

This term is smaller in magnitude than is the DOS one, and therefore makes 
a relatively small contribution to the overall fluctuation conductivity and to 
the temperature at which the c-axis resistivity is a maximum. In the 3D 
regime below Tld, it is proportional to J^, and in the 2D regime above Tld, 
it is proportional to J^. 

The anomalous part of the MT diagram gives rise to a contribution to 
the transverse conductivity of the form 

^MT(an) VTcV^sV f (fq COS q^S 



zz 



(2vr)3 [i/r^ + Dg2][e + ,^q2 + r(i _ cosg.s 
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Tie s f d^q 



4(5 - 7^) i (27r) 



'y^ + T]q^ + r/2 e + r]q^ + r/2 



In examining the limiting cases of (|87|) , it is useful to consider the cases 
of weak (7^ << r, = J^rr^ >> 1/2) and strong (7^ >> r, = J^rr^ << 1/2) 
pair-breaking separately. For weak pair-breaking, we have 



"rh^Y^"^ for e << 7,^ << r 



16?7 



r/{2e) if 7,^ << r << e. 



In this case, there is the usual 3D to 2D dimensional crossover in the anoma- 
lous MT contribution at T^d, for which e(Ti£)) = r. There is an additional 
crossover at Ti (where Tc < Ti < Tld), characterized by e(Ti) = 7,^, below 
which the anomalous MT term saturates. Below T^d, the MT contribution 
is proportional to J, but in the 2D regime above T^d, it is proportional to 

For strong pair-breaking, 

^Mr(an) _^ ^ir/'y^ for e « r « 7,^, 

32?7 \ r V(47^£:) for r << min(£:, 7,^). ^ ^ 

In this case, the 3D regime (below Tld) is not singular, and the anomalous 
MT contribution is proportional to J^, rather than J for weak pair-breaking. 
In the 2D regime, it is proportional to for strong pair-breaking, as op- 
posed to for weak pair-breaking. In addition, the overall magnitude of 
the anomalous MT contribution with strong pair-breaking is greatly reduced 
from that for weak pair-breaking. 

Let us now compare the regular and anomalous MT contributions. Since 
these contributions are opposite in sign, it is important to determine which 
will dominate. For the in-plane resistivity, the situation is straightforward: 
the anomalous part always dominates over the regular and the latter can be 
neglected. The case of c— axis resistivity requires more discussion. Since we 
expect > r, strong pair-breaking is likely in the dirty limit. When the 
pair-breaking is weak, the anomalous term is always of lower order in J than 
the regular term, so the regular term can be neglected. This is true for both 
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the clean and dirty limits. The most important regime for the regular MT 
term is the dirty limit with strong pair-breaking. In this case, when t^T ~ 1, 
the regular and anomalous terms are comparable in magnitude. In short, 
it is usually a good approximation to neglect the regular term, except in 
the dirty limit with relatively strong pair-breaking and only for out-of-plane 
conductivity. However, we include it for generality. 

As discussed in greater detail in the previous section, when Jr << 1, 
the effective interlayer tunneling rate is of the order of J^r. When « 
.Pt « 1/r, the quasiparticles scatter many times before tunneling to the 
neighboring layers ||112|| , and the pairs live long enough for them to tunnel 



coherently. When J^r << l/r^, the pairs decay before both paired quasi- 
particles tunnel. 

It is interesting to compare the anomalous MT contribution with the 
DOS and AL contributions to transverse fluctuation conductivity. This com- 
parison is best made in the 2D regime above T^^). For weak pair-breaking, 
the anomalous MT and DOS terms are proportional to J^, but opposite in 
sign, and the former has a stronger temperature dependence than the lat- 
ter. With strong pair-breaking, the anomalous MT and AL contributions 
are proportional to J^, but the former is less singular in e << 1 than is the 
latter. Hence, the transverse MT contribution is in some sense intermediate 
between the transverse DOS and AL contributions. Nevertheless, as we will 
show in the following, the MT contribution can be important in the overall 
temperature dependence of the transverse resistivity, eliminating the peak 
for weak pair-breaking. 

6.5 MT anomalous contribution as the precursor phe- 
nomenon of Josephson effect 

Now, after the cumbersome calculations, let us speculate about the nontrivial 
results obtained for the MT anomalous contribution to the c-axis current. 



Comparing the results of section 5.2 (|]8T|, |28[|) and the section 6.4 one can 



notice the richness of the approach based on spectrum (Rl) with that based on 



tunneling Hamiltonian without the momentum conservation [^. The latter 
does not take into account such delicate effects as the fluctuation pairing 
through the barrier, coherent hopping etc., while the former provides these 
possibilities. One can ask: does a precursor phenomenon of the Josephson 
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effect exist and, if so, witli ffuctuation process is it? 



I.O.Kulik replied positively on this question [118, 119|. He demonstrated 



that although the average Josephson current above in the absence of the an 
applied voltage is zero, precursor radiation of the Josephson junction can be 
expected above at the frequency ~ T — Tc. He associated this radiation 
with the AL paraconductivity. 

The analysis of the expression for contradicts to this association. In 
spite of the fact that the related current is due to the motion of Cooper pairs, 
it is evidently proportional to the square of transparency of the junction 
(~ J^) and can in no way be matched with the Josephson component of 
current flowing through the junction at temperatures below Tc. The DOS 
contribution occurs in the first order of the barrier transparency , but is 



evidently due to the quasiparticle branch of the tunnel current below Te(p9 
pip. So the last candidate is the MT contribution and it is easy to see that 
under specified conditions it really represents the precursor phenomenon of 
the Josephson effect. 

One can see that in the weak pair-breaking limit (r^ ^ thop) the af^'^ 
is proportional to the first order of the barrier transparency and so the hy- 
pothetical current associated with the MT process can be matched in this 
parameter with the Josephson current below the transition. In the opposite 
case (r</, -C thop) there are no traces of the Josephson effect above Tc. 

These results can easily be understood from the picture of self-intersecting 
trajectories proposed in section 6.3. The hopping character of the electron 
motion along the c-axis leads to the minimal self-intersecting trajectory in 
this case consisting of intralayer diffusion, followed by scattering to a neigh- 
boring layer, diffusion in this layer, and hopping back to the same point in 
the original layer (see Fig. 10). The Josephson current below Tc is related to 
the phase coherence of the Cooper pairs condensate of both electrodes. Since 
the MT contribution appears as the pairing of two electrons moving along 
the trajectory described one can see that the condition (r^ ^ thop) is evi- 
dently necessary to permit the required phase coherence along the minimal 
self-intersecting trajectory. 

And there is a final argument, formal but convincing. The Josephson 
current is represented diagrammatically as the correlator of two Gorkov's 



F-functions ( [p.20| , |69[). Approaching from below the F-functions can be 
expended in powers of \1/ (see Fig. 11). Above T^ = but the correlator 
is nothing else as the fluctuation propagator and we pass to the MT 
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diagram. 

The message of these speculations is the following: the layered super- 
conductor with low phase-breaking could, in principle, irradiate or react to 



external irradiation at a frequency ~ T — Te[ |118|| in the vicinity of T^. This 
radiation can be related both to order parameter and pancake vortex fluctu- 
ations (at T < Tc). Nevertheless, the estimate for HTS compounds are very 
pessimistic: the necessary condition S> jj^y Jr -C 1 is unrealistic. 



6.6 The total fluctuation conductivity 

From the previous considerations, the total zero-field in-plane and out-of- 
plane fluctuation conductivities are found to be 
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We note that the second term in (|9TD contains both the AL and the regular 
MT contributions. Writing the AL term in this fashion, it is easy to see that 
the AL term is generally larger in magnitude than the regular MT term, 
except when r << e ^ 1. 

Although equations (pOD and ( |9TD contain both positive (AL and anoma- 
lous MT) and negative (DOS and regular MT) contributions, the behavior 
of the overall in-plane and out-of-plane fluctuation conductivities is qualita- 
tively different. 
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In fact, for a^., the negative contributions are less than the positive ones 
in the entire temperature range above the transition, leading to a total cor- 
rection which is always positive. On the other hand, in the case of cr^], both 
positive terms (AL and anomalous MT) are suppressed by interlayer trans- 
parency, leading to a competition between positive and negative terms. This 
can lead to a maximum in the c-axis fluctuation resistivity (minimum in the 
c-axis fluctuation conductivity), whilst the in-plane resistivity is expected to 
be monotonous. Since the temperature dependencies of r, k, and k are weak 
compared with that present in e, to obtain the position of the Pc{T) max- 
imum, it suffices to extremise (^) with respect to e. Using the restriction 
Jr << 1 for the validity of our theory, it is sufficient to consider the cases 
in which the resistive maximum occurs in the 2D regime. Setting e = at 
T = Tm, we have 

1 1 r~ 1 1 

which is valid for rn << land 7^^: > 1. 

We see that the regular MT term decreases the position of the maximum 
somewhat, but the anomalous term increases it somewhat. It is then qual- 
itatively correct to neglect the MT terms altogether (as confirmed e.g. by 
[ |121|| ), but quantitatively, they can change the overall shape of the fluctua- 
tion resistivity. We remark that there may be cases in which could occur 
in the 3D regime, but such cases cannot be addressed by our theory, as they 
would require a proper treatment of non-local effects, as well as the removal 
of the restriction Jr << 1. 

In Figs. 12 and 13, we have plotted Pxx/Pxx Pzz/p^z versus T/T^o for 
various values of the scattering lifetime tT^q and the pair-breaking lifetime 
r^T,o. We have taken af, = \N{fi)e^vlT and = |Ar(0) J^e^sV here. 

In these flgures, the solid curves are plots of Pzz/ P^z- The dashed curves 
are plots of Pxx/ Pxx^ same sets of parameters. In Fig. 12, we have 

chosen rT^o = 1, which is relevant for the high-TcCuprates. In each of these 
flgures, curves for t^Tcq = 1, 10, 100 are shown. We have r(Tco) = 0.1,0.01, 
and 0.001 in Figs. 12(a), (b),(c), respectively. These values correspond roughly 
to those expected for YBCO, BSCCO, and Tl2Ba2CaCu208+5 (TBCCO), 
respectively. In order that the overall fluctuation conductivity not give a large 
correction to the normal state conductivity at temperatures well above the 
transition (i. e., for T > 1.03Tco), we have chosen Ep/Tco = 300 for r(Tco) = 
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O.land 0.001 [Figs. 12(a), (b)], but Ep/T^o = SOOfor r(Tco) = 0.001(Fig. 
12(c)). As can be seen from each of these figures, for a fixed amount of pair- 
breaking and intralayer scattering, strong pair-breaking (e. g., t^T^q = 1) 
gives rise to a peak, or maximum in Pzz/pfz- small, broad peak in Pxx/Pxx 
can also occur, but only for weak anisotropy (r(Tco) = 0.1) and for such 
strong pair-breaking (r^Tco = 1). Increasing the anisotropy (or decreasing 
r(Tco)) greatly enhances the magnitude of the peak in p^z/p^z- Decreasing 
the amount of pair-breaking decreases the amplitude of the peak, as seen in 
each figure. 

In Fig. 13, plots with the same parameters as in Fig. 12(b) are shown, 
except that the intralayer scattering lifetime has been decreased to tTco = 
0.1, which is in the dirty limit. It can be seen that the magnitude of the 
peak in Pzz/p^z is reduced by interlayer hopping, by interlayer scattering, 
and by pairbreaking. For highly anisotropic materials, no peak in pxx/ Pxx 
expected for any amount of pair-breaking shown in these figures. 

Another important issue which should be discussed is the role of the 
DOS term in the interpretation of in-plane resistivity data within fluctuation 
theory. As we have mentioned above, the DOS term in a[l. cannot result in 
a change of sign, but is able to change the magnitude of the total correction 
which is important for the quantitative comparison with experimental data. 
The DOS correction gives rise to the term proportional to k in (0). In the 
dirty limit, since k — > 0.5k, the DOS contribution is not much larger than 
the (relatively small) regular MT term, and was therefore neglected by all 
previous workers. However, when rT^o = 1, h{Tcq) = 14.3123, which is much 
larger than k(Tco) = 0.5578, and the DOS contribution cannot be neglected 
relative to the other terms. Hence, fits to data in which the DOS term has 
been neglected can only be trusted for systems which are in the dirty limit. 
Since the cuprates are thought to have tTco ~ 1, it is necessary to include 
the DOS contribution in the fits. This term dramatically alters the shape of 
the overall parallel resistivity, even for zero magnetic field. This change in 
Pxx due to the inclusion of the DOS contribution is pictured in Fig. 14. In 
this figure, we have plotted Pxx/Pxx '^'^co = 1, t^Tcq = 10, Ep/Tco = 300, 
for both r(Tco) = 0.1 (dashed curves), and r(Tco) = 0.01 (sohd curves), both 
with and without the DOS contribution. As is easily seen from Fig. 14, the 
DOS contribution greatly alters the overall resistivity, with the main aspect 
of the alteration being an overall increase in the resistivity. 

However, as we will discuss in section 8.3.1, no experiment revealed up 
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to now a significant contribuion of tlic DOS fluctuations. Indeed, this small 
contribution can be masked by a change in the parameters of fits which do 
not include the DOS term. 
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7 Experimental observations of fluctuation con- 
ductivity in HTS 

7.1 Introduction 

We have already remarked that the normal state electrical transport prop- 
erties of HTS are very peculiar in many aspects. Particularly interesting 
are the differences between the in-plane and c-axis resistivities. While the 
quantitative differences among them (up to a factor 10^) are obviously due to 
the layered nature of these compounds, the explanation of their qualitative 
behavior is related to the quasi-2D character of HTS in a much less straight- 
forward way. The apparently opposite behaviors of pahiT) and pciT) close 
to Tc in HTS (decrease of resistivity for the former, increase for the latter as 
temperature is decreased) could not be explained for a long time. Although 
the peak in pdT) observed in all high Tc compounds [|T], 0, ^, 0, |122| , |123 



appeared to be strictly connected with the anisotropy of the sample, being 
very pronounced in highly anisotropic samples (BSCCO, YBCO annealed 
in reducing atmosphere) but almost absent in samples with low anisotropy 
(fully oxygenated BSCCO or YBCO), every attempt to explain it through 
normal state conductivity models in highly anisotropic systems failed to give 
satisfactory results 124, 125 , 126 , 127 , 12S, 129 1. This failure extends to 



recent models ||130| , p7[ which, as we will see, do not satisfactorily match the 
experimental behavior in the vicinity of T^. 

However, it has been shown in section 6 how the different behaviors of 
Pab{T) and Pc{T) close to Tc can be explained by a single physical mechanism 
(namely fluctuations), provided that all fluctuation contributions are taken 
into account, the hierarchy of the various fluctuation contributions being 
different for in-plane and c-axis conductivity. As we have seen, the suppres- 
sion of the positive fluctuation paraconductivity along the c-direction by the 
square of the interlayer transparency together with the decrease of the nor- 
mal state conductivity due to the fluctuation decrease of the density of states 
at the Fermi level leads to an increase of resistivity in c-axis measurements 
in samples having a sufficiently high anisotropy. 

In the following sections we will analyze experimental data in the frame- 
work of the theory presented in section 6 and show how the theory can 
quantitatively describe the experimental data. Section 7.1 will be devoted to 
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a brief review of the role of fluctuations in the in-plane conductivity, where 
the AL paraconductivity dominates, section 7.2 will deal with the origin of 
the resistivity peak in transverse measurements, where competition between 
the positive AL paraconductivity and the negative DOS contribution to con- 
ductivity takes place. A comparison between the two methods employed to 
explain the transverse resistivity peak (normal state conductivity or fluctu- 
ation conductivity) concludes this section. 



7.2 In-plane resistance: Crossover phenomena obser- 
vations 

Soon after the discovery of HTS superconductivity in the YBCO compound, 
the observation of a large in-plane excess conductivity above the supercon- 



ducting transition in measurements was reported by Freitas et al. ||131|| and 
Dubson et al. [|13|]. This excess conductivity was attributed to thermo- 
dynamic fluctuations, which in HTS are expected to be much larger than 
in conventional superconductors because of their short coherence length and 
high transition temperature. An early review of the fluctuation effects on the 



electrical transport properties was given in ||133|| . As pointed out in ||134| , |135| 



earlier measurements, especially when carried out on bulk single crystals, 
were sometimes affected by (oxygen-content) inhomogeneities. Such dif- 
ficulties were surmounted mostly by using monocrystalline films with well 
defined geometries. A typical Pab{T) curve in HTS is shown in Fig. 15. 

These measurements generated great interest, especially because the na- 
ture of thermodynamic fluctuations is related to the important topic of the 
dimensionality of superconductivity in HTS JlS^, |l38l |139|, |140|, |141], [l42 



|143| , |144|| . In fact, since the cuprate superconductors have layered structure, 
both 2D and 3D behavior can be observed, depending on the relative values 
of the temperature-dependent coherence length perpendicular to the layers 
C,c(T) and of the spacing s between superconducting layers. When ^dT) 3> 
s (i.e. close to T^) the behavior is 3D, while as the temperature is increased 
and ^c(^) decreases, a Lawrence-Doniach crossover between the 3D and 2D 
regimes occurs at the temperature Tld defined by the condition ^c(^ld) ~ •s- 



For the YBCO compound [|131|| the excess conductivity has been described 
by the 3D AL fluctuation contribution (Aa ~ e^^^"^) in a wide temperature 
range —4 < Ine < —2, while at higher temperatures [Ine > —2) a breakdown 
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of this simple theory was observed as shown in Fig. 16. These conclusions 
were confirmed by later experiments [|137| , |138 , 139 , |140| , 141, 142 1 



On the other hand, the much more anisotropic BSCCO compound showed 
a 3D fluctuation behavior only in the close proximity of T^, while at higher 



temperatures a clear 2D behavior (Act ~ e"^) was observed ||145 |. 

At sufficiently high temperatures the experimental behavior of fluctua- 
tion conductivity deviates from the simple AL theory, which is indeed valid 



only for T — Tc <^ Tc. Reggiani et al. [|117|| generalized the 2D AL theory 



for the high temperature region by taking into account the short wavelength 
fluctuations and obtained the following universal formula for paraconductiv- 
ity 



Here f{e) is a function calculated in ||117|| which, for clean 2D superconduc 



tors, tends to /(e) = 1/e in the GL region of temperature {e -C 1), so that 
the result coincides with the well known AL one, while in the opposite case 
(e > 1 ), tends to the asymptotic f{e) ~ l/e^ = l/ln^(T/Tc). 



In 1 136 the validity of this formula was carefully verified. Clearly, when 
dealing with the small fluctuation effects measured far from Tc extreme care 
must be taken in the method of measurement. Small sample-dependent de- 
viations of the R{T) curves from the optimal, "intrinsic" behavior could 
severely affect the fluctuation effects deduced from the R{T) measurements. 
Therefore, from a batch of 20 epitaxial Bi2Sr2C aCu20s+x films, only three 
were selected after checking their compositional homogeneity, structural qual- 
ity and electrical transport properties (low extrapolated resistivity at K, 
narrow transition). The R{T) curves of these three films were directly com- 
pared by plotting their resistance normalized to a reference value (namely 
i?(1.33Tc)) vs. the normalized critical temperature T/T^ as shown in Fig. 
15. The curves for all three films superimpose very well, so that they can be 
assumed to represent the "intrinsic" resistive transition for the 2212 BSCCO 
compound in spite of the small spread of their critical temperature and re- 
sistivity values. For these films the excess conductivity was analyzed. A 
very good fitting with the formula ( p3|) was found in the temperature region 
0.02 <e< 0.14 (i.e. 1.5 < T - < 11 K), while the original AL theory 
fits the data only in a much narrower temperature range (Fig. 17). For tem- 
peratures below Tc + 1.5 Eq. p3| fails because the sample is no longer in 
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the 2D region, while for T > Tc + 11 K the usual choice of the normal state 
resistance as the linear extrapolation from the high-temperature behavior 
artificially forces the extracted fluctuation conductivity to go to zero as the 
temperature approaches the range used for the linear extrapolation. 

As already mentioned, formula (|93|) predicts the asymptotic behavior 
f{e) ~ 1/e'^ = l/ln^(T/Tc) at high enough temperatures. In a recent paper 
[0 a careful analysis of the higher temperature region (above the edge of the 
region investigated in [ p.36|| ) permitted observation of this asymptotic regime, 
although at a surprisingly low reduced temperature (e* ~ ln(T*/Tc) ~ 0.23) 
in YBCO, 2212 BSCCO and 2223 BSCCO samples. The background of the 
normal state conductivity ctn = 1/pn was evaluated with particular accuracy 
by starting the linear extrapolation at temperatures higher than about 150 
K and checking that in the range from 150 K to 330 K pn did not change by 
shifting the interpolation temperature region. Therefore the upper limit of e 
at which the excess conductivity could be analyzed was ~ ln(160/92) = 
0.55 for YBCO, e„p ^ 0.46 for 2212 BSCCO and e^p ^ 0.51 for 2223 BSCCO. 

In Fig. 18 {16hs/e^) is plotted for the three samples as a function of e; 
the solid line represents 1/e, the dashed line 1/e^ and the dotted line 3.2/ y/e. 
The value of the interlayer distance s is adjusted so that the experimental 
data follow the 1/e behavior in the temperature region where the AL behavior 
is expected. Obviously, the extension of the region where the 2D AL behavior 
(l/e) is followed depends on the sample anisotropy. The less anisotropic 
YBCO compound asymptotically tends to the 3D behavior (1/e^/^) for e < 
0.1, showing the LD crossover at e ^ 0.07; the 2223 BSCCO curve starts 
to bend for e < 0.03 while the most anisotropic 2212 BSCCO shows a 2D 
behavior in the whole temperature range investigated. All three compounds 
show a universal high temperature behavior of in-plane conductivity in the 
2D regime, above the LD crossover. At £ ~ 0.24 all the curves bend down 
and follow the same asymptotic 1/e^ behavior. Finally at the value e ~ 0.45, 
close to the values of e^^p reported above, all the curves fall down indicating 
the end of the observable fluctuation regime. 

A further refinement of the theory of the in-plane fluctuation conductivity 
was carried out considering two different interlayer distances and different 



strengths of the tunneling coupling between adjacent layers for YBCO [|146 



147|| and BSCCO 2212 ||148|| . This approach led to the same qualitative 



results as those obtained by the conventional single layer approach, but to a 
better quantitative agreement between the theory and experiments. 
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In conclusion, the analysis of in-plane fluctuation conductivity has shown 
that in the whole temperature range from Tc to temperatures high enough 
that the fluctuation contribution becomes lower than the experimental reso- 
lution, the AL theory, corrected to take into account short wavelength fluc- 
tuations, is able to explain all the experimentally observed features in HTS, 
correctly describing the crossovers between different regimes as the temper- 
ature is increased. Therefore in most cases there is no need to introduce 
other contributions to fluctuation conductivity (DOS, regular and anoma- 
lous MT), since they are small with respect to the AL one as predicted by 
(pOl). However very accurate measurements have shown that for YBCO, a 
logarithmic contribution to a[l exists at temperatures as high as 180 K which 



could originate from in-plane DOS fluctuations |149 



7.3 Out-of-plane resistance: fluctuation origin of the 
c-axis peak 

The quantitative agreement of the transverse fluctuation theory outlined in 
section 6.3 with the experimental data was proved shortly after its proposal 
3l|, |32|, by fitting the resistivity peaks of BSCCO and YBCO samples. 



This shows good metallic behavior far from the transition, and thus had 
relatively small resistivity peak (Fig. 19). For strongly oxygen deficient 
samples, the increase in the c-axis resistivity begins so far from and the 
peak has such a large magnitude [0] that it cannot be due to fluctuation 
effects only: in this case the effect is probably due to some metal-insulator 
transition. 



In the theory was fit to data from 2212 BSCCO films grown on 
misaligned substrates (to allow measurement of c-axis resistivity on epitaxial 
films) in the temperature region 93—110 K. The fit used only the DOS contri- 
bution since the temperatures are far enough from Tc for the AL contribution 
to be negligible. In this way, using a single fitting parameter representing the 
amplitude of the DOS fluctuation correction to conductivity, the agreement 
of the logarithmic increase of resistivity with experimental data was proven. 



This analysis was later completed by enlarging the fitting region to in 



elude the peak and considering the AL contribution besides the DOS one. 
Their different temperature dependencies allow them to be separated and 
therefore to extract the values of the physical parameters involved. Two fit- 
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ting parameters were necessary in this case, the Fermi velocity Vp = 1-4 • 10^ 
cm/s and the electronic elastic scattering time r ~ 5 • 10"^'* s. In the same 
paper the carrier concentration and the anisotropy of a BSCCO film were 
changed by means of reducing and oxidizing annealing treatments. As the 
AL contribution is heavily dependent on the interlayer coupling than the DOS 
one, a more pronounced peak is expected for materials with higher anisotropy. 
The carrier density also affects the magnitude of the peak, since a higher car- 
rier concentration means a lower fluctuation contribution as compared to the 
normal-state conductivity. These facts strongly reduce the relative change in 
conductivity for samples having high oxygen content. The evolution of the 
resistivity peak under redox treatments confirmed these predictions. 

We have just seen that the main features of the resistivity peak in HTS 
can be explained by the fluctuation theory including the DOS contribution. 
However, this is not the only possible approach. Indeed, both before and 
after this theory was proposed and experimentally checked, many attempts 
were carried out to describe such a peak by means of non metallic normal 
state conduction mechanisms. Although early attempts in this direction 
failed to give a satisfactory description of the steep increase of resistivity 
just above [[124| , |125| , p7| , |126| , |127| , p.28|| , recent approaches involving an 



"activated" behavior of the normal state resistivity seem more promising. It 
is therefore important to discuss these mechanisms in order to see whether 
the less traditional description in terms of fluctuations is really necessary. 
We will here analyze the resonant tunneling model proposed by Abrikosov 
[|13CI|| and a phenomenological model proposed by Yan et al. . 

According to the model of Abrikosov the ratio between c-axis ad ab-plane 
resistivity is 

Pc{T) ^ cosh(ro/r) + cosh(ri/T) 

p,,(T) sinh (Ti/T) ^ ^ 

where Tq and Ti are respectively the mean energy and the half width of the 
energy spread of the resonant defects relative to Fermi level. The agreement 
of this formula with experimental data for both YBCO ||13CI|| and BSCCO 



is very good in the whole temperature range between the temperature at 
which the ratio of the resistivities shows a maximum, and room temperature. 
Even when the oxygen concentration in the samples is varied by annealing 
over a fairly wide range the quality of the fit remains good. However this 
large-scale agreement of the theory with experiment is not conclusive: re- 
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ally, the important temperature region for the comparison of normal-state 
and fluctuation theories of the resistivity peak is very close to Tc. It is here, 
indeed, that the two approaches are fundamentally different, since the fluc- 
tuation theory predicts a weak (logarithmic) divergence of Pc{T) at T^, while 
any normal state theory can at most provide a divergence a.tT = OK. If the 
PciT) peak is fluctuation induced, therefore, a normal state theory should 
not be able to reproduce the divergence close to whilst we know that our 
fluctuation theory can. The comparison of the two approaches must therefore 
be carried out by comparing with experimental data the calculated "trends" 
of Pc(T) close to Tc, rather than the absolute values of pciT). 

A closer look at the data on to BSCCO shows indeed that the fit 
with (^) tends to slightly overestimate the experimental data at both high 
and low temperatures, while underestimating it at intermediate temperatures 
(except for the oxygen annealed sample, which has a weak peak). Appar- 
ently, the fit is forced to increase as much as possible the curvature of the 
theoretical function to try to match the steep increase in Pc{T) just above 
Tc. It is likely therefore that a good description of the region just above 
Tc (where fluctuation theory predicts a singular behavior) cannot be given 
by this model. To better understand this point in |Q the curvatures of the 
experimental data and of the fitting functions above Tm were compared. The 
result was that in spite of the apparently good fit of Pc{T) in this temperature 
region, the curvature was strongly underestimated (up to a factor of about 
3) by the fitting function for the as grown and argon annealed samples, even 
if the experimental curvature is partly depressed by the AL fluctuation con- 
tribution. Moreover, this underestimation disappears (within experimental 
resolution) for the sample having the highest oxygen content, which shows 
only a very weak peak (Fig. 20). 

This seems to be a clear indication that the experimental data show 
a divergent trend at Tc which cannot be reproduced by any normal state 
theory. Such a divergence is predicted by the DOS contribution and is coun- 
terbalanced by the AL fluctuation contribution very near (a few K) to T^. 
The gradual reduction in the curvature as the oxygen content is increased 
is consistent with the higher AL fluctuation contribution expected in less 
anisotropic samples. For the oxygen annealed sample the divergence is so 
weak that, within the experimental error, it is indistinguishable from a steep, 
non divergent behavior and the theory of Abrikosov works well even close to 
T 

c- 
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It has been pointed out that when Tq, Ti >> T as is the case for the 
fits on both BSCCO and YBCO, (0) reduces to an activated behavior: 



p,{T)/pab{T) =A[1 + exp (A/T)] (95) 

Here A = Tq—Ti is the energy of the lowest resonant impurity level relative to 
the Fermi level. Only the difference between Tq and Ti is therefore important 
and not their separate values, to which the fit is insensitive when To,Ti >> 
T. Eq. ( ^5]) is very similar to the semi-phenomenological formula for pc 
proposed in P7| : 

Pc{T) = A + BT + {C/T) exp(A/T) (96) 

(where A is some kind of pseudogap). Therefore, in spite of the additional 
1/T dependence and the presence of four phenomenological parameters which 



slightly reduce the discrepancy with experimental data, ( |96|) faces the same 
difficulties as (|9^ in describing the curvature of pc{T) just above Tc in 2212 
BSCCO. 

The above discussion leads to the conclusion that in models based on an 
"activated" behavior for Pc{T) the strong divergence of the exponential at 
T = K simulates the weaker divergence of the experimental data at Tc, 
except in the temperature region very close to Tc. It is then to be expected 
that these models will give a curvature which is too high when applied to the 
low Tc compound of the BSCCO family (2201 BSCCO) for which Tc is only 
about 15 K. Applying (Q) to data [0 on 2201 BSCCO crystals indeed gives 
a very bad fit, as shown in [Q. Moreover, while Eqn. (|9^) predicts a lower 



curvature at Tc than that measured in 2212 BSSCO, the opposite is true for 
2201 BSSCO. 

In conclusion, when Tc ^ 100 K, the curvature of the exponential, which 
diverges at T = i^', is not high enough to account for the apparent di- 
vergence of the experimental pc just above Tc. On the other hand, when Tc 
is closer to K, the exponential divergence becomes much higher than the 
experimental one. 

The analysis carried out in the previous paragraph suggests that it will 
be very difficult for any normal state theory, with a divergence of any kind at 
T = , to reconcile both situations. The same is not of course true of course 
if the peak in pc has a superconducting origin such as the DOS fluctuation 
contribution, whose divergence shifts with Tc, while remaining logarithmic 
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in shape in all cases. Although the semi-phenomenological formula for 
proposed in i.e. ( 0) correctly describes the pdT) curves of the low Tc 



2201 BSCCO phase, this is only true for A = i^, when the exponential 
divergence is canceled and substituted by the 1/T one. In this case, however, 
the formula loses its significance. 

If the DOS fluctuation correction to conductivity is responsible for the 
pc{T) peak just above it is clear that the goal of a normal state theory 
of electrical conductivity along the c axis in layered superconductors is not 
to give a detailed description of this peak. The normal state resistivity 
pNciT) curve must instead lie somewhere below the measured pdT) one. To 
find out what a reasonable behavior of pnc{T) would be, compatible with 
the fluctuation origin of the transverse pdT) peak, in the calculated 



contribution due to fluctuations was subtracted from the experimental Pc{T), 
using ( pi] ) with values for the parameters taken from literature ^ 
I50| , [T5T1 (r = 2 ■ 10-^^ s,T^ = 2- lO'^^ s, and J = AO K), while the critical 



temperature Tc is taken from experimental data. Three simulated pNciT) 
curves were calculated for three different values of the Fermi energy Ep (i.e. 
0.8 eV, 1.0 eV and 1.25 eV). The Fermi energy is just a scale factor for 
the global fluctuation contribution to conductivity cr//, and in order to keep 
the latter within the limits of validity of the theory underlined in section 
6 {afi -C ctat) it must be assumed to be of the order of 1 eV , somewhat 
higher than expected in these materials. These curves are plotted in Fig. 
21 together with the experimental Pc{T) for T < 150 K (for T > 100 K 
the theory is however no longer very accurate since the limit e <^ 1 is not 
fulfilled). It is interesting that because of the less divergent behavior of the 
simulated pNciT) curve as compared to pdT), simpler functional dependences 
for pNdT) could be compatible with it. Some of the theories for normal state 
transverse conductivity which failed to describe the c-axis resistivity peak 
could be in this context reconsidered. 

We conclude that using the fluctuation theory to describe the transverse 
resistivity peak in HTS in zero external magnetic field is well justified. A 
further check of the theory must be sought by adding another parameter, 
besides temperature, on which the resistivity depends. This can be done by 
applying an external magnetic field, as described in the next section. 
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8 The effects of magnetic field 



8.1 Introduction 

The behavior of the resistivity peak under an external c-axis oriented mag- 
netic field is certainly one of the intriguing anomalies of HTS. As the 
field intensity is increased, the position of the peak in pdT) is shifted to- 
wards lower temperatures. However above the peak temperature Tm{B), the 
resistivity curve pc{T) retains the temperature dependence shown in zero 
field above the peak temperature Tm{B = 0). As a result, the magnitude of 
the peak in PciT) increases, and a very strong positive magnetoresistance is 
observed below Tm{B = 0), as shown in Fig. 22. 

The c-axis magnetoresistance shows an even more characteristic behavior 
above Tcq. In contrast to the ab-plane magnetoresistance which is positive at 
all temperatures, along the c-axis the magnetoresistance has been found to 
have negative sign not too close to Tco in many HTS compounds - i.e. BSSCO 
13, |18|, I, H, LSSCO 0, YBCO and TIBCCO ]152[ - and turn to 



positive values at lower temperatures. We will show how these behaviors find 
their explanation within the fiuctuation theory in the magnetic-field induced 
suppression of the AL contribution along the c-axis. 

In this section we will mainly discuss the relevance of fiuctuations in the 
explanation of c-axis magnetoresistance data in HTS. The in-plane magne- 
toresistance is less interesting in the framework of this review, since it has 
already been shown in section 7.2 that in zero magnetic field the DOS con- 
tribution to in-plane transport properties is small, at least not too far from 
T,. Many authors (see e.g. Refs. |5|, |159|, [ISg, |6g, |162|, [151], |15|, [154 



155| , |156|| ) have therefore successfully explained in-plane magnetoresistance 



data in HTS (including BSCCO) using the AL and MT contributions only 



[|64j, 11631, 11651, 11661 . 

After a review of the theoretical predictions in section 8.2, the experimen- 
tal study of magnetoresistance above T^o and the role of the DOS fiuctuation 
contribution in it will be discussed in section 8.3. In section 8.3.1 a short 
survey of in-plane magnetoresistance results is given, then the peculiarities of 
the c-axis magnetoresistance are analyzed in section 8.3.2. Section 8.4 will be 
focused on magnetoresistance effects below Tcq. A summary and discussion 
of the informations given by both zero-field and in-field electrical transport 
measurements presented in sections 7 and 8 will conclude this section. 
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8.2 Theory of c-axis conductivity in magnetic field 
8.2.1 General expressions 

The full theoretical treatment of the effect of magnetic field on the fluctuation 
conductivity of layered superconductors above Tc has been given in |3^. The 



AL, DOS, regular and anomalous Maki-Thompson fluctuation contributions 
to the c-axis conductivity are considered there in detail. We recall here 
the qualitative aspects of the problem and present the results necessary to 
understand the following sections. 

The effect of a magnetic field parallel to the c-axis is considered. For 
this particular field direction the current vertices do not depend upon the 
magnetic field, and both the quasiparticles and the pairs form Landau orbits 
within the layers. The c-axis dispersion remains unchanged from the zero- 
field form. For this simple field direction, it is elementary to generalize the 
zero-field results reported in section 6.3 to finite field strengths. One replaces 

riq^ riC^/i - 2eA)2 (97) 

in each of the integral expressions for the contributions to the fluctuation 
conductivity. The two-dimensional integration over q is replaced by a sum- 
mation over the Landau levels (indexed by n), taking account of the Landau 
degeneracy factor in the usual way ||3^ , 



(27r)2 $0 ^ Attt] 

where n = 0, 1, 2, . . .. So the general expressions for all fluctuation corrections 
to c-axis conductivity in magnetic field can be simply written in the form: 



DOS 

^ zz 



V = f99) 



MT(reg) 
^ zz 



(jMT(an) 



IQr] ^0 [{sb + M {^B + Pn + r)] 1/2 
e'sklS^f eB + l3n + r/2 \ 



eh(3 -fB + pn + r/2 



IQvie - 7^) to V [ilB + HilB + + r)] 
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eB + l3n + r/2 



[{eB + l3n){eB + l3n + r)Yn 



(102) 



where /3 = ArfeB = B/[2T^ \dH^2/ dT\Tc], eB = e + /?/2, and 7^ = 7,^ + p/2. 

For the in-plane component of the fluctuation conductivity tensor the 
only problem appears in the AL diagram, where the matrix elements of the 
harmonic oscillator type, originating from the B\\ (gy) blocks, have to be 
calculated. The other contributions are essentially analogous to their c-axis 
counterparts: 



°° / 1 

4^ h^"" ^ V [{es + Pn){es + /5n + r)] " (^O^) 
2 

+ 



{[es + /3(n + 1/2)] [es + /5(n + 1/2) + r]}i/2 
1 

{[eB + /3(n+l)][eB + /3(n + l) + r]}V2 



cr. 



e2/5(/t + ^) ^ 



4s ^o[(eB + /3T2)(eB + /3n + r)]V2^ 

and 



cr. 



MT{an) 



XX 



8s (e - 7^) ^oM(7B + /3n)(7B + /3n + r)] V2 

1 



[(eB + /3n)(eB + /3n + r)]Vs 



Note that we have included the DOS and regular MT terms together, 
as they are proportional to each other for this conductivity tensor element, 
and the DOS diagram differs from that for the transverse conductivity by the 
factor (vf/sJ)"^, which measures the square of the ratio of the effective Fermi 
velocities. The same factor enters into the normal state conductivity, leading 
to the standard cr^ = lN{0)e'^VpT = Epre^ / {2t:s). We note that (|103|) was 



given previously in [|111|| , using standard procedures |6^, and rederived in 
the formula (|106|) was also given previously ^ . 
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These results can in principle be already used for numerical evaluations 
and fitting of the experimental data. Resistivity curves calculated by (p9|) 
( |102D using reasonable parameter values are shown in Fig. 23. It can be seen 
that the simulated behavior is similar to the experimental one reported in 
Fig. 22. 

Nevertheless, it is useful to manipulate formula ( |102| ) and (|106| ) alge- 
braically, in order to remove the apparent (but spurious) singularity at e = 
7<^. The low field expansions of all results (|99D - (|106|) can be calculated in a 
straightforward way. The analysis of the "strong" field regime e << (3 << 1 
turns out much more sophisticated. In this case, the sums over n for the 
AL and MT terms converge rapidly, and it is enough to keep only the n = 
term in the sums. For the DOS contribution however the formal logarithmic 
divergence of the sum requires a little bit more careful treatment. We will 
deal with these expansions in the next sections. 



8.2.2 Weak magnetic field 

In the weak field regime (/? << e), we expand the various conductivity 
contributions in powers of (3 0]. Such expansions are simplified by using 
the Euler-Maclaurin approximation formula. 



^ /-^ 1 1 , , 

Y: fin) = / f{x)dx + -[/(iV) + /(O)] + -[/ (N) - f (0)] + 
n=0 ^0 2 12 



(106) 



If one writes the expressions in terms of sb, terms linear in B will appear. 
However, writing the expressions in term of the zero-field e, all terms linear 
in B vanish identically, leaving leading terms of order B^. To order B'^{(3'^), 
we find 



AL 
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DOS 



a 



a 



MT(reg) 



e^STK 



16?7 



e STK 



e + r/2 
£(£ + r)]i/2 

ln( 



-1 - 



/jV(£ + r/2) - 
32[£(e + r)]5/2_ 

2 /32(£ + r/2) 



ri/2 + (e + r)iA 



16?7 



(e + r 



,1/2 _ ^1/2 



MT(an) 



16?7 



24[£(£ + r)]3/2 

48[£(£ + r)]3/2 
1 



(107) 
(108) 
(109) 
(110) 
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96He + r)7^(7^+r)]3/2([e(e + r)] V2 + [7 j7,+r)] V^) 



In ( p,lCI| ), we have assumed e « \ and /3 << 1. Typically, 1 T corresponds 
to /5 ~ 10~^ so that the low field expansion is often realized in practice. 
Using (|110|) , one can find the position of the resistive maximum: 



3r 



8k 



167^K 



(111) 



Note that the weak magnetic field reduces by an amount proportional to 

We now present the low-field expansions for the contributions to the fiuc- 
tuation conductivity parallel to the layers. Using the same Euler-Maclaurin 
approximation formula, we obtain 



AL 



e 

16l 



1 



f3^[8e{e + r) + 3r^ 
[e(e + r)]i/2 32[e(e + r)]5/2 



(112) 
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MT(reg) 

XX 



and 



a. 



Mr(an) 



e-7) 

e + r/2 



2 In 



.^ (" + '^ [2 In 

f3\e + r/2) ' 
■24[e(e + r)]3/2. 



.i/2 + (e + r)i/5 



;ii3) 



e^/^ + (e + r)^/^ \ _ / 7 + r/2 
7V2 + (^ + r)i/2j 24V[7(7 + r)]3/2 



(114) 



[e(e + r)]3/2 

The zero-field term in a^^ was first given by Lawrence and Doniach ||115 



and the term of order f3^ was first obtained explicitly in |64] by inverting the 
order of the summation over n and the integration over . 
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8.2.3 The AL and anomalous MT contributions in intermediate 
and strong fields 

We start the discussion of the non-weak magnetic field (/? >> e) from the 
simplest case of the well converging AL contribution. Performing the sum- 



mation in (|102|) - (|106|) one finds 

128r/ 13' 



ai^iP » max{5, r}) = (115) 



in the case of strong field and 



AL/ n , e^s Fr ^ 1 4.57eS Fr 

'-'^ « ^ « - l^lhh l(n + 1/2)13/^ = 

for the intermediate regime which can be realized in the 3Z) case. 
An analogous treatment for the in-plane components results in: 

a4^(/3»max{£,r}) = ^ (117) 

and 

a4^(e«/3«r)^f (118) 

The anomalous MT contribution can be analyzed in the same way, but the 
situation is a little bit more cumbersome. Here is necessary to distinguish the 
cases of strong (J5 » max{e,r, 7^}) and several intermediate field regimes 
(e, 7<^ -C /9 -C r; e <^ /5 <^ 7^,r;) which can be realized in the 3D situation. 
The first limit can be studied in the same way as was done above: the 
expansion of ( |102| ) over /3~ ^permits to evaluate the sum and results in the 
high field asymptotic: 

-f/'-\e,^.,r« /3) = ^ ■ ii:±^, (119) 
In the intermediate cases one finds: 

a^/'-He, 7,, « P«r) = ^ (120) 
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a^^ie « 7^ « /3 « r) = ^ (121) 

and 

^(^"^^ (6 « /3 « 7^ « r) ~ ^ (122) 
The evaluation of the in-plane component gives: 

^.T^^'^He,7.,r«/?) = ^^, (123) 

c^.T^^'^He, 7., « /5 « r) = (124) 

los VP^ 



and 



8.2.4 Renormalization of the DOS contribution divergency in in- 
termediate and strong fields 

As mentioned in H, the fit of experimental data with the theory based on 



the fluctuation renormalization of the one-electron density of states [29, 34 
is excellent for weak magnetic fields but meets noticeable difficulties in the 
region of strong fields. This is due to the formal divergence of the DOS 
contribution to conductivity and to the dependence of the cut-off parameter 
on the magnetic field itself . In this section we clarify the problem of the 
regularization of the DOS contribution in an arbitrary magnetic field 

To avoid the problem of the ultraviolet divergence of the DOS contribu- 
tion with the badly defined cut-off depending on the magnetic field we 
calculate the cut-off independent difference p2| : 



Aa,T^(/5, e) = a^^^P. e) - (0, e). (126) 

For this purpose the zero-field value cr£'^'^(0, e) |Q may be rewritten in 
the form: 
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iu,ej iim 



e srK ^ /-V/J+i/s dn 



/3-0 IQr] 7-1/2 + + l/2)^e + r + p{n + 1/2) 



16r/ „=o-^-i/2 .^e + /5(n + x + l/2)y^e + r + + x + 1/2) 

e^srfi:!/^, Ve + /5n + /3 + Ve + + + /9 
= lim > m , = , 127 

p^o IQn (jn + ^e + r + (3n 

Substituting expression ( |127| ) in ( |126| ), for the difference (T^^^{[3,e) — 
(j£*^'^(0, e) we may write the following formula, where the summation may 
be extended up to — > oo since the sum is convergent (the n-th term of 
the sum is proportional to n^^^^ for large n): 

^(^zz (P> e) = a^, {fj, e) - a^, (0, e) = 

IQV hii P ^ y/e + Pn + y/e + r + pn 

^ =} (128) 



This expression is very suitable for numerical calculation to analyze ex- 
perimental data. It permits to obtain easily the asymptotic behavior of 
magnetoconductivity in the case of non-weak fields (note the inaccuracy 



of the analysis of this asymptotic in [34]). The case of very strong fields 



h ^ max{e, r}, in contrast to |3^], becomes now trivial: it is mainly deter- 
mined just by the logarithmically large contribution of the n = term in 
([128D ( the contribut ion of n > 1 terms is the second in the parentheses below 



Aa£0^(/i » max{e,r}) = ^^{ln— — ^^^^^^=-. + 0,02} 



if] e(A/e + A/e + r) 



-,2, 



'^''"■In^-^^, (129) 



8r] ^ + y/eTr) 
where in the logarithm e = 2, 71 
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Further analysis of ( p.28| ) shows that for the intermediate fields in the 
temperature range of three dimensional fluctuations (e <^ /3 <^ r) the DOS 

contribution shows the magnetic field dependence contrast to ./-^ in 





the above discussed AL and MT contributions (nevertheless the magnitude 
of the former remains smaller that the latters up to /5 ~ r): 

.DOS. - l + yjv/^ 



Aago'ie « /? « r) = ^ Y^i^n ^ ^ ^ - (130) 



:} = 0.428 



2r y/2n + 1 ' ' 16r] 

In addition to the DOS contribution it is necessary to take into account 
the regular Maki-Thompson one, which in the case of in-plane component, 
as we know, has the same sign and functional dependence as the overall 
AcT^'^'^and differs only by the impurity concentration dependent factor k 
(instead of k). 

In weak fields Aaf^'^^^'^^^P <^ r, e) becomes comparable with ( |127] ) only 
in 3D case (e <C r) and for the dirty or intermediate case (Tr ^ 1) when k 
is of the order of k. In the dirty limit (Tr <Cl) k = 2k and Aaf/^'^^s^P <C 
r, e) reaches a half of Aa^f^{P<t^ r, e). 

The evaluation of the regular Maki-Thompson contribution to magneto- 
conductivity may be done by the same procedure as in ( |127| ) for the analysis 
of non-weak fields : 



„=o \/€ + 0{n + l/2)^e + r + /?(« + 1/2) 

-^[^e + f3{n + l)^e + r + f3{n + 1) - (131) 

^Je + (3n ^e + r + (3n]} 
For the 3D case in the region of intermediate fields it leads to 

Aaf/^'-''\e <^ h <^ r) = ^-f^Jjlf^{y'2^rT2-V2^- (132) 



2r .=0 



, ^ } = 0.428 
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One can see that in this region too the contribution Aa^'^^^^^\e <^ P <^ r) 
has the same sign and differs from Aa^J^^{e -C /3 -C r) (see (|130|) ) only by 
the substitution of k by 2k . 

This means that for the 3D situation in the case of strong or intermediate 
(Tr ^ 1) impurity concentration Aa^'^^^'^^\P) has to be taken into account 
side by side with Aaj^^^dS) for all fields up to /5 ~ max{e, r} = r. In the 
case of strong scattering (Tr -C 1) both contributions coincide, while as the 
impurities concentration increases the role of the regular part of the MT 
contribution diminishes. 

The analysis of ( |131| ) in the case of strong field ( (3 3> max{e, r} ) leads to 



A af/^'^^s^ (/5>max{e,r}) = (133) 



16r] ^ ^ 12877 'J 



(134) 



The comparison of ( |134| ) with ( |129| ) demonstrates, that in the 3D dirty case it 



can at least contribute as a constant of the order of 1 in comparison with the 
large field dependent logarithm of (|129|) . In the 2D case its contribution is 
negligible at all both in the absolute value and the magnetic field dependence. 



8.2.5 Discussion 



The results obtained for the c-axis magnetoconductivity are collected in Table 
1. 

Let us start the analysis from the 2D case (r ^ e) (to visualize them 
is enough to skip the third column in Table 1). One can see that the posi- 
tive DOS contribution in the magnetoconductivity turns out to be dominant. 
It growth as i?^ up to -Bc2(e) and then the crossover to a slow logarithmic 
asymptote takes placed. At i? ~ -Bc2(0) the value of Aa^f^{P ~ l,e) = 
— (t£*^'^(0, e) which means the total suppression of the fluctuation correction 
in such a strong field. The regular part of the Maki-Thompson contribu- 
tion does not manifest itself in this case while the AL term can compete 
with the DOS one in the immediate vicinity of T^, where the additional small 
anisotropy factor r can be compensated by the additional in denominator. 
The MT contribution can contribute in the case of small pairbreaking only, 
which is an opposite to the expecting one in HTS. 
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In the 3D case (e <^ r) the behavior of the magnetoconductivity is more 
complex. In weak and intermediate fields the main, negative, contribu- 
tion to magnetoconductivity occurs from the AL and the MT terms. At 
B ~ i?c2(e)(/9 ~ e) the paraconductivity is already considerably suppressed 
by the magnetic field and the dependence of the magnetoconductivity 
changes through the tendency to the high field asymptote —a^/J\0,e). 

In this intermediate region of fields (e <^ /5 <^ r), side by side with the de- 
crease (~ ^J^) of the main AL and MT contributions, the growth of the still 
relatively small DOS term takes place. At the upper border of this region 
{l3 ~ r) its positive contribution is of the same order as the AL one and 
at high fields (r <^ /3 <^ 1) the positive DOS contribution determines the 
slow logarithmic decay of the fiuctuation correction to conductivity which is 
completely suppressed only at -B ~ i?c2(0). The regular part of the Maki- 
Thompson contribution is not of special importance in 3D case. It remains 
comparable with the DOS contribution in the dirty case at fields /3 ^ r , but 
decreases rapidly (~ ^) at strong fields ( /? ^ r), in the only region where 

the robust Aa^'^^{/3, e) ~ In ^ shows up surviving up to /? ~ 1. 

The analogous formulae for the in-pane magnetoconductivity are pre- 
sented in Table 2. 

Analyzing this table one can see that almost in all regions the negative 
AL and MT contributions determine the behaviour of in-plane magnetocon- 
ductivity. Nevertheless, similarly to the c-axis case, the high field behavior 
is again determined by the positive logarithmic A{a^^^ + cr^f^*-^^^-*) contri- 
bution, which the only one to survive in strong field. 

It is important to stress that the suppression of the DOS contribution 
by magnetic field takes place very slowly. Such robustness with respect to 
the magnetic field is of the same physical origin as the slow logarithmic 
dependence of the DOS-type corrections on temperature. This difference 
between the DOS and the Aslamazov-Larkin and Maki-Thompson contribu- 
tions [p.64|| makes the former noticeable in a wide range of temperatures (up 
to ~ 2 — 3Tc) and magnetic fields (~ -Bc2(0)). The temperature scale of the 
suppression of DOS contribution is given by the value of the experimentally 
observed pseudogap. It is of the order of Apseudo ~ 2 — ST^ for magneto- 
conductivity and NMR, Apseudo ~ ttTc for tunneling and Apseudo ^ for 
optical conductivity. 
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8.3 Magnetoresistance above Tc 
8.3.1 In-plane magnetoresistance 

Soon after the discovery of HTS many investigations of the in-plane para- 
conductivity have been performed. It turned out that the early results 
obtained on bulk ceramic samples could be well reproduced in highly ori- 
ented thin films and single crystals ||133|| . A major drawback in the analysis 
of superconducting fluctuations from the paraconductivity is the need for 
separating fluctuation and normal-state conductivity contributions in the 
analysis. Usually, a linear temperature dependence of the resistivity in the 
normal state is postulated, but such assumption is lacking an undisputed 
theoretical background. To overcome the above-mentioned problems with 
the unknown transport properties of the normal-state in HTS, an analysis 
of the in-plane fluctuation magnetoconductivity above Tc has been proposed 
0, |168|| . Early studies of the in-plane fluctuation magnetoconductivity on 



YBCO ceramics, thin films, and single crystals |6|, |T69|, |17|, |T70|, |7T|, [156 
irg, |15|, |T59|, [n^, |I7|, [m], |ri|, in oxygen-deflclent YBCO with = 55 K 

[iTfi , 2212-Bscco ini Q izi B B m], 2223-Bscco pg, ^ 

Tl2Ba2CaCu20x 



were interpreted within the 
1631 , |179|| and clean-limit 



and La2_xSrxCu04 
then available phenomenological dirty-limit |6^, |63t 
theories | 165 | and reanalyzed including non-local corrections [WE . 

Within the phenomenological theoretical context, four different contri- 
butions to Aaxx{B) were predicted. As discussed previously, the orbital AL 
and anomalous MT contributions result from the suppression of the paracon- 
ductivity by orbital interaction with a magnetic field. Due to the anisotropic 
nature of the cuprate superconductors, those effects are prominent with -B || c 
(transverse orientation), but are severely suppressed with B \\ j c (lon- 
gitudinal orientation). In addition, an isotropic Zeeman effect on the spins 
of the fluctuating pairs has been proposed [|163|| , which could dominate the 
magnetoconductivity of the anisotropic HTS with the magnetic field ori- 
ented parallel to the planes and result in two additional contributions, the 
AL-Zeeman and the anomalous MT-Zeeman terms. 

The problem of the amplitude of the anomalous MT contribution in HTS, 
which could not be satisfactory resolved with paraconductivity analysis, was 
re-addressed by several authors with the analysis of the in-plane magnetocon- 
ductivity. Most authors observed that Aa^xiB) with i? || c is dominated by 
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the orbital AL process for e < 0.1, but in optimally-doped YBCO and 2212- 
BSCCO an additional contribution appears at higher temperatures which 
was at first associated with a considerable weight of the anomalous MT con- 
tribution. An analysis within both the dirty and clean-limit phenomenolog- 
ical models revealed several inconsistencies regarding the resulting value for 
r<^ [p.65| , |166| , p.56| , |121|| . This additional magnetoconductivity is suppressed 



in Zn-doped YBCO [|, [ml, partially Pr-substituted YBCO m and in 



the more impure compounds 2223-BSCCO ||160| , |151|| and oxygen-depleted 
YBCO [|7|]. The results obtained by various authors from the analysis of 
the in-plane magnetoconductivity of YBCO vary for ^afe(O) from 11 to 18 
A, for ^c(O) from 1 to 4.6 A and, in those cases, where the anomalous MT 
contribution was detected, for from 2.2 x 10~^^ s to 5.7 x 10~^^ s. For 
YBCO in the 60-K phase the respective values are: ^afe(O) = 25 A, ^c(O) = 
0.9 A. In 2212-BSCCO several authors found ^(0) = 10 to 38 A, ^^(0) = 
0.1 to 2.3 A, and = 0.13 to 7.5 x 10~^^ s. Only few results are available 
for other compounds, like 2223-BSCCO, with ^^^(O) = 16 A, ^^(0) = 1.4 A, 
TlsBaaCaCusO^ with ^^^(O) = 11.8 A, ^^(0) = 0.12 A, = 2.9 x 10"^^ s, 
and Lai.85Sro.i5Cu04 with ^afe(O) = 28 A and ^c(O) = 1.4 A. In general the 
values obtained from the in-plane magnetoconductivity analysis agree well 
with those from other experimental methods, but the scatter in the results 
for T(j, is quite large. 

On the other hand, measurements of the in-plane magnetoresistance in 
YBCO single crystals near and above room temperature revealed an un- 
conventional Ap/po = {pxx{B) — Pxx(0))/pxx(0) ~ T~'^ temperature depen- 
dence, which the authors attributed to the cyclotronic motion of normal-state 
quasiparticles. They proposed that in a fourfold-symmetric, two-dimensional 
metal the orbital magnetoresistance and the Hall effect measure the variance 
and the mean of the local Hall angle along the Fermi surface, respectively. 
Hence, the temperature dependence of the normal-state magnetoresistance 
is Ap/po = A tan^^H? where tan^H is the Hall angle in the material and 
A = 1.7 for YBCO ||180|| . However, superconducting fluctuations were en- 
tirely neglected in this analysis. Using these results and magnetoresistance 
data from YBCO thin films it was shown that neither the fluctuation pic- 
ture alone nor the normal-state magnetoresistance can account for the data 
from about Tc to above room temperature ||187|| and it was proposed that 
the quasiparticle magnetoresistance resembles the anomalous MT contribu- 
tion ||156| , |154|| . In Fig. 24 it can be noticed that in fact the combination of 



87 



the four fluctuation contributions, the orbital AL and anomalous MT, and 
the respective Zeeman effects (ALO+ALZ+MTO+MTZ) fails to fit the data 
above 2Tc. It should be noted that using the full theory outlined in the pre- 
vious chapters (including the DOS and regular MT expressions) instead of 
the phenomenological model would not resolve this discrepancy. As an al- 
ternate scenario it was proposed that the magnetoconductivity close to is 
dominated by the orbital AL term, with a crossover to the normal-state mag- 
netoconductivity Aaxx — —Ap/p^ at e > 0.2 (see Fig. 24). In this case the 
AL Zeeman and the DOS contributions, which are expected to gain impor- 
tance relative to the orbital AL term at higher temperatures are masked by 
the normal-state magnetoresistance. Similar results were recently obtained 
on YBCO single crystals ||121[| . In this context the large scatter of the results 
for and the apparent absence of the anomalous MT contribution in several 
materials rather reflect different shapes of the Fermi surface, resulting in a 
variation of the parameter A, i.e. a different ratio between the normal-state 
magnetoresistance and the Hall angle. 

As can be seen from the previous paragraph, the DOS contribution is 
unlikely to have significance for the in-plane magnetoresistance of HTS as 
far as 5 || c is concerned. In the longitudinal orientation however, the 
orbital terms are suppressed due to the anisotropy of the cuprates and the 
normal-state magnetoresistance due to the absence of the Lorentz force on the 
quasiparticles, respectively, and, thus, the Zeeman terms become dominant. 
The results available for this geometry [[16|, [l7l, HM HHl, Ull, HH, |153 



175|| are not entirely conclusive, probably due to the required very accurate 



orientation of the magnetic field parallel to the Cu02 planes. Some authors 
report a good accordance with the AL Zeeman term only. The anomalous MT 
Zeeman effect can be supposed to be negligible if the corresponding orbital 
anomalous MT contribution is not strong, a fact which is now well agreed 
for the cuprates. As mentioned in section 8.3.2, a DOS Zeeman process 
is needed to explain the longitudinal out-of-plane magnetoresistance data 
above T^. Accordingly, the DOS Zeeman term can be expected to induce a 
sign change from negative to positive magnetoconductivity at temperatures 
considerably above T^. Apart from a single report ||153|| this has not been 
observed so far. Certainly the application of the DOS effect to the in-plane 
magnetoresistance needs more attention in future. 

Finally, two other effects commonly associated with the in-plane magne- 
toresistance are worth mentioning. It was shown that the sample shape, in 
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particular the usually almost square geometry of single crystals, has consid- 
erable influence on the in-plane magnetoresistance at higher temperatures 



due to bending of the carrier trajectories | |181| . Thus, evaluations of the 
magnetoresistance in this temperature region have to be regarded with some 
caution. On the other hand, a large increase of the in-plane magnetoconduc- 
tivity above theoretical predictions close to T^, which was observed by several 



authors, can be naturally explained by a non- uniform T^. in the samples [182 



It was shown that measurements of the magnetoresistance in fact can serve 



as a very sensitive probe for Tc inhomogeneities in HTS [p. 83 



8.3.2 Out-of-plane magnetoresistance 

We will now address the problem of the origin of the negative out of plane 
magnetoresistance above T^o, discussed in the Introduction. The negative 
c— axis magnetoresistance of BSCCO single crystals was first observed in 
[185] and then carefully measured in [3^ (see Fig. 25) at temperatures above 



95 where it was qualitatively interpreted in terms of the holon and spinon 
model by Anderson ||124|| . However, this analysis is not quantitative and is 



based on the phenomenological assumption of an "activated" behavior for 
pc{T), which in the previous section has been shown to be unsatisfactory. 
On the other hand, from the data reported in it can be seen that 



the effect becomes significant below approximately 140 K and its magnitude 
increases dramatically as the temperature goes down to 95 K. This temper- 
ature range is the same as that observed for the fluctuation induced positive 



ab-plane magnetoresistance, which naturally leads [^ to the hypothesis that 
fluctuations are also responsible for the negative c-axis magnetoresistance. 
Since the DOS fluctuation contribution is held responsible for the peak in 
Pc(T), its contribution to magnetoresistance cannot be neglected and may de- 
termine its sign. This contribution is indeed expected to be negative in sign, 
since a suppression of the DOS contribution by the magnetic field would give 
a decrease of the resistivity. In the temperature region where the DOS con- 
tribution dominates over the AL one, a negative fluctuation induced c-axis 
magnetoresistance is therefore conceivable. All the features of the observed 
magnetoresistance are therefore consistent with its attribution to fluctua- 
tions, with a key role played by the DOS contribution. Starting from ( |110| ) 
the following expression for the fluctuation c-axis magnetoresistivity close to 
Tc in the presence of weak magnetic fields (this assumption is fulfilled in the 
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experiment reported in |]3^) has been found ||40|] : 

2.92 ■ 10^^ Pc{B, T)f{T)B\ (135) 



p,(5,T)-p,(0,T) 



Pc(0,T) 

The temperature dependence is mainly included in the factor /(T): 



f(T) = f 3(^ + r/2) 

•^^^ [£(£ + r)]3/2 \ e(£ + r) 




+ 



(136) 



2 

+ - 



(e + 7^ + r) {^(e + r) + 7^(7^ + r) + [e{e + r)7^(7^ + r)]^/^} 
Ml^ + r)f' [{e{e + r)y/' + (7^(7^ + r))V2] 

Here cgs units are used except for the magnetic field B (measured in 
Tesla) and the resistivity pc{B,T) (measured in Qcm). 

The first term in (|136|) represents the AL contribution, the second is the 
sum of DOS and regular MT contributions and the third is the anomalous 



MT one. The fit of ( |135| ) with the experimental data was performed in 
using as adjustable parameters vp, t and the phase pair-breaking lifetime 
r^. The values of the interlayer spacing s ~ 10~^cm and of the hopping 



integral J ^ AO K were taken from literature data ||3^, since they are not 
likely to vary strongly from sample to sample (at least for BSCCO samples 
with metallic behavior far from Tc), while Pc{B, T) and Tc = 85 K have been 
extracted from the experimental curves. 



The results of the fit performed using ( |135D for the magnetoresistance 



curves are shown in (Fig. 26). The curves measured at T = 95 i^' and 
T = 100 K were fitted simultaneously (i.e. using the same values of the 
fitting parameters for both curves) to put more constraints on the fitting 
procedure. Those measured at T = 105 and higher temperatures were 
not considered in the fit because they lie outside the temperature region 
e <^ 1 in which the theory is quantitatively accurate (at 105 K,e = 0.21). 
However, the theoretical curve at 105 K has been drawn in Fig. 26 using 
the values of vp, t and r<^ found for the curves at 95 K and 100 K in 
order to show that, even at higher temperatures, the calculated temperature 
dependence of the transverse magnetoresistance is in substantial agreement 
with the experimental data. 
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The values of the fitting parameters were vp = 3.1 x 10^ cm/s, r = 
1.0 ■ 10-1^ s and = 8.7 ■ 10"^^ s. 

While the field dependence of the magnetoconductivity is simply B"^, 
at least for moderate fields, its behavior with temperature, given by ( |136| ), 
is much more interesting. On the basis of (|136|) and using the above val- 
ues for the fitting parameters, the expected temperature dependence of the 



transverse magnetoresistance has been calculated in |^ (see Fig. 27). The 



existence of an inversion temperature at which a reversal of the sign of 
magnetoconductivity occurs is predicted. The physical origin of this change 
of sign is the same as for the appearance of the peak discussed in Section 6.1 
: relatively far from Tc the AL negative magnetoconductivity is suppressed 
by its dependence on the square of the transparency and the positive DOS 
contribution dominate, while very close to the very singular temperature 
dependence of the negative AL contribution (~ e^^) makes it prevail over the 
less singular DOS one (~ e^^), despite the latter's linear dependence upon 
transparency. 

Last years the problem of the out of plane magnetoresistivity has been 
studied experimentally very carefully [|, |186| , |187| , Hashimoto et al. 
measured the c-axis magnetoresistance of BSCCO single crystals in fields up 



to 30 T and found results similar to ^T^, except that at very high fields/low 



temperatures the field behavior is no longer concave, as can be expected from 
[0]). They fitted the experimental data to the original theory of loffe et al. 
[p9|l (which is valid only in the clean case) but including magnetic field and 
renormalization effects on the mass term. The agreement between theory 
and experiment is good for fields up to about 10 T , although at higher fields 
the measured magnetoresistance is higher than expected. Fitting parameters 
were the ratio Tc/Tco = 0.84 (Tco being the mean field critical temperature), 
the in-plane coherence length ^||(0) = 16 A and the effective mass rriab = 3.6, 
T having been fixed to 1.2 ■ 10~^^ s to match the clean limit assumption. 

In | |186| | the negative c-axis magnetoresistance of BSCCO single crystals 
was again observed at temperatures higher than 100 K. The fit of experi- 
mental data was performed using AL and DOS terms only, with parameters 
Tc = 93.5 K , T = 1.5- 10"^^ s, ^||(0) = 14 i and ^±{0) = 0.12 A. The 
unphysically low value of ^_l(0) was attributed to the possible overestimation 
of the intrinsic resistivity due to the influence of microcracks perpendicular 
to the crystal c-axis. The fit was anyway very good, except for a slight over- 
estimation of the effect at the lowest temperature used (T = 100 K). In a 
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later paper [|187]| the authors extended their analysis to include both longitu- 
dinal and transverse in-plane magnetoresistance data measured on the same 
sample in the fit. With the addition of an AL Zeeman term (necessary to de- 
scribe in-plane longitudinal magnetoresistance, in which the other terms are 
suppressed), the experimental data were described very well for 5 = 12 T in 
an extremely wide temperature range (up to about 2Tc) using as parameters 
vf = 2.2 X 10^ cm/s, r = 1 ■ lO'^^ s and J = AK. 

In the transverse magnetoresistance was measured down to tempera- 
tures close to Tc in two YBCO single crystals in different oxygenation states 
and with different twin densities. After corrections for the inhomogeneous 
current distribution the measured magnetoconductivity for B = 12 T was 
fitted with ( |102| ), employing a weighted cutoff for the sum in the DOS 
contribution to smooth the field dependence of the calculated values. Its 
temperature dependence was found to follow the behavior shown in Fig. 
28 and the existence of a sign reversal temperature was unambiguously 
confirmed at about 10 K above Tc. Assuming vp = "2 ■ 10^ cm/s and 
the temperature dependence r = ~ 1/T the fitting parameters were 
r(100 K) = T^{100 K) = {A±1)- 10^^^ s, J = (215 ± 10) K. 

We point out that although YBCO generally shows a weak transverse 
resistivity peak because of its relatively small anisotropy, the different tem- 
perature dependences of the AL and DOS terms allow for a negative total 
magnetoresistance. This occurs at higher temperatures than for BSCCO, 
making the observation of the change of sign easier. 

In the sign reversal of magnetoresistance in HTS was proved, but 
the smallness of the negative magnetoresistance did not allow for a precise 
quantitative check of the predictions of the fluctuation theory. Very recently, 
similar measurements have been performed on single crystals of the very 
anisotropic Tl 2223 compound [|152|| . A strong negative magnetoresistance 
is observed a few degrees above initially increasing as temperature is 
decreased and then turning to positive as is approached. In this extremely 
anisotropic compound this effect is very pronounced, so that the relative 
experimental errors on the data are very small. It turns out that the fit of the 
data with the fluctuation theory including the DOS contribution is excellent, 
the theoretical behavior being perfectly reproduced with quite reasonable 
values of the fitting parameters. 

There is therefore now increasing evidence that the negative c-axis mag- 
netoresistance observed in many HTS is really connected with the DOS fiuc- 
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tuation contribution. The competition among the DOS and AL fluctuation 
contributions leads to the negative c-axis magnetoresistance observed above 
Tc and to its change of sign in the vicinity of T^. 



8.4 Fluctuation magnetoresistance below Tc 

In this section we analyze the role of the DOS fluctuation contribution in the 
behavior of pc{H,T) below the zero field critical temperature Tco, and the 
associated increase of the c-axis resistance peak with magnetic field. This 
effect was first measured by Briceno et al. ^ and then analyzed by several 
authors. 

The main features observed are: 

1) a well pronounced maximum at zero field in the pc{T) behavior; 

2) a very strong shift of the zero resistance temperature in external mag- 
netic field (ATc ~ 30 isT already for 5 ~ 1 T) 

3) a large broadening of the peak for nonzero magnetic fields with the 
appearance of a long resistivity tail for higher fields and a large associated 
positive magnetoresistance. 

Kim and Gray |^ explained the broadening of the peak with increasing 
magnetic field in terms of Josephson coupling, describing a layered supercon- 
ductor as a stack of Josephson junctions. The growth of the resistance above 
the transition was attributed to the DOS fluctuation contribution. The re- 
sults obtained by the Kim and Gray model are interesting, and its agreement 
with experimental data is good, as seen in Fig. 29, but their approach intro- 
duces another element in the discussion beyond fluctuations and lacks a full 
microscopic foundation. The field dependence of pc is described assuming 
that phase slips in layered superconductors are identical to phase slips in a 
single Josephson junction having area ^q/ B. It has been remarked ||188|| that 



this model is inadequate to describe the field dependence of the activation 
energy extracted from resistivity curves in the low temperature region. 

In [^] the first attempt to describe the experimental pc{T, H) curves for 
perpendicular magnetic fields up to 8 T using only the fluctuation theory was 
undertaken. The first step was to fit the zero-field data in the temperature 
range from 92 K to 115 K (which corresponds to 0.02 < e < 0.25 ) 



with the fluctuation theory using ( p.02| ). The values obtained for the fitting 
parameters were r = (5.6 ± 0.6) • 10"^'^ s, = (8.6 ± 1.4) • 10~^^ s, Ep = 
(1.07±0.12) eV and J = {A3 ± A) K (T^o = 89.8 K was taken as the midpoint 
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of the transition). Keeping these parameters fixed, the data measured in 
a magnetic field were then fitted using the field parameter /? as a fitting 
parameter. It was found that to reproduce the shape of the peak the critical 
temperature Tc{B) also had to be used as a fitting parameter, because the 
mean-field Tc{B) incorporated with ( |102| ) appeared to be unable to describe 
the strong shift of the peak in magnetic field. With this correction the theory 
satisfactorily describes the resistivity behavior for weak fields B < 1 T, while 
at higher fields the strong broadening of the peak cannot be reproduced. It 
appears from Fig. 30 that (|102|) are able to fit the experimental data only for 
temperatures above the temperature Tm{B) corresponding to the maximum 
in pc- For T < Tm{B) the experimental data decrease much more slowly than 
expected from the theory, and this discrepancy increases with increasing field 
strength. 

This lack of agreement is not surprising. The first reason is the absence 
of critical fiuctuations in the above treatment. As it was shown in Refs. 
189| , |190|| , the extension of fiuctuation theory beyond the Gaussian approxi- 



mation results in a shift of Tc{B) to zero temperature with a corresponding 
broadening of the transition, which increases strongly with increasing mag- 
netic field strength in qualitative agreement with experimental data. Another 
reason for peak broadening is the presence of a complex vortex structure 
which leads to an additional dissipation in the mixed state. However, it can 
be noticed that the broadening of the superconducting transition in HTS 
in the presence of an external magnetic field is quantitatively very similar 
for in-plane and c-axis experiments. This in spite of the quite pronounced 
qualitative difference, i.e. the presence of the resistivity peak in c-axis mea- 
surements as opposed to the monotonic decrease of in-plane resistivity with 
decreasing temperature. This suggests a common, intrinsic origin for this 
broadening, which in view of the results reported above and in the previous 
sections, could be attributed to fiuctuations. 

The possibility of describing the shape of the transverse resistivity peak 
in presence of an external magnetic field within the fiuctuation theory using 
the Hartree approximation instead of the Gaussian one was analyzed in . 
Here, (|102|) were modified within the self-consistent Hartree approximation 



developed by Ullah and Dorsey ||19CI|| , in which is renormalized according 
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to self-consistent equation [p.9CI| 



SB = £b- Gf t 



[(£B + /5n)(£B + /3(n + l)+r)]V2 



1 



(137) 



n=0 



where Gi = 2nTc{0)/ Hc{0)sC,lij{0) is the Ginzburg number (-f^c(O) is the zero- 
temperature thermodynamical critical field and C,ab{^) is the zero-temperature 
Ginzburg-Landau coherence length in the afe-plane). 

In this paper the authors report simultaneous measurements of both the 
in-plane and the out-of-plane resistivities of BSCCO films in magnetic fields 
up to 1 T. Measuring both pab(T,B) and pc(T,B) on the same films allows 
them to use the same set of microscopical parameters for all curves, thus 
putting more constraints on the fitting procedure. They found that both 
Pab{T, B) and Pc{T, B) are in remarkable agreement with the fiuctuation the- 
ory in a very wide range of temperatures both below and above Tc. This 
gives a strong evidence in favor of the fiuctuation origin of the transverse 
resistance peak. 

Figures 31(a) and 31(b) show the PabiT) and pdT) dependences respec- 
tively for several applied magnetic fields in Ref.(||36||). The curves were fitted 
with the formula p = 1/{<Jn + o"//). Here aji = a"^^ + a^^^ is the fiuctua- 
tion conductivity obtained replacing sb with in ( |102| ) and neglecting the 
Maki-Thompson contribution, while Pn = l/c"7v is the normal-state resistiv- 
ity which was assumed to be linear and extrapolated from the temperature 
range 150-200 K for both components of resistivity. The zero-field curves 
PabiT, 0) and pciT, 0) were fitted by the above equation with two free param- 
eters, s and r. This procedure gave s = (15 ± 2) A and r = (5.0 ± 0.6) ■ 10~^ 
for both directions. These parameters were then kept fixed while the fit was 
performed for finite-field data using as free parameters /? (field dependent) 
and Gi (field independent). This procedure provides a uniquely determined 
set of parameters for both components of the resistivity tensor. It was found 
that Gi = 0.12 ± 0.01 for all curves, in agreement with estimates from mi- 
croscopic parameters, while P = B /[2Tc\dHc2/dT\Tc] (Eq-(35) in must 
be corrected) was correctly found to vary linearly with the applied field, as 
expected, yielding a slope of the upper critical field of about IT/K at Tc. 
The agreement between theory and experimental data is remarkable, given 
the number of curves fitted the small number of free parameters used and 
the considerable qualitative and quantitative differences among in-plane and 
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transverse resistivity behaviors. The fit is unexpectedly good even down to 
temperatures at which the resistivity falls below about 10% of normal state 
value, i.e. a region where vortex dynamics is generally assumed to be the 
main source of dissipation. 

Soon after similar conclusions were reached using a different approach 
[|191|| . It was reported that the ab-plane resistivity of YBCO films and crys- 
tals in external magnetic field could be well fitted down to low resistance 



values using the fluctuation theory developed by Ikeda et al. ||192|| in the 
Hartree approximation. In the same region a scaling law of the resistivity 
curves for different orientations of the magnetic field is valid. This, and the 
comparison of the similarities and differences among the behavior of films 
and crystals, leads to the conclusion that in a wide region of the H-T plane 
the dissipation is due only to intrinsic properties (fluctuations) and only at 
lower temperatures pinning-related dissipation processes take place. 

There is convincing evidence, therefore, that the fluctuation theory in 
the Hartree approximation is able to describe very well the experimental 
Pab{T, B) and pdT, B) curves, using consistent values for the fitting parame- 
ters for both current flow directions. There are clues that highly anisotropic 
layered superconductors can be described as a normal, strongly fluctuating 
phase rather than a superconducting mixed phase down to very low temper- 
atures, possibly leading to a reconsideration of the role of vortex dynamics 
in these compounds. 

To conclude this section, we want to mention again the paper by Nakao 



et al. [185] in which the resistance of a BSCCO single crystal was measured 
by sweeping an external magnetic field up to 40 T at a constant temperature. 
Below Tco a peculiar behavior was observed, the resistance initially increasing 
steeply with field, then levelling off and finally decreasing in very high fields 



(Fig. 32(a)). This behavior (later observed also in [^, |193|| ) can be ascribed 
to the combined effect of the magnetic field induced destruction of the phase 
coherence between layers and the suppression of the DOS. Indeed, curves 



simulated using the fluctuation theory in ||185|| show a behavior very similar 
to the observed one (Fig. 32(b)). 

8.5 Final remarks 

At the end of this survey of conductivity and magnetoconductivity measure- 
ments explained in terms of the DOS fluctuation contribution, we want to 
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make some comments about the values of the parameters extracted from the 
fits by the various authors. 

The variety of experiments on c-axis electrical transport properties which 
can be quantitatively explained by fluctuation theory alone convinces one 
of the substantial correctness of this approach. It must also be stressed 
that in this fluctuation theory all parameters have a well defined physical 
meaning, no phenomenological parameter ever being used. The reasonable 
values obtained for the physical parameters involved in the theory are an 
important check of the consistency of the theory itself. However, the values 
of the microscopic parameters found in the literature are not always fully 
compatible. 

There are many obvious reasons for a certain scatter of the values of mi- 
croscopic parameters found by comparing the results of several papers, apart 
from the intrinsic uncertainties. It can be ascribed both to differences among 
the samples used (single crystals vs. films, different oxygenation states) and 
to differences in the methods used to compute pc from voltage drop measure- 



ments (see e.g. g). Some uncertainty in the values of the fitting pa- 
rameters is due to the rather arbitrary choice of additional parameters which 
are simply deduced from experimental data or taken from literature data 
and, apart from magnetoresistance measurements, the linearly extrapolated 
normal state resistance behavior. Also, the weak temperature dependencies 
of r and in the narrow temperature range considered have been neglected 
by some authors but estimated as ~ 1/T by others, and formulae containing 
different approximations have been used in different papers. Finally, other 
physical effects can partially contribute to the observed phenomena. 

In spite of this, we generally find a good overall agreement among the sets 
of parameters. For BSCCO (there are not enough papers concerning YBCO), 
r is the range l-^5■ 10^^^ s, while varies over one order of magnitude from 
8 ■ 10~^^ s to 8 ■ 10^^'^ s. However, only enters in the small Maki-Thompson 
term, so that the uncertainty on this parameter is likely to be rather large. 
Ep and vp have a similar role in the fits (they define only the scale for 
the fluctuation conductivity) but cannot be directly compared. To provide 
indicative reference values for Ep we repeated the fits performed in Refs 
[^, 40] using the same formulae as in [^|. Ep is introduced in the formulae 



instead oi vp through the expression for the normal state conductivity, thus 
avoiding this source of error too. It turns out that for ^2|] (resistivity data 
on films) r = 3.0 ■ 10"^^ s, = 3.6 ■ lO'^^ s and Ep = 1.07 eV, while for the 
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magnetoresistance data of single crystals analyzed in |^ r = 0.9 ■ 10"^*^ s, 
= 7.8 ■ 10"^"^ s and Ep = 0.25 eV, values to be compared with those found 
in A big discrepancy is instead found among the values for J. However, 
it must be noticed that the estimate J = 40 K used by some authors was 
deduced from in-plane measurements in overdoped samples ||32|, p.94|| , while 
the fitting value J = 4 K refers to c-axis measurements on samples optimally 
doped for the highest Tc and showing a very pronounced peak. 

We have collected the parameters following from the analysis of the c-axis 
magnetoresistance in Table 3 0. 

It can be seen that r is almost always close to 2 • 10~^^ s, which is a 
widely accepted value in the intermediate region between clean and dirty 
limits, while the ratio r^/r is always about 10, in good agreement with the 



expected one The values imply that the pair-breaking in BSCCO 
has a moderate strength, which is consistent with several papers (see sec- 
tion 8.2) where it was shown that the AL term alone describes adequately 
the fluctuation contribution to the in-plane conductivity in BSCCO. The 
regular MT contribution is negligible in the clean and intermediate cases, 
while the relative importance of the anomalous MT contribution depends 
upon the relationship between 7^ and r. For the parameters found above 
'^^p '0 r which still corresponds to moderate pair-breaking, but closer to 
the weak limit. This is the reason why it usually suffices to consider only 
the AL contribution. As for the Fermi energy value, which being a scale 
factor is most influenced by the experimental conditions (oxygenation state, 
methods for calculating resistivity, etc.), it lies within a factor of 2 around 
0.5 eV , in good agreement with the estimate given in |Q and not too far 
from that calculated using band structure calculations and photoemission 
measurements. 

Comparing the values of parameters extracted using ( |102| ), it must also 
be noticed that these equations are used to describe both the zero-field cor- 
rections leading to the existence of the resistivity peak (i.e. a change in 
conductivity of the order of 20% in typical BSCCO films) and the small 
magnetoresistance corrections (i.e. an effect of the order of 0.1%), thus span- 

means that the parameters where recalculated by the authors of the review on the 
basis of the complete theory presented here; (x) means that this parameter was not used 
as a fitting one but was taken by authors from literature;x means that this parameter has 
been recalculated by the authors of the review from the original parameters presented by 
the authors of the article 
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ning more than two orders of magnitude. The consistency of the values for 
the fitting parameters, and especially those which only give the scale of the 
corrections (e.g. Ep) is much better than that. 

We can conclude that the number of unexpected experimental facts con- 
cerning the electrical transport properties of HTS for which the fluctuation 
theory provides a quantitative explanation without any other additional as- 
sumption is impressive. Moreover, the theory uses no phenomenological pa- 
rameter and, apart from a reasonable scatter of the values of the physical 
parameters involved, these values are certainly acceptable. No other theory 
at present has achieved so many successes in this field. 
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9 The fluctuation induced pseudogap in HTS 



9.1 Introduction 

The observations of pseudo-gap like phenomena in the normal state of HTS 
is currently one of the most debated problems. Among these an important 
place is occupied by measurements ^ of the c-axis reflectivity spectra, 
in the FIR region on YBa2Cu40s single crystals. With the decrease of 
temperature the c-axis optical conductivity decreases showing around ISOi^' 
a transition from a Drude-like to a pseudogap-like behavior. The value of 
gap is Uq ~ 180cm~^ and it seems weakly dependent on temperature. The 
decrease of temperature does the gap structure more pronounced without any 
abrupt change at the superconducting transition temperature Tc = 80K. 

There are many hypotheses concerning the origin of this pseudo-gap. We 
will show below that it can also be explained by the suppression of the 
one-electron density of states on the Fermi level due to the interelectron 
interaction. It will be demonstrated that the relatively low electromagnetic 
wave frequencies suppress the main positive AL and MT terms, while leaving 
almost unaffected the negative DOS contribution. This means that noticeable 
changes will occur in the electromagnetic wave reflectivity at frequencies well 
before the value u) ~ r~^expected for a Drude term. The positive AL and 
anomalous MT contributions show frequency dependence when uj ~ min{T — 
Tc,r~^}[|^, |195|| , whilst the negative DOS contribution shows dependence 



on scale u ~ min{T, r~^}. This competition results in the rapid decay of 
the dissipative processes at frequencies of the order of u ^ T — Tc and in 
the appearance of a transparency window up to u; ~ min{T, r^^}. The high 
frequency behavior of Re [cr(to')] is mostly governed by a"'{u!) which decreases, 
in agreement with the Drude law, for uj > t~^. 

Below we will study the a.c. fluctuation conductivity tensor for a layered 
superconductor taking into account all contributions and paying attention 
to the most interesting case of c-axis polarization of the field. Nevertheless, 
for completeness, the a6-plane results of the old paper of Aslamazov and 
Varlamov will be re-examined and discussed in application to the novel 
HTS layered systems. 
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9.2 Paraconductivity 



The optical conductivity of a layered superconductor can be expressed by 
the same analytical continuation of the current- current correlator (electro- 



magnetic response operator) Q 



(fi) 

al3 



as in the section 6 but in contrast to 



d.c. conductivity case, calculated without the assumption uj ^ Q: 



Re [cra/3(u;)] 



Im 



:i38) 



One can see that the calculations are analogous to those for the d.c. fluc- 
tuation conductivity, discussed in section 6.4, but here we cannot treat the 
external frequency as a small parameter. We will use the same model and 
definitions introduced in section 6. 

Let us start with the AL contribution (diagram 1 of Fig. 9) to the a.c. 
fluctuation conductivity. The general expression for the appropriate contri- 
bution to the electromagnetic response operator as a function of the Matsub- 
ara frequencies of the external electromagnetic field lOy is defined by (|67|). 

In the vicinity of Tc, for frequencies u -^T, the leading singular contribu- 
tion to the response operator Q^p^^^ arises from the fluctuation propagators 
in (^7|) rather than from the frequency dependencies of the Ba blocks, so it 
suffices to neglect its frequency dependencies |55|, as was done for the d.c. 
case. Carrying out the same calculations as in section 6.4 but without the 
limit C(J — >■ 0, one can find the explicit expression for the imaginary part of the 
retarded electromagnetic response operator |^ for real frequencies u <^T: 



Im 
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(139) 
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2/ V l6Tc 



^^It is necessary to mention that the direct calculation of the expression ( leads to 
the appearance of divergent expressions in the ReQ. Nevertheless, as it was shown in |55| , 
the thorough summation of all diagrams from the first order of the perturbation theory 
for Q before momentum integration leads to the exact cancellation of Re [(5^(0)] . This 
fact justifies the possibility of the further calculation of different diagrammatic terms to 
separately. 
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where (2) = 2 In [^/i + + r)J , Z^al^) = -^z{z + r), ADi{z) = Di {z)- 
Di (e), AD2 {z) = D2 (z) - D2 (£)and r = AriJ^v^p. 

The expressions presented above solve the problem of the frequency de- 
pendence of the paraconductivity tensor in the general form for e <C land 
uj < Tfor an arbitrary relation between e, rand uj, but they are too cumber- 
some. 

We concentrate therefore on the most interesting case of 2D fluctuations 
where ^ciT) -C s (r <C e). af-^ turns out to be suppressed by the necessity 
of the independent tunneling of each electron participating in the fluctuation 
pairing from one CUO2 layer to the neighboring one. The approximation 
r -C e simplifles considerably the expressions ( |139| ) and ( |140| ) while still 
giving valid results for frequencies comparable to T^. 
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The second expression coincides with that one obtained in 195 1 



(142) 
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Two facts following from the expressions obtained should be stressed. 
First, the paraconductivity begins to decrease rapidly with the increase of 
frequency for u > T — T^ (the critical exponents of this power decrease coin- 
cide with those of the e-dependence of d.c. -conductivity tensor components: 
z/|| = 1{2D fluctuations) and z/_l = 2{0D fluctuations)). Secondly the assump- 
tion that one can neglect the cu-dependence of the Green's functions blocks 
evidently breaks down at frequencies u >T and has the effect of accelerating 
the decrease of paraconductivity. 



9.3 Density of States contribution 



The four main diagrams for the DOS contribution to the electromagnetic 
response tensor are diagrams 5-8 of Fig. 9. The general expression for the 
DOS contribution to Qapi^^) from diagram 5 is given by expression ([73|) . 
The external frequency enters in this expression (|73|) only by means of 
the Green's function G(p, ujn+u) and it is not involved in q integration. So, 
near Tc, even in the case of an arbitrary external frequency, we can choose 
the propagator frequency u;^ = 0. The contribution of diagram 7 of Fig. 9 
can be treated in the same manner. 

The diagrams 5 and 6 of Fig. 9 are topologically equivalent and this 
fact would lead one to believe that they give equal contributions to a{uj). 
Nevertheless the thorough analysis of the analytical continuation over the 
external frequency shows that their contributions differ slightly and for the 
total DOS contribution to conductivity one can find: 
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In contrast to ( |139| ) and ( p.4CI| ) , this result has been found with only the 
assumption e <^ 1, so it is valid for any frequency, any impurity concentration 
and any dimensionally of the fluctuation behavior. The function k{uj,T,T) 
can be easily used to fit experimental data. Nevertheless we also present the 
asymptotics of the expression ( p.44| ) for clean and dirty cases. In the dirty 



case 



and in the clean case 
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9.4 Maki- Thompson contribution 

The total contribution of the MT-like diagrams to aa/sioo) has been analyzed 
in for the case of zero frequency and the frequency dependence of (T|j (u) 
has been studied in |5^. In section 6.4 we have demonstrated that the 



regular part of the MT diagram can almost always be omitted. So we will 
not discuss it in this section and will concentrate instead on the anomalous 
MT contribution only. 



In 1^ it was demonstrated that in the case of quasi two-dimensional elec- 
tron motion (|3^ there is no formal necessity to introduce the pair-breaking 
time because the Maki-Thompson logarithmic divergence is automatically 
cut off by interlayer hopping. Nevertheless, all evidences show that the intrin- 
sic pair-breaking in HTS is strong (at least one of its sources may be identified 
as thermal phonons) and an estimate of the appropriate ~ 2 5 ■ 10~^^s 
is only several times larger than T~^.So we are actually in the overdamped 
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region of the MT contibution, and the latter does not noticeably affect the 
e dependence of conductivity (the major part of experimental results is 
explained in terms of AL or AL and DOS contributions). However we are 
still interested in the MT contribution because of its frequency dependence 
which evidently determines another characteristic scale in addition to the 
previous three (T — T^, T, r~^) we have introduced: ujmt ~ ■ 

We start, as usual, from the general expression for the anomalous MT con- 
tribution to the electromagnetic operator tensor (|80|) and after the integra- 
tion over momentum p and the summation over ujn in the range uJn G [—^v, 0[ 
{anomalous part), one finds: 
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+ Dg^j is + riq^ + r sin {q±s/2) 



At this stage of calculation we artificially introduce the phase-breaking 
time in the "Cooperon" 

vertices. Carrying out the integration and separating the real and the 
imaginary parts we have: 
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. In the two-dimensional overdamped regime {r <^ e < '-f^) 



the expression ( |148| ) gives the following limits: 
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We point out that the expression ( p.48|) has been obtained without any 
hmitation on frequency. Nevertheless we have made the assumption rDg^ ^ 
1 as in section 6 and it turns out that, while this condition doesn't restrict 
our results for the AL and the DOS contributions over the full range of 
frequency, temperature and impurity concentration, this is not so for the 
MT contribution. In fact, as already mentioned in section 6.4, in the ultra- 
clean (or non-local) limit, when Tr > l/\/e, the assumption rDg^ ^ 1 is 
violated for the MT contribution and the results obtained are not valid there. 
Nevertheless one can see that this non-local situation can be realized only in 
the clean case (Tr ^ 1) for temperatures in the range l/(Tr)^ <^ e <^ 1. 
We suppose Tr ~ 1, as in the case of HTS, so we postpone the discussion of 
the non-local limit up to section 11. 



9.5 Discussion 

We will now analyse each fluctuation contribution separately and then dis- 
cuss their interplay in Re [crx(c(j)]. Because of the large number of parameters 
entering the expressions we restrict our consideration to the c-axis compo- 
nent of the conductivity tensor in the region of 2D fluctuations (above the 
Lawrence-Doniach crossover temperature). In purpose to make the discus- 
sion more transparent we present the asymptotics of the results obtained in 
the Table 4. The in-plane component will be discussed in the end of this 
section. 

The AL contribution describes the fluctuation condensate response to the 
applied electromagnetic field. The currents associated with it can be treated 
as the precursor phenomenon of the screening currents in the superconducting 
phase. Above T^ the binding energy of virtual Cooper pairs gives rise to 
a pseudo-gap of the order of T — Tc, so it is not surprising that the AL 
contribution decreases with the further increase of u^. Actually o;^^ ~ 
T — Tcis the only relevant scale for a^^: its frequency dependence doesn't 



13 



experimentally such decrease was observed recently on the in-plane conductivity mea- 



surements in the vicinity of Tc on YBCO samples [196| 
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contain T, r<^ and r. The independence from the latter is due to the fact that 
elastic impurities do not represent obstacles for the motion of Cooper pairs. 
The interaction of the electromagnetic wave with the fluctuation Cooper 
pairs resembles, in some way, the anomalous skin-effect where its reflection 
is determined by the interaction with the free electron system. 

The anomalous MT contribution also involves the formation of fluctuation 
Cooper pairs, this time on self-intersecting trajectories. Being the contribu- 
tion related with the Cooper pairs electric charge transfer it doesn't depend 
on the elastic scattering time but it turns out to be extremely sensitive to 
the phase-breaking mechanisms. So two characteristic scales turn out to be 
relevant in its frequency dependence: T — Tcand t^^. In the case of HTS, 
where r^^ has been estimated as at least O.lTc for temperatures up to 5-j-lOi^ 
above T^, the MT contribution is overdamped, is determined by the value of 

and is almost temperature independent. 

The DOS contribution to Re [o"(ci;)] is quite different with respect to those 
above. At low frequencies {u ^ r^^) the lack of electron states at the Fermi 
level leads to an opposite effect in comparison with the AL and MT contribu- 
tions: Re a^'^^{u;) turns out to be negative and this means the increase of 
the surface impedance, or, in other words, the decrease of reflectance. Nev- 
ertheless, the applied electromagnetic field affects the electron distribution 
and at high frequencies uj ~ the DOS contribution changes its sign. It 
is interesting that the DOS contribution, as a one-electron effect, depends 
on the impurity scattering in a similar manner to the normal Drude conduc- 
tivity. The decrease of Re a^'^^{u;) starts at frequencies u; ~ min{T', r"^} 
which for HTS are much higher than T — Tc and r~'^. 

The cu-dependence of Refa^*] with the most natural choice of parameters 
(r <^ £ < T^^ ^ min{T, r^^}) is presented in Fig. 33. 

The positive AL and MT contributions are pronounced at low frequencies 
on the background of the DOS contribution which remains in this region a 
negative constant. Then at a; ~ T — Tc the former decays and the Recrj, 
remains negative up to ~ min{T, r~^}. The DOS contribution changes its 
sign at ~ r^^ and then it rapidly decreases. The following high frequency 
behavior is governed by the Drude law. So one can see that the characteristic 
pseudo-gap-like behavior in the frequency dependence of the optical conduc- 
tivity takes place in the range u E [T — Tc,t~^]. The depth of the window 
increases logarithmically with e when T tends to Tc, as shown in Fig. 34. 



107 



In case of a6-plane optical conductivity the two first positive contributions 
are not suppressed by the interlayer transparency, and exceed considerably 
the negative DOS contribution in a wide range of frequencies. Any pseudo- 
gap like behavior is therefore unlikely in (Tj|^'(u;): the reflectivity will be of 
the metallic kind. The comparison between (|142|) and (|143|) shows that the 
compensation of the two contributions could only take place at c^o ~ T/ Ine 
which is out of the range of validity of the AL contribution. 

Let us compare now the results of these calculations with the experiments 
available. The recent measurements p|, |^ of the c-axis reflectivity spectra, in 
FIR region on YBa2Cu40s single crystals, show the response of a poor metal 
with the additional contributions from IR active phonon modes (which we 
do not discuss here). 

With the decrease of temperature the c-axis optical conductivity decreases 
showing around 180K a transition from a Drude-like to a pseudogap-like 
behavior. The value of gap is ~ 180cm~^ and it seems weakly dependent 
on temperature. The decrease of temperature does the gap structure more 
pronounced without any abrupt change at the superconducting transition 
temperature Tc = 80K (see Fig. 35). 

Such experimentally observed behavior of the optical conductivity is in 
qualitative agreement with our results. The suppression of the density of 
states due to the superconducting fluctuations in the vicinity of Tcleads to a 
decrease of reflectivity in the range of frequencies up to ~ r~^. The magni- 
tude of this depression slowly (logarithmically) increases with the decrease of 
temperature but at the edge of the transition it reaches some saturation be- 
cause of the crossover to the 3D fluctuation regime (where instead of ln(l/e) 
one has ln(l/r) — y/e, see (|143| )). So no singularity is expected in the value 
of the minimum even in the first order of perturbation theory. Below 
the fluctuation behavior of (^j;) is mostly symmetrical to that one above 
Tc (see IQ) and, with the further decrease of temperature, the fluctuation 
pseudo-gap minimum in the optical conductivity smoothly transforms itself 
into the real superconducting gap, which opens very sharply in HTS. We 
stress that the temperature independence of the pseudo-gap threshold ap- 
pears naturally in the theory: it is determined by uq ~ r~^(see ( |145|) , ( |146| ) 
and Fig. 36). Comparing Fig. 33(b) and Fig. 36, one can easily see that 
the threshold doesn't move when the temperature changes but varies when 
the inverse of the scattering rate changes. As far as far as the numerical 
value of ujq is concerned, assuming T^r = 0.35 (which is the value for the 
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scattering rate of the sample used in the experiment under consideration [^j, 
that is also in the experimental range of the inverse of the scattering rate 
Tt ^ 0.3 ^ 0.7| |197| , |198| ), one can see that the pseudo-gap threshold is pre- 
dicted to be of the order of 200cm~^, in quantitative agreement with the 
experimental data 

The outlined theory is, strictly speaking, valid only in the vicinity of the 
critical temperature, where e -C 1. Nevertheless the logarithmic dependence 
on e of the result obtained gives grounds to believe that qualitatively the 
theory can be valid up to e = ln(T/Tc) ~ 1, i.e. for temperatures up to 
200K in the experiment discussed. So the theory is again in agreement with 
the experimental value of temperature 180K up to which the pseudo-gap is 
observable. 

In conclusion we have calculated the optical conductivity tensor for lay- 
ered superconductors. A pseudo-gap-like minimum of its c-axis component in 
a wide range of frequencies for temperatures in the vicinity of is found. It is 
due to the fluctuation density of states renormalization which can be treated 
as the opening of a fluctuation pseudo-gap. These result are quahtatively, 
and in some aspects quantitatively, in agreement with recent experiments on 
YBa2Cu40s samples. Further experiments, with more anisotropic samples 
like BSSCO single crystals, would be useful, because the effect should be 
more pronounced. 
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10 Thermoelectric power above the supercon- 
ducting transition 

10.1 Introduction 

Thermoelectric effects are difficult both to calculate and measure if com- 
pared to electrical transport properties. At the heart of the problem lies the 
fact that thermoelectric coefficients in metals are the small resultant of two 
opposing currents which almost completely cancel. In calculating thermo- 
electric power one finds that electrons above the Fermi level Ep carry a heat 
current that is nearly the negative of that carried by electrons below Ep. 
In the model of a monovalent metal in which band structure and scattering 
probabilities are symmetric about Ep, this cancellation would be exact; in a 
real metal small asymmetry survives. In calculating this small effect one can- 
not with impunity ignore any possible correction or renormalization merely 
on the grounds that it has been shown to be negligible in calculations of 
the conductivity alone. This small effect is proportional to the electron-hole 
asymmetry factor fas, which is defined as the ratio of the difference between 
numbers of electrons and holes to the total number of particles. Apart from 
this factor, to estimate the thermoelectric coefficient, one should consider 
also the characteristic energy involved in thermoelectric transport, e*. The 
thermoelectric coefficient may then be obtained from electrical conductiv- 
ity cr, if transport coefficients are defined through the electric current fiowing 
on the system as a response on electric field and temperature gradient: 

Je = aE + i?Vr, i9 - {e*/eT)fasa 

For non-interacting electrons e* T and fas ~ T/Ep, and therefore ■j?o ~ 
{T/eEp)ao- Because of their compensated nature, thermoelectric effects are 
very sensitive to the characteristics of the electronic spectrum, the presence 
of impurities and peculiarities of scattering mechanisms. The inclusion of 
many-body effects, such as electron-phonon renormalization, multi-phonon 
scattering, a drag effect, adds even more complexity to the problem of cal- 
culating the thermoelectric power. Among such effects, there is also the 
influence of thermodynamical fluctuations on thermoelectric transport in su- 
perconductor above the critical temperature. This problem has been at- 
tracting the attention of theoreticians for more than twenty years, since the 
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paper of Maki |^ appeared. After this pioneering work about ten contra- 
dictory papers appeared and up to now not only the magnitude, but also the 
temperature dependence and even the sign of the fluctuation correction to 
thermoelectric coefficient ("i?) are under discussion. 

The main question which should be answered is whether or not the correc- 
tion to ■}} has the same temperature singularity in the vicinity of the critical 
temperature as the correction to electrical conductivity cr. In the paper 
of Maki |^2[ the only logarithmically divergent contribution was predicted in 
the two-dimensional case and its sign was found to be opposite to the sign 
of the normal state thermoelectric coefficient "i^o- Later on, in a number of 
papers |199| , pO0| , |201| , p02|| it was claimed that the temperature singular- 
ity of the fluctuation correction to is the same as it is for a (oc (T — Tc)~^ 
in 2D). Finally, Reizer and Sergeev have recently revisited the problem 
using both quantum kinetic equation and linear response methods and have 
flnally shown that, in the important case of an isotropic electronic spectrum, 
strongly divergent contributions [^, |199| , pOCI|] are canceled out for any di- 
mensionally, with the flnal result having the same logarithmic singularity as 
found by Maki, but of opposite sign. We should emphasize that in all papers 
cited above only the Aslamazov-Larkin (AL) contribution was taken into ac- 
count, because the anomalous Maki-Thompson (MT) term was shown to be 
absent in the case of thermoelectric transport. It was mentioned that 
the reason for discrepancy between different authors lies in the difficulties 
connected with the introduction of many-body effects in the heat-current 
operator. As a matter of fact, the incorrect evaluation of interaction cor- 
rections to heat-current operator can produce erroneously large terms which 
are really canceled out within the consistent procedure. Due to this strong 
cancellation the AL term turns out to be less singular compared with the 
corresponding correction to conductivity Pi |. 

On the other hand, in every case where the main AL and MT fluctuation 
corrections are suppressed for some reason, the contribution connected with 
fluctuation renormalization of the one-electron density of states (DOS) can 
become important. In this section we show that the analogous situation 
also occurs in the case of the thermoelectric coefficient |E03, fflH. In what 



follows we study the DOS contribution to the thermoelectric coefficient of 
superconductors with an arbitrary impurity concentration above T^. We will 
be mostly interested in 2D case, but the generalization to the case of layered 
superconductor will be done at the end. We show that, although the DOS 



111 



term has the same temperature dependence as the AL contribution [Q, it 
turns out to be the leading fluctuation contribution in both the clean and 
dirty cases due to its specific dependence on electron mean free path. 



10.2 DOS contribution to thermoelectric coefficient 

We introduce the coefficient in the framework of linear response theory as: 

where Q^^'^^'^{uj) is the Fourier representation of the retarded correlation func- 
tion of two current operators 

g(eh)R(^ - X') = -e{t - t'){{[j''{X), r{X')])). (152) 

Here and are heat-current and electric current operators in Heisenberg 
representation, X = (r,t) and ((•■■)) represents both thermodynamical av- 
eraging and averaging over random impurity positions. We use units with 
h = c = k-B = 1- 

The correlator in the diagrammatic technique is represented by the bubble 
with two exact electron Green's functions and two external field vertices, the 
first, ev, associated with the electric current operator and the second one, 
|(e„ + en+u)^, associated with the heat current operator (e„ = 7iT{2n + 
1) is fermionic Matsubara frequency and v = dC,{p)/dp with ^ being the 
quasiparticle energy). Working to the first order of perturbation theory in 
Cooper interaction and averaging over impurity configurations one finds the 
ten diagrams presented in Fig. 37. 

As usual, the solid lines represent the single-quasiparticle normal state 
Green's function averaged over impurities. The shaded objects are the vertex 
impurity renormalization A which, neglecting the q -dependence (see Section 
6) that is unimportant here, are given by 

A(q = 0, €„,€„/)= A(e„,e„/) = -p-—^ (153) 

Finally the wavy line represents the fiuctuation propagator L(q, flk). 

The first diagram describes the AL contribution to thermoelectric coef- 



ficient and was calculated in |Q with the electron-hole asymmetry factor 
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taken into account in the fluctuation propagator. Diagrams 2-4 represent 



the Maki-Thompson contribution. As was mentioned in Refs. ||43|, neither 
anomalous nor regular parts of these diagrams contribute to in any order 
of electron-hole asymmetry. In what follows we will discuss the contribution 
from diagrams 5-10 which describes the correction to due to fluctuation 
renormalization of one-electron density of states. 

Let us start from the contribution of the diagrams 5 and 6. We have 
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A^(e„, -e„)G^ (p, e„) G (q - p, -e„) G (p, e„+i.) + (154) 



{d^)v^ X 



-tn^v)G'^ (p, tn^v) G (q - p, -tn^v) G (p, 



(We use the shorthand notation (rfq) = d^q/{2T[)^, where D is dimension- 
ality). Evaluating Eq. ( |154| ) one naturally obtains a vanishing result if 
electron-hole asymmetry is not taking into account. The first possible source 



of this factor is contained in the fluctuation propagator and was used in [44 



for the AL diagram. Our calculations show that for the DOS contribution this 
correction to the fluctuation propagator results in non-singular correction to 
d in 2D case and can be neglected. That's why we ignored such correction 
in Eq. (|64D. Another source of electron-hole asymmetry is connected with 
expansion of energy- dependent functions in power of ^/ Ep near Fermi level 
performing p -integration in Eq. ( |154| ) [Ep is the Fermi energy). In the case 
under discussion we have to perform such expansion of the following product: 



Ar(Ov2(O=iV(0)v2(0)+e 
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J 5=0 



Only second term in Eq. ( |155| ) contributes to thermoelectric coefficient. 
Contribution of diagrams 7 and 8 is described by the expression: 
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(156) 
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+ A^(e„+^, -e„+^)G^ (p, e^+u) G (p, e„) J {dk)G^ (k, e„+,,) G (q - k, -e„+j,) 

The calculation of this expression is analogous to Eq. ( |154| ), but the ex- 
pansion of iV(^)v(^) should be performed either in p or in k momentum 
integration. Diagrams 9-10 do not give any singular contribution to ther- 
moelectric coefficient due to the vector character of external vertices and 
as a result an additional factor appears after p-integration. The same 
conclusion concerns the MT-like diagram 3-4. 

Performing integration over ^ we find the contribution of the important 
diagrams 5-8 in the form 
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9(iV(0v2(0) 



y'(dq)L(q,0)(Si + S2 + E3) (157) 



In the last equation we have taken into account that near only the term 
with = is important and we have separated sums over semi-infinite 
(] — oo, —u — 1], [0, oo[) and finite {[—i^, —1] ) intervals : 




^2 = ! E (2e„ + co.r - I) (158) 

S3 = (1 + u^t) E (2en + 

n=—u 

El and S2 are associated with diagram 5-6, while S3 with diagram 7-8. 
Calculating the sums ( |158| ) we are interested in terms which are linear in 
external frequency u,^. The sum Si turns out to be an analytical function 
of uJi, and it is enough to expand it in the Taylor series after analytical 
continuation u^, — —iu. The last two sums over finite intervals require more 
attention because of their nontrivial tUiy-dependence and before analytical 
continuation they have to be calculated rigorously. As a result: 

sf = ^; s^ = -^; = (159) 
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Finally, we perform integration over q. In ID case the total contribution 
associated with density of states renormalization takes the form 
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In 2D case: 
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(162) 



The generalization of this result to the important case of layered supercon- 
ductor is straightforward. One has to replace ln(l/e) — * ln[2/ (v^+ \/e + r)] 
and to multiply Eq. ( |162| ) by l/pps, where s is the interlayer distance. In 
the limiting case of 3D superconductor (r ^ e) both the AL and the 
DOS contributions are not singular. 



10.3 Discussion 



Comparing Eq. ( |162|) with the results of for the AL contribution, we 
conclude, that in both limiting cases of clean and dirty systems the decrease of 
i) due to fluctuation the DOS renormalization dominates the thermoelectric 
transport due to the AL process. Really, the total relative correction to 
thermoelectric coefficient in the case of 2D superconducting film of thickness 
s can be written in the form: 
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where the first term in square brackets corresponds to the DOS contribution 
( 162 ) and the second term describes the AL contribution from (G^ is 
the Debye temperature). Assuming ln(6j)/Tc) ~ 2 one finds that the DOS 
contribution dominates the AL one for any value of impurity concentration: 
K* as a function of Tr has a minimum at Tr ^ 0.3 and even at this point 
the DOS term is twice larger. In both limiting cases Tr -C 1 and Tr ^ 1 
this difference strongly increases. In one- dimensional case 
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for Tr < 1 



and the DOS correction again turns out to be dominant. 

The temperature and impurity concentration dependencies of fluctuation 
corrections to i) in important 2D case can be evaluated through a simple 
qualitative consideration. The thermoelectric coefficient may be estimated 
through the electrical conductivity a as -i? ~ {e* /eT)fa_^a, where e* is the char- 
acteristic energy involved in thermoelectric transport and /as is the electron- 
hole asymmetry factor, which is defined as the ratio of the difference between 
numbers of electrons and holes to the total number of particles. Conductiv- 
ity can be estimated as a ~ e^A/'r*/m, where N , r* and m are the density, 
lifetime and mass of charge (and heat) carriers, respectively. In the case 
of the AL contribution the heat carriers are nonequlibrium Cooper pairs 



T-Tc 



and characteristic 
Thus in 2D case 



8(T-Tc)- 

jSj^- Oiie can easily get that the 



with energy e* T — Tc and density A/" 
time, given by Ginzburg-Landau time r* 

A^?AL _ _ Tj/(eTe) /asAa^L _ ^f^^ 1^ 

fluctuation correction due to the AL process is less singular (logarithmic in 
2D case) with respect to the corresponding correction to conductivity and 
does not depend on impurity scattering |^ . 

The analogous consideration of the single-particle DOS contribution (e* ~ 
T, r* ~ r) evidently results in the estimate 'd ~ e/asTcr In j^^jt which coin- 
cides with ( |162| ) in clean case. The dirty case is more sophisticated because 
the fluctuation density of states renormalization strongly depends on the 
character of the electronic motion, especially in the case of diffusive motion 
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. The same density of states redistribution in the vicinity of Fermi level 
directly enters into the rigorous expression for ■(9 and it is not enough to 
write the fluctuation Cooper pair density A^.p. but is necessary to take into 
account some convolution with 6Nq{E). This is what was actually done in 
the previous calculations. 

Experimentally, although the Seebeck coefficient 5* = —d/a is probably 
the easiest to measure among thermal transport coefficients, the comparison 
between experiment and theory is complicated by the fact that S cannot be 
calculated directly; it is rather a composite quantity of electrical conductivity 
and thermoelectric coefficient. As both d and a have corrections due to 
superconducting fluctuations, total correction to the Seebeck coefficient is 
given by 

AS^S„(^-^) (165) 

Both these contributions provide a positive correction A^9, thus resulting in 
the decrease of the absolute value of S at the edge of superconducting transi- 
tion (At^/t^o < 0). As for fluctuation correction to conductivity Acr/cxo > 0, 
we see from Eq. (|165|) that thermodynamical fluctuations above always 
reduce the overall Seebeck coefficient as temperature approaches Tc. So the 
very sharp maximum in the Seebeck coefficient experimentally observed in 
few papers ||17| , |18|, p02|| seems to be unrelated to the fluctuation effects within 
our simple model even leaving aside the question about the experimental reli- 
ability of these observations. This conclusion is supported by recent analysis 
of temperature dependence of thermoelectric coefficient close to transition in 
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11 DOS fluctuations on NMR characteristics 
in HTS 



11.1 Introduction 

In this section we discuss the contribution of superconducting fluctuations to 
the spin susceptibihty Xs and the NMR relaxation rate 1/Ti. We base on the 



work of M.Randeria and A.Varlamov [461 were the effect of fluctuations on 



Xs was examined first time. The fiuctuation contribution to NMR relaxation 



rate 1/Ti has been previously studied in the set of papers ||204| , p05| , |206 



focusing mainly on the most singular contribution, the anomalous Maki- 



Thompson (MT) term in the dirty limit Tr <^ 1. Nevertheless in [pOT] , ^ 
the problem was reexamined and the important role of the DOS contribu- 
tion was underlined. The matter of fact that in conditions of the unusually 
strong pair breaking in high materials, one might expect the main MT 
contribution to be suppressed and 1/Ti to be dominated by the less singular 
DOS contributions. 



It has recently been suggested [ p08|| that dynamic P09|| pairing correla- 



tions beyond the perturbative weak coupling regime are responsible for the 
spin gap anomalies ||210|| observed well above Tc in the underdoped cuprates. 



The analysis presented here, which only treats static fluctuations very close 



to Tc, nevertheless constitutes the flrst correction to Fermi liquid be- 
havior, to order m.ax.{l/ EpT^Td Ep), arising from pairing correlations above 

The main results of this section, valid for e -C 1, can be summarized as 
follows: 

(1) Fluctuations lead to a suppression of the spin susceptibility Xs , due 
to the combined effect of the reduction of the single particle density of states 
(DOS) arising from self energy contributions, and of the regular part of the 
Maki-Thompson (MT) process. 

(2) "Cooperon" impurity interference terms, involving impurity ladders 
in the particle-particle channel, are crucial for the Xs suppression in the dirty 
limit. 

(3) The processes which dominate the results in (1) and (2) above have 
usually been ignored in fluctuation calculations (conductivity, 1/Ti, etc.). 
The spin susceptibility is unusual in that the Aslamazov-Larkin, and the 
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anomalous MT terms, which usually dominate, are absent. 

(4) For weak pair-breaking (I/t,^ -C Tc) , we find an enhancement of 1 /TiT 
coming from the anomalous MT term. We recover known results | |205|| in the 
dirty limit, and extend these to arbitrary impurity scattering. 

(5) In the clean limit (T^r ^ 1) we find a different asymptotic behavior 
of 1/TiT depending on whether one has T^r is greater or smaller than 

(6) Finally, strong dephasing suppresses the anomalous MT contribution, 
and 1/TiT is then dominated by the less singular DOS and the regular MT 
terms. These contributions lead to a suppression of spectral weight and a 
decrease in 1/TiT. 

11.2 Definitions 

We begin with the dynamic susceptibility x*|^(k, uj) = x+_(k, iuj,y uj+iO'^) 
where 

XH-(k,cu.) = dre'-^^iif {S4k,T)S4-k,0)))) (166) 

with the Bose frequency = 27rz/T. Here S± are the spin raising and 
lowering operators, T is the time ordering operator, and the brackets denote 
thermal and impurity averaging in the usual way. The uniform, static spin 
susceptibility is given by Xs = X+^(k 0, u; = 0) and the NMR relaxation 
rate by 

-L = hm - / (dk) QxfliKuj) (167) 

III i^^O 00 J 

where A is a positive constant involving the gyromagnetic ratio. 

For non-interacting electrons x+-(k, i^i/) is determined by the loop dia- 
gram presented in Fig. 38. Simple calculations lead to the well known results 
for T < Ep: Xs = ^(0) (Pauh susceptibility) and (l/TiT)° = Att[N{0)]'^ 
(Korringa relaxation), where N{0) is the DOS at the Fermi level. We shall 
present the fluctuation contributions in a dimensionless form by normalizing 
with the above results. 

To leading order in ma:x{l / Eft,Tc/ Ep) the fluctuation contributions to 
X+- are given by the diagrams shown in Fig. 38. The diagrams are con- 
structed from fermion lines, fluctuation propagators (denoted by the wavy 
lines) and impurity vertex corrections (represented by shaded objects). It is 
important to note that the two fermion lines attached to the external ver- 
tex have opposite spin labels (up and down) for x+-- Consequently, the 
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Aslamazov-Larkin diagram (1) does not exist since one cannot consistently 
assign a spin label to the fermion line marked with a '?' for spin-singlet 
pairing. 

The next set of diagrams to consider is the Maki-Thompson (MT) di- 
agram (2a and 2b), and the MT with the Cooperon impurity corrections 
(3) and (4). While the MT diagrams for x+- appear to be identical to the 
well known MT diagrams for conductivity, there is an important difference 
in topology which arises from the spin structure. It is easy to see, by draw- 
ing the fluctuation propagator explicitly as a ladder of attractive interaction 
lines (diagram (2b)), that the MT diagram is a non-planar graph with a sin- 
gle fermion loop. In contrast the MT graph for conductivity is planar and 
has two fermion loops. The number of loops, of course, affects the sign of 
the contribution. 

The diagrams (5) and (6) represent the effect of fluctuations on the one- 
particle self energy, leading to a decrease in the DOS. The DOS diagrams 
(7) and (8) include impurity vertex corrections. (Note that these have only 
a single impurity scattering line as additional impurity scattering, in the 
form of a ladder, has a vanishing effect.) Finally (9) and (10) are the DOS 
diagrams with the Cooperon impurity corrections. 

The fermion lines represent the one electron Green function G(p, a;„) = 
(ie„ — ^(p))""^, where e„ = -|- sgn(en)/2r describes the self-energy effects 
of impurity scattering. The momentum relaxation rate 1/r -C Ep , however, 
TcT is arbitrary. For simplicity, we will discuss the two-dimensional (2D) case 
for the most part, indicating at the end how the results are modified for the 
layered system, and the 2D to 3D crossover. 

Pairing fluctuations above Tc are described in the usual way by means of 
the vertex part of electron-electron interaction in the Cooper channel [^2|, ^], 
or the fluctuation propagator (shown in Fig. 38 by the shaded wavy line) 
L(q, Vtfj) where the arguments refer to the total momentum and frequency 
of the pair. The full (q, VL^) -dependence of L is important far from Tc; see 



[ ^ . We restrict our attention here to the vicinity of transition and thus it 
suffices to focus on long wavelength, static {VL^ = 0) fluctuations. In this 
regime we have 

(q, n^ = 0) = -N{0) [e + r/^g'] , (168) 
We now turn to vertex corrections due to impurity scattering. First, it 
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can be shown Q that the external vertices do not need to be renormalized 
by impurity lines. Next, the (three-legged) impurity vertex A(q, —en, ^n+i) is 
defined as the sum of impurity ladders dressing the bare vertex consisting of 
two fermion lines, with frequencies — e„ and en+u = ^n + ^v, and a fiuctuation 
propagator with momentum and frequency (q, uo^). Having in mind to discuss 
in this section the non-local limit (/ ^ Cafe(^) ) and to show how this modifies 
the MT contribution (the same situation takes place in conductivity, but we 
skipped this in the view of cumbersome calculations of the section 6) we use 
below the purely 2D presentation of the vertex, calculated for an arbitrary 
value of r~^Dg^. The result differs from (|59| ) (but naturally coincides with 
the latter in the local limit) 

A(q, -e„, e.+.) = | 1 - , Q(^»^»+-) \ . (159) 



where Q{x) is the Heaviside step function. 

Finally, the Cooperon C(q, —en,en+y) is defined as the sum of impurity 
ladders in the particle-particle channel where — e„ and En+v are the frequencies 
of the two particle lines and q is their total momentum. It is given by 



C(ct) ~Cn) ^n+u) 



27tN{0)t 



~To ; 1 , -Tk 2 ^n^n+u 

T |2e„ + uj^l + Dq^ 



(170) 



Note that the diffusion pole involves the sum of the momenta and the differ- 
ence of frequencies. 



11.3 Spin Susceptibility 

We note that, when the external frequency and momentum can be set to zero 
at the outset, as is the case for Xs, there is no anomalous MT piece (which 
as we shall see below is the most singular contribution to l/Ti). The MT 
diagram (2) then yields a result which is identical to the sum of the DOS 
diagrams (5) and (6). We have: 

Xs5 + Xse ^_Z!L\n(l/e) (m) 
Xs^ Ef 



^"^The external vertex correction equals zero for the case of Xs since the external fre- 
quency LUi, = 0. For 1/Ti, the large external momentum k leads to a suppression of this 
correction by the parameter 1/EpT. 
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and this result is valid for any impurity concentration. 

We first discuss the clean limit, where the fluctuation contribution is given 
by Xs = Xs2 + Xs5 + Xse! S'll other diagrams are negligible for T^t ^ 1. The 
final result is 

^ = --^ln(l/£) forT,r»l. (172) 

In the dirty limit (TcT <^ 1), the DOS diagrams (5) and (6), together with 
the regular part of MT diagram (2), yield the same result (|172|) of the order 
Tc/Ep. One can see, that this contribution is a negligible in comparison 
with the expected for the dirty case dominant one, of the order 0{l/EpT). 
The thorough study of all diagrams shows that the important graphs in dirty 
case are those with the Cooperon impurity corrections MT (3) and (4), and 
the DOS ones (9) and (10). This is the first example known to us where 
the Cooperons, which play a central role in the weak localization theory, 
give the leading order result in the study of superconducting fluctuations. 
Diagrams (4) and (5) give the one half the final result given below; diagrams 
(9) and (10) provide the other half. The total fluctuation susceptibility = 
Xs3 + XsA + Xs9 + Xsio, is 

- '^^^^^ ^ ln(l/£) forT,r<l. (173) 



It is tempting to physically understand the negative sign of the fluctuation 
contribution to the spin susceptibility in eqns. ( |172| ) and ( |1 72| ) as arising from 



a suppression of the DOS at the Fermi level. But one must keep in mind that 
only the contribution of diagrams (5) and (6) can strictly be interpreted in 
this manner; the MT graphs and the coherent impurity scattering described 
by the Cooperons do not permit such a simple interpretation. 



11.4 Relaxation Rate 

The calculation of the fluctuation contribution to 1/Ti requires rather more 
care than Xs because of the subtleties of analytic continuation. We define the 
local susceptibility Kiujy) = /((ik)x+_(k, cj,^). The "anomalous" MT contri- 
bution, denoted by a subscript (an), comes from that part of the Matsubara 
sum which involves summation over the interval —{uj^/2t[T) = —u < n < — 1. 
The reason this piece dominates, and has to be treated separately, is that 
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S'(an) has a singular q- dependence as we shall see. The remaining terms in S, 
i.e., those with n < —u and n > 0, are called the "regular" MT contribution. 

We evaluate S'(an) using standard contour integration techniques and then 
make the analytic continuation iuj^ ^ uo + iO^ to obtain 



lim l^Ki'll^iuj) = J (rfq)L(q, 0)/C(q), (174) 

dz 1 



. COsh2(z/4Tr) (^^I2q2 _ _ ,)2 _ (^y^/2g2 _ + ,)2 _ 

(175) 

where I = vt is the mean free path. 

The first simple limiting case for ( |175| ) is Iq <^ 1, for which /C(q) = 
27r/Dg^. Since ( |1 74| ) involves a q-integration, we need to check when the 
above approximation is justified, the characteristic g- values being determined 
from the fluctuation propagator L. In the dirty limit, we have Dg^ ~ eT^ 
, thus leading to /^g^ ~ eTcT ^1. In the clean case, on the other hand, 
v^gVTc ~ eT^ and Pq^ ~ siT^ry < 1 only when K T^r < 

For the above conditions (either TcT <^ 1 or 1 ^ TcT ^ y/^) we obtain 
the singular MT contribution / d'^q [(Dg^ + l/r(^)(e + ?72g^)] where we have 
introduced the pair breaking rate 1/t<^ as an infrared cutoff. We define the 
dimensionless pair-breaking parameter 7^ = ?72/Dr^ ^ 1; in the dirty limit 
lifi = T^/^TcT^ while for the clean case 7^ = 7({3)/16tt'^T^tt^ . The "bare" 
transition temperature T^q is shifted by the pair breaking, so that e = eo+T^^; 
with Eq = (T — Tco)/TcQ, and we obtain the final result 

The other limiting case of interest is the "ultra-clean limit" when the 
characteristic g- values satisfy /g ^ 1. This is obtained when T^r ^ 3> 
1. From ( |175|) we then find /C(q) = 41n(/g)/t'g, which leads to 

(l/TiT)fl vr^ 1 

(i/rir)o yi4C 3 Ef ve 



We note that in all cases the anomalous MT contribution leads to an 
enhancement of the NMR relaxation rate over the normal state Korringa 
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value. In particular, the superconducting fluctuations above Tc have the 
opposite sign to the effect for T <^Tc (where 1/Ti drops exponentially with 
T). One might argue that the enhancement of 1/TiT is a precursor to the 
coherence peak just below T^. Although the physics of the Hebel-Slichter 
peak (pile-up of the DOS just above gap edge and coherence factors) appears 
to be quite different from that embodied in the MT process, we note that 



both effects are suppressed by strong inelastic scattering [^, pl2| , |213 |. 

We now discuss the DOS and the regular MT contributions which are 
important when strong dephasing suppresses the anomalous MT contribution 
discussed above. The local susceptibility arising from diagrams (5) and (6) 
can be easily evaluated. The other remaining contribution is from the regular 
part of the MT diagram. This corresponds to terms with n < —u and n > 
in the Matsubara sum. It can be shown that this regular contribution is 
exactly one half of the total DOS contribution from diagrams (5) and (6). 
All other diagrams either vanish (as is the case for graphs (7) and (8)) or 
contribute at higher order in 1/EpT (this applies to the graphs with the 
Cooperon corrections). The final results are given by 



(l/TiT)o E, 



ln(l/e) (178) 



for TcT ^ 1, and 



(V^i^)^ _ 21C(3) 1 

Ji/jyif - ^ ^^^^/"^ ^^^^^ 

for T^r < 1. 

The negative sign of the result indicates a suppression of low energy 
spectral weight as in the Xs calculation, however in contrast to Xs we note 
that the same graphs (2), (5) and (6) dominate in both the clean and dirty 
limits. 

It is straightforward to extend the above analysis to layered systems by 
making the following replacement in the 2D results given above: 
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11.5 Discussion 



The negative DOS contribution to the NMR relaxation rate is evident from 
the Korringa formula and it sign seems very natural while the sign of the 
positive Maki-Thompson contribution can generate a questions about its 
physical origin. It is why we present here the qualitative consideration of the 
result obtained above from the microscopic theory in the spirit of the section 
6.2. 

In the case of nuclear magnetic relaxation rate calculations the electron 
interaction with nuclei with spin flip is considered. If one would try to imag- 
ine the AL process of this type he were in troubles, because the electron-nuclei 
scattering with the spin-flip evidently transforms the initial singlet state of 
the fluctuation Cooper pair in the triplet-one, what is forbidden in the scheme 
discussed. So the formally discovered absence of the AL contribution to re- 
laxation rate is evident enough. 

The negative density of states contribution in ^ has the same explanation 
given above for the conductivity: the number of normal electrons decreases 
at the Fermi level and as a result the relaxation rate diminishes with respect 
to the Korringa law. 

The positive Maki-Thompson contribution can be treated in terms of the 
pairing on the self-intersecting trajectories like this was done with that one in 
conductivity (see section 6.2). This consideration [ 214 | clears up the situation 
with the sign and explains why the MT type of pairing is possible when the 
AL one is forbidden. Nevertheless the principle of the electron pairing on the 
self-intersecting trajectory in the case of the NMR relaxation rate has to be 
changed considerably with respect to the case of conductivity. 

Let us consider a self-intersecting trajectory and the motion of the elec- 
tron along it with fixed spin orientation (e.g. "spin up"). If, after passing 
of the full turn, the electron interacts with the nucleus and changes its spin 
state and momentum to the opposite value it can pass again the previous 
trajectory moving in the opposite direction (see Fig. 39). The interaction of 
the electron with itself on the previous stage of the motion is possible due to 
the retarded character of the Cooper interaction and such pairing process, 
in contrast to the AL one, turns out to be an effective mechanism of the 
relaxation near T^. This purely quantum process opens a new mechanism of 
spin relaxation, so contributes positively to the relaxation rate 

In conclusion of this section we want to say several words about the ex- 
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perimental difficulties in the observation of fluctuation effects in NMR. Much 
of the existing work on fluctuations as probed by NMR has been restricted to 
small particles (zero-dimensional limit) of conventional superconductors; see 
pi5|| for a review. There has been resurgence of interest in superconducting 
fluctuations since the high Tc cuprates show large effects above due to their 
short coherence length and layered, quasi-two-dimensional structure. How- 
ever, in order to extract the "fluctuation contributions" from experiments 
one needs to know the normal state backgrounds, which in a conventional 
metal would simply be the Pauli susceptibility for Xs and the Korringa law 
for 1/TiT. The problem for the high Tc materials is that the backgrounds 
themselves have nontrivial temperature dependencies [0 above T^. for exam- 
ple, the non-Korringa relaxation for the Cu-site in YBa2Cu307_5, and the 
spin-gap behavior [piO|| with dxs/dT > and the O and Y 1/TiT ~ Xs(T). 
The attempt to overpass these difficulties was done in |^ which we discuss 
in the next section. 
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12 D versus S pairing scenario: pro and contra 
from the fluctuation phenomena analysis 

12.1 Introduction 

One of the hottest debate in sohd state physics is the question of the pairing 
state symmetry in HTS. The determination of the order parameter symmetry 
is the crucial first step in the identification of the pairing mechanism in HTS 
and subsequent development of a microscopic theory of high temperature 
superconductivity. Recently much of the attention of both theoreticians and 
experimentalists has been focused on a pairing state with ci— wave symme- 
try, although many other pairing symmetries are also possible. The challenge 
for investigators is to derive tests that are capable to distinguish the differ- 
ent possible pairing states and to make the determination of pairing state 
symmetry. The direct way to study the anisotropy of the phase of the or- 
der parameter is investigation of phase coherence of Josephson and tunnel 
junctions incorporating HTS. Recent phase coherence experiments support 
a scenario of superconducting pairing state with d^2_y2 symmetry in HTS 



216| , |217|| . Nevertheless the studies of the amplitude of the order parameter 
may be also fruitful and the simplest experiments are studies of fluctuation 
effects above superconducting transition temperature. In view of the fact 
that the theory of fluctuation effects in different physical properties of super- 
conductor is now well established, the main question of the present section 
is how the interpretation of experimental data on fluctuation effects in HTS 
alters if one assumes d— rather than s— pairing state and which experiment 
on fluctuation effect would be a test for a pairing state in HTS. 

In the theoretical studies of fluctuation effects in c?— wave superconductors 
two models were used. First is a model spherically symmetric superconductor 
with the bare electron-electron interaction with the interaction constant 

V(k,k') =^?oP2(kk') (181) 

(P2 is Legendre polynomial) pi8| ]. The second model considers tight-binding 
electrons with nearest-neighbor transfer with underlying 2D square lattice 
mmpq. in this case 



V(k, k') = (7o(cos /c^^a — cos A;j^a)(cos /i;^a — cos A;' a) (182) 
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with a being a lattice constant. Both models lead essentially to the same 
results. Further we will use the second model which provides a more simple 
calculations. The main difference between s— and d— wave superconductors 
lies in the different manifestation of non-magnetic impurities, which are pair- 
breaking for the last. Further, in order to avoid cumbersome equations, we 
restrict ourselves to establish only dependencies of fluctuation contributions 
on temperature and impurity scattering thus ignoring numerical factors of 
order unity. 



12.2 The fluctuation propagator 

In the case of rf— wave superconductor the fluctuation propagator depends 
upon directions of k and k', but using the ansatz: 

-^kk'(^) = L{q){cos k^a — cos kyo) {cos k'^a — cosk'yO) (183) 

(q = k — k'), one can easily solve the Dyson equation for a quantity L[c\). 
One has: 

^(q)"' = go-^ - P{^,ook), (184) 

f <Pk 

P(q,cjfc) = j ——;^{coskxa-coskyafG{\^,Un)G{(l-\^,Uk-Un) 

Here 

G(k, cij„)~^ = iu;„+isgncij„/2r — , C,k = —t{coskxa+coskya)—fi (185) 

{t and /i are the transfer integral between the nearest-neighbor sites and 
the chemical potential, respectively). Note, that in contrast to s— wave 
pairing state, the impurity renormalization of quantity P(q, uj^) is simply 
absent ||218|, |205|] due to dependence of bare interaction on the momentum 



directions. 

For temperatures close to transition point, one can solve (|184D for small 
q and cjfc. In this case 

^q-k ~ 6 + ^ , ^ = at{qx sin k^a + qy sin kyO) (186) 

Then one can replace the k by v / {2Ti)d^{...) , where v is 2D density of states 
and (...) means averaging over angles of k taken under condition = 0. 
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After straightforward calculations one has: 

n,_, 1^-1 A^l^^ \n 



"2(g) 



Q;i47rT2 



Here we defined: 



CKi = ((cos/cajO — coskyaY) ~ 1 
'^2(?) = {{coskxa — coskyofrf) Pi a^t^q^ 

The divergence of L(q, 0;^) at q = 0, cufc = determines the critical temper- 
ature. Finally: 



-1 



47rr V2 47rTT/ Q;i47rT2'^ V2 AttTt 



(188) 



where 



is a shifted by impurities critical temperature and 

Teo^^expf —] (190) 

(where 7 is the Euler constant) is the hypothetical critical temperature in 
the absence of impurities. One can see that usual non-magnetic impurities 
result in the reducing of critical temperature for (i— wave superconductor in 
the same manner as magnetic impurities in the case of s— wave supercon- 
ductor. It is interesting to note, that taking into account typical for HTS 
values of parameter Tt ^ 1 one obtains a huge suppression of Tcq. Thus 
Tco/Tc = 2.86, 1.85, and 1.43 for Tt = 0.5, 1, and 2, respectively, result- 
ing in embarrassing prediction of critical temperature in pure HTS materials 
at room temperature. It seems that such prediction rules out completely 
the d— wave scenario. Nevertheless, one should realize that HTS materials 
are obtained by the isovalent substitutions from the parent compounds and 
structural disorder here is an intrinsic feature. Thus considering the "pure" 
limit for these materials has no sense. 
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12.3 d.c. conductivity 



The contribution of AL process to in-plane conductivity was calculated by 
Yip [ pl8|| . Current-current response function, ignoring factors of order unity, 



is: 



Q 



AL _ 2 
XX ~ ^ 



(fq 



'9P(q, Uk) 



dqx 



As a result: 



AL 



1 1 

2 ^ AttTt 



e 

de 



(191) 



(192) 



Here t/;— function represents the pair-breaking effect of impurities which re- 
duces the magnitude of AL conductivity in the dirty case with respect to 
s— wave superconductor. As impurity vertex corrections do not contribute, 
the anomalous MT contribution is simply absent ||218|| . Other terms (DOS 
and regular MT) give: 



- e^T/^L(q,0)TE/c^6-(^)'G(k,c.„)3G(-k, 



-UJr. 



a, 



DOS 



-— In- 

d e 



(193) 



In contrast to what happens in the s— wave case, the DOS contribution does 
not depend on r. It is easy to extend ( p,92| ) and ( |193D to the case of layered 
superconductors. In this case one can calculate also out-of-plane component 
of conductivity tensor. The results are similar to those found above: both 
relevant contributions have the same temperature dependence as for s— wave 
case, but other dependence on r. Namely, AL contribution is suppressed for 
dirty case, while the DOS contribution does not depend upon r. 



12.4 NQR-NMR relaxation 

For s— wave the main contribution to spin relaxation originates from MT pro- 
cess, if the latter is not suppressed by strong intrinsic pair-breaking. Instead, 
the hierarchy of fluctuation contributions changes essentially if one consider 
the (i— wave superconductor. As both AL and MT processes are absent, the 
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negative DOS term becomes the only present. Corresponding results for the 
d-wave superconductor with quasi-two-dimensional spectrum are 



X 



(0) (l/TiT)(o) t eV2 + (£ + r)i/2 



By comparing ( |194|) with the appropriate result for s-wave case from P6 
one can see that the type of pairing does not effect on the magnitude and 
temperature dependence of the DOS contribution. 

Finally, one can conclude that the only essential difference between s— 
and d— pairing states in context of fluctuation theory consists in the ab- 
sence of anomalous MT process in the latter. Therefore, studies based on 
measurements of thermodynamical character do not give the possibility to 
distinguish the type of pairing on the basis of fluctuation effects experi- 
ments. In-plane and out-of-plane conductivities, the experimental data in 
zero field do not manifest the signs of MT term. Nevertheless the measure- 
ments of fluctuation conductivity cannot provide reliable tests for possible 
d— or s— pairing because the high values of determine strong pair-breaking 
(at least due to the electron-phonon scattering the expected minimal value 
of r^^ ~ i'(^)^)- Thus the MT process is ineffective even in the case of s— 
pairing scenario. Additionally, as the MT contribution to conductivity tem- 
perature dependence is similar to AL one, they can hardly be distinguished. 
The additional comparison with s— wave case function gives the numeri- 
cal factor of order unity, which cannot be tested experimentally. Therefore, 
the only physical property related to superconducting fluctuations which is 
due to the MT contribution is the NMR relaxation rate, in which AL process 



does not contribute at all [46 



Within the context of the s— wave scenario both the positive MT sin- 



gular contribution ||205|| and the negative DOS |^ (independent on phase- 
breaking) fluctuation renormalization exist. The concurrence of these two 
effects should be observable in the relaxation measurements. In the case of 
c?— pairing, vice versa, the MT anomalous process in accordance with the 
consideration presented above does not contribute at all and in the NMR re- 
laxation rate above Tc a monotonous decrease with respect to Korringa law 
has to be observed. Consequently, from point of view of fluctuation theory, 
the sign of the correction to NMR relaxation rate above could be a test for 
the symmetry of the order parameter. Further, as the form of background 



131 



(normal state) relaxation is unknown, the better choice for experimental veri- 
fication of the existence of the MT contribution is the study of relative change 
of 1/Ti in the vicinity of transition induced by external magnetic field. The 
main results for the fiuctuation contributions to the relaxation rate in exter- 
nal field in the case of s— wave superconductor are the following. The DOS 
contribution has the same form as calculated above. The MT contribution 
is given by 

(l/riT)(o) ^ E^re- 7, + (7^ + r) V2 ^ ^ 

While the MT contribution is very sensitive to the presence of pair-breaking, 
the simplest way to discriminate DOS and MT contributions is to apply an 
external magnetic field, a relative small value being expected practically to 
suppress the MT contribution. 

12.5 Experiment and discussion 

As it is clear from discussion above the most appropriate physical property 
those measurement would be used as a test for a pairing state symmetry 
in HTS is NQR-NMR relaxation rate. Thus below we present the result 
of recent experimental work performed in order to get a clear evidence 
of the role of superconducting fiuctuations and meantime to achieve some 
conclusion about the pairing mechanism in HTS. 

At first, we discuss experimental details which are important for interpre- 
tation of a such kind of experiments in HTS. From the recovery of the ^^Cu 
signal after RF saturation in both type of NQR and NMR experiments one 
arrives at 1/Ti = 2W given by 

2W = — J (/i+(t)/i_(0))e-*^«*rft (196) 
where F is the gyromagnetic ratio of ^^Cu nucleus. In (|196|) h± are the 



components of the field at the nuclear site transverse with respect to the c 
axis both for NQR where the quantization axis is the Z one of the electric field 
gradient tensor as well as for NMR when the external field is along the c axis 
itself. ujR, the resonance frequency, is ci;_r(0) = 31 MHz in NQR (zero field) 
and ujr{H) = 67 MHz in NMR for H = 5.9 T. In the following this difference 
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will be neglected. The fictitious field h can be related to the electron spin 
operators S± through the electron-nucleus Hamiltonian and the relaxation 
rate can formally be written in terms of a generalized susceptibility. One can 
write 

k 

where Ak is a term involving the square of the Fourier transform of the ef- 
fective field h, which can be averaged over the Brillouin zone. In a Fermi 
gas picture in (|197|) one can introduce, in the limit ujr — > 0, the static 
spin susceptibility x°(0,0) and the density of states A^(0), namely W cx 
T[x°(0, 0) / {1 — a)]hN (0) (a is a Stoner-like enhancement factor) and the Ko- 
rringa law 1/Ti oc T is thus recovered. Having to discuss only the effect of the 
external field HonW around T+, in the assumption that the field does not 
change appreciably the electron-nucleus Hamiltonian (as it is proved by the 
equality W{0) = W{H) for T ^ T^) we will simply write W/T oc Xab = X-i 
where x is then the /c-integrated, cj/j contribution. The fact that in 
YBCO the Korringa law is not obeyed above Tc for ^^Cu NQR-NMR relax- 
ation rates (most likely because of correlation effects) should not invalidate 
the comparison of M^(O) to W{H) in the relatively narrow temperature range 
of 10 K above T^.. 

The ^^Cu NQR and NMR relaxation measurements have been carried 
out in an oriented powder of YBa2Cu306.96. From the ^^Cu NMR line-width 
(FWHI~ 40 kHz) the spread in the direction of the c axis was estimated 
within 1 — 2 degrees. The superconducting transition was estimated Tc = 90.5 
K in zero field and T^. = 87.5 K in the field of 5.9 T used for NMR relaxation. 
In ^^Cu NQR measurements the recovery of the amplitude s{t) of the echo 
signal at the time t after complete saturation of the ibl/2 —>■ ±3/2 transition 
was confirmed of exponential character, thus directly yielding the relaxation 
rate. The temperature was measured with precision better than 25 mK and 
the long term stabilization during the measurements was within 100 mK. 
In the presence of the magnetic field the sample was aligned with c || H. 
The echo signal for the central transition was used to monitor the recovery 
of the nuclear magnetization after fast inversion of the ±1/2 populations. 
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From the solution of the master equations one derives for the recovery law 
y{t) = 0.9exp{—12Wt) + 0.1exp{—2Wt) This law was observed to be very 
well obeyed and the relaxation rate was extracted. The experimental error 
in the evaluation of W was estimated well within 5 %. In Fig. 40 the 
experimental results for 21^(0) and 2W{H) around are reported. It is 
noted that W{0) = W{H) for T>T, + 15 K while the effect of the field, 
namely a decrease of the relaxation rate, is present in a temperature range 
where paraconductivity and anomalies in the c axis transport are observed. 

Let us discuss the interpretation of the experimental results in terms of the 
contributions related to superconducting fluctuations. The relative decrease 
of W around Tc{H) and Tc(0) induced by the field is about 20 %. We 
can evaluate if this decrease is quantitatively consistent with the picture of 
superconducting fluctuations. First one should estimate the strength of pair- 
breaking. By using the reasonable values, r ~ 10~^^ s and ^ 2 x 10^^^ s, 
one finds that the pair-breaking parameter 7^ is about 0.15. The effective 
anisotropy parameter r in YBCO can be estimated around 0.1. Finally one 
can observe that in zero field the relative MT contribution is larger than the 
absolute value of the DOS contribution by a factor 1.5, thus providing the 
positive fluctuation correction to W. The effect of the field on these two 
contributions can be expected as follows. In the case of strong pair-breaking 
7B > {£B,r} {eb = e + /3/2, 7^ = 7^ + f3/2, l3 = 2B/H,2iO)) there are 



two different regimes for the fluctuation corrections in magnetic field [34 



Low-field regime {(3 <^ e) corresponds to a decrease quadratic in (3 of the 
fluctuation correction. In the high- field regime {e <^ f3) one can use lowest 
Landau level approximation. In our experiment e ~ 0.05, while f3 is about 
0.2. Since relative corrections to W coincide with relative corrections to 
conductivity p4i, one easily finds 



MT 
B 



Direct calculations according to these equations show that the MT contri- 
bution is much more affected by the magnetic field. In fact from (|136|) the 
modification in W^'^^ induced by the field is small, of the order of 15 %. On 
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the contrary, according to ( |199| ) the field reduces the MT term to 1/4 of its 
zero field value. In the case of d-wave pairing symmetry, as it was already 
noted above, the MT contribution is absent and the applied magnetic field 
results only in the slight reducing of the DOS contribution in accordance 
with the equation similar with ( |136| ). If the decrease of W by magnetic field 
is due to the reduction of the MT term, then one deduces that zero-field total 
fiuctuation correction to NMR relaxation rate is positive and is within 10 % 
of background, due to the partial cancellation of MT and DOS contribution. 
In a field of 6 T the total fiuctuation correction becomes negative with an 
absolute value about 15 % of the background. The lack of detailed informa- 
tion on the normal-state NMR relaxation rate in HTS does not allow one to 
achieve more quantitative estimates. However, it should be remarked that 
the observation of the decrease of in a magnetic field cannot be accounted 
for in the case of c/-wave pairing. Finally we would like to emphasize the 



following. It has been recently pointed out by Miiller [p21|| that two type of 
condensates, with different symmetry but the same transition temperature, 
can exist in oxide superconductors. It is conceivable that also the spectrum 
of the fiuctuations of the order parameter above Tc could refiect both com- 
ponents, if present. Since the effect of the magnetic field discussed in our 
paper works only on the component of s symmetry, our conclusion is not in 
contrast with the experimental evidences indicating (i— wave pairing and it 
could be considered a support to the hypothesis of order parameter having 
simultaneously s and d symmetry. 

Summarizing, from the accurate comparison of ^^Cu NQR and NMR re- 
laxation rates in zero field and in a field of 5.9 T around in YBa2Cu30x, 
we have provided evidence of a contribution to W related to the supercon- 
ducting fiuctuations, in a temperature range of about 10 K. The experimental 
observation of a decrease in W induced by the field is consistent with the 
hypothesis of a strong reduction of the MT contribution |^. Since the MT 
contribution does not exist in the case of a rf-wave scenario, the interpreta- 



Concerning the behaviour of W{H) there is still some controversy in literature. Results 
analogous to those presented above, were obtained by the Urbana group ]222( | and indicate 
that the magnetic field reduces the relaxation rate in accordance with our consideration. 
However a comparison of NQR and NMR measurements in the YBa2Cu40x phase yields 



an opposite behaviour of W{H), which turns out to be larger than W{0) |223, 224|. The 



reason of this discrepancy could be due either to the presence of a negligible s-component 



in the orbital pairing of the YBa2Cu40x or to the role of diamagnctic terms | 214 | 
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tion of the experimental finding is an indication in favor of the presence of 
an s symmetry component in the orbital pairing. 
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13 Conclusions 



Several comments should be made in the conclusion. We have demonstrated 
that the strong and narrow in energy scale renormalization of the one-electron 
density of states in the vicinity of the Fermi level due to the Cooper channel 
interelectron interaction manifests itself in experiments as the wide enough 
pseudo-gap-hke structures. The scale of these anomalies, as we have seen 
above, is different for various phenomena {eV = ttT for tunnel conductance, 
uj ~ r^^ for optical conductivity, T ~ Tc in NMR and conductivity measure- 
ments). 

The results presented above are based on the Fermi liquid approach which 
is formally expressed through the presence of a small parameter of the theory 
Gi2 In ^ ~ In ^ <^ 1 in obtained results. Moving along the phase dia- 
gram of HTS from the metal region (overdoped or optimally doped samples) 
to poor metals (underdoped compounds) one can see that the small parame- 
ter of the perturbation theory grows {Ep 0) causing the effects discussed 
to be more pronounced. Nevertheless, analysing the rapid growth of the 
normal state anomalies with the decrease of the oxygen content below the 
optimal doping concentration one can notice that it strongly overcomes our 
theoretical prediction. We can attribute this discrepancy to the simplicity 
of the Fermi surface model supposed above to be isotropic in the ab-plane. 
The ARPES study of HTS shows the presence of the strong anisotropy of 
the Fermi surface of such type and even the existence of two characteristic 
energy scales Ep ^ O.SeV and A ^ O.OleV (extended saddle point of the 
spectrum). So one can suppose that as the oxygen concentration decreases 
below the optimal one, by some reasons the massive part of the Fermi surface 
is " obliterated" and the crossover in properties related with the special role of 
"slow" electrons of extended saddle points takes place. Formally in this case 
the large value Ep ^ O.SeV in the denominator of the Ginzburg-Levanyuk 
parameter has to be substituted by small A ^ O.OleV making rapidly the 
perturbation approach to be unapplicable. 

The existence of the nonequilibrium Cooper pairs in the normal metal 
phase of HTS (the state with < >^ 0, < ^ >= 0, < >= 0) resembles 
the state of preformed Cooper pairs in the underdoped phase |p26| , |227|| (< 

>7^ 0, < \l/ >7^ 0, < >= 0). Both of them are determined by the 
presence of interelectron interaction and it would be interesting to study the 
plausible "condensation" of the fluctuation pairs in the preformed ones in the 
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unknown land "hie sunt leones" where Gi2 ~ 1. Such crossover quahtatively 



was discussed by M.Randeria [ P25|| , but its systematic study requires the 
formulation of the appropriate model. 

In conclusion, we have presented here what we hope to be a compre- 
hensive overview of the theoretical and experimental facts which allow us 
to attribute an important role in the behaviour of HTS (the metal part 
of its phase diagram) to the fluctuation theory, and especially, to the fre- 
quently neglected DOS contribution. We did show that a large number of 
non-trivial experimental observed anomalies of HTS normal state properties 
can be qualitatively, and often quantitatively, described under this model. 
Of course, as a rule, this is not the only explanation which can be given 
for these behaviours. However, in our opinion, alternative approaches often 
lack the internal coherence and self consistency of the picture based on the 
fluctuation theory. In particular, it is important that all the experimental 
facts reviewed here can be explained by the fluctuation theory alone, without 
need for any other ad hoc assumption or phenomenological parameters, and 
that the values of all physical parameters extracted from the fits are always 
in a good agreement with independent measurements. Even more important 
is the fact that within the same approach one can explain so many different 
properties of HTS, in the wide temperature range from tens of degrees below 
to tens of degree above the critical temperature. It is very unlikely that these 
circumstances are purely fortunate. 

Therefore we believe that, although almost for all experiments discussed 
the alternative explanations can be proposed, no other model, besides the 
fluctuation theory, has the same appeal embracing in a single microscopic 
view the wide spectrum of unusual properties which have been discussed in 
this review. 
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15 Appendix A: Calculation of the impurity 

vertex \[q^uJi,uj2] 



Let us demonstrate the details of calculations of the impurity vertex A(g, c^i, UJ2) 
which appears in the particle-particle channel as the result of the averaging 
over the impurities configuration. These calculations can be done in the 
frameworks of the Abrikosov-Gorkov approach to the diagrammatic descrip- 



tion of the alloys [p2S 



In the assumption of relatively small impurity concentration ppl ^ 1 
(what in practice means the mean free path up to tens of interatomic dis- 
tances) and for the spectra with dimensionality D > 1 one can neglect the 
contribution of the diagrams with intersecting impurities lines. In these con- 
ditions the renormalized vertex A(g, cui, to'2) is determined by the graphical 
equation of the ladder type (see Fig. 41): 

Analytically the last equation can be written in the form 

A"^(g, ^1, 0^2) = 1 - o Arrn\ ^1' ^2), (199) 

27rA(Ujr 

where 

V(q,u:u^2) = j(dp}G{p + q,uli)G(-p,ul2)= (200) 

= iv(o){ r 7 ^^^^^ = 



= 27rAr(0)e(-a;i^2)( p . . ] w,.^ ^| )f.s.- 

Here 6(x) is the Heaviside step-function and 

Aaq,p)\\p\=PF = [ap+ <i) - a-m\p\=PF- m) 

Now one has to reduce the formal averaging of the general expression 
over the Fermi surface (( ...) f.s.), to the particular one of the expression 



(|20 1|) . For an arbitrary g*it is possible to accomplish this program for some 
simple spectra types only. Below we will show such exact solution for the 
important for HTS case of the 2D isotropic spectrum, but let us start before 
from the analysis of the practically important calculation of X{q,uj 1,002) for 
small momenta in the case of an arbitrary spectrum. 
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The angular averaging of ( |20CI| ) over the Fermi surface can be carried out 
in general form for values of 



A^iq,p}\\p\=p^<^\0i-U2\, (202) 

what in the case of corrugated cylinder spectrum ( p4D means q\\ <^ mm{^^f^, l^^}. 
Expanding the denominator of (|20CI|) over the appropriate small parameter 
one can find: 

V{q,uJi,uj2) = — — — — '-^—r^ e(-u;iu;2) (203) 



27rA^(0) ' ' Icui— a;2| \iUi — lu2 

and 

= \ :■ (204) 

1^1 - + TWT^ Q{-UJlUJ2) 

It is easy to see that the restriction ( |202| ) is not too severe and almost 
always can serve as a reasonable approximation for temperatures near (the 
exclusion is the very clean case, when the non-local fluctuation effects take 
place (see section 11)). For instance, from the calculations of section 6 one 
can see that not far from the effective propagator momenta are determined 
by \q\eff ~ [^ab^{T)]~^ = CbVe < while in Green functions the q- 
dependence becomes important for much larger momenta q ~ min{^afe,/~^} 
(what is equivalent to the condition (P02|) ). 



The last average in (|204|) can be easily calculated for some particular 
types of spectra. For example in the cases of 2D and 3D isotropic spectra it 
is expressed in terms of the diffusion coefficient : 

{{^mP)\\p\=P.y)F.s.D = t-'Bm' = (205) 
while in the case of the quasi- two- dimensional electron motion (0) expression 



( |201| ) may be presented by means of the action of the diffusion operator D 



on the momentum q: 

{{Aaq,p)\\p-i=p,?)F.s. = ^(4q' + 4J'sin2(g,s/2)) = T~'t)q\ (206) 
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The vertex A(g, cji, C(J2) in this case can be written as 



and namely this expression was used through the review. In the dirty case 
(Tr -C 1) this expression is reduced to (|59D . 

Finally we discuss the exact calculation of the vertex X{q,u!i,u!2) in the 
case of 2D electron spectrum which turns out to be necessary for the con- 
sideration of the Maki- Thompson contribution in the extra-clean limit / 3> 
Cab'iT). In this case the q -dependence of the vertices in the domain of small 
momenta becomes of the same importance than that one which comes from 
the propagator, but because of the large value of I the condition ( p02| ) is 
violated. This is why we express (|201|) in the exact for 2D case form 



^^{q,P)\\p\=pf = [^{P + Q) - ^i-P)]\\p\=PF- = VFq cos yD, (208) 

which is valid for any vpq ^ Ep , and average ( |200|) without any expansion. 
The angular average in this case is reduced to the calculation of the integral 

^(^'"^'"^^ = w7q Jo ^^^^7riliEMEE#209) 

r- vpq 

which can be easily carried out. Really, the integral in (|209| ) is calculated by 
means of the substitution 2; = e**^: 

2^ 1 _i-sign{a) ^^^^^ 



lo 271 COS (p. — ia vT+o^ 
what gives for V{q,uJi,ui2) : 

-vi- \ 27riV(0)e(-^i^2) 

V{q, ^1, U2) = I ■ (211) 

\JVpq^ + [uJi - UJ2Y 

The proper exact expression for A(g, cui, u;2) can be written as 

m -1, = f 1 - , ^^-^y^'^ ^ ] ' . (212) 

which was used for the study of non-local limit in section 11. One can see that 
this expression can be reduced to ( |204|) - ( |205| ) in the case olvpq <C |a)i — a;2|. 
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16 Appendix B: Calculation of the fluctua- 
tion propagator L{q,uj^) 

In this section we discuss the calculation of the fluctuation propagator L{q, u^) 
which is nothing else that the two-particle Green function of electrons inter- 
acting in the Cooper channel. Graphically it is presented by the sum of the 
diagrams with two entrances and two exits. In the BCS theory it is demon- 
strated that the ladder type diagrams with the accuracy of In-^(^) turn 
out to be of first importance in the calculation of L(g, a;^). This statement 
can be applied for any type of weak interaction with characteristic energy 
Eq S> Tc in the range of temperatures Tc<T <^ Sq, so in purpose to calculate 
the fluctuation propagator L(g, cu^) one has to solve the ladder type Dyson 
equation presented in the Fig. 42: 

Analytically this equation can be written in the form 

L-\q,uj^)=g-^-Ii{q,uj^) (213) 

where g is the effective constant of the electron-electron interaction in the 
Cooper channel and the polarization operator Il{q. lo^) {uj^ > 0) is determined 
as the two one-electron Green's functions correlator averaged over impurities 
positions, which can be expressed in terms of the introduced above functions 
X{q, ujn+i^, uj^n) and V{q, ujn+^, uj^n) (see Appendix 1) : 



/d'^p 
J^^G{p + q, UJn+n)G{-p, UJ-n) (^14) 

= T^\{q,Un+,j,,UJ-n)V{q,Un+t,,UJ-n), 

One can easily calculate L{q,ujp) for the most important case of small 
momenta A^(g,p)||p|=p^ <^ \oJn+fjt — ^-n\ ~ max{ T,t~^}. Using the proper 
expressions for A(q', and V{q,u!n+iu,i^-n) it is possible to flnd: 



1 - 



{{^aq,p)\\p\=ppf)F.s. 



n(g>,) = 2.TAr(o)Ee(--n+.-_„)- , - 

= AnTN(0)[j: (215) 
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-m{q,P)\\P\=P.)")F.S. E 



The calculation of sums can be accomplished in terms of the logarithmic 
derivative of F-function ipi^x) what leads to the explicit expression for Il{q, uj^) 



rn(g,c.,) = + + + (216) 



Ar(0) ^' ^'2 47rT 27rr' ^'2 AnT' 

{i^C{q,P)\\p[=PF?)F.s. 



.^4 inTr' ^^2' At:Tt ^2\ 



The term proportional to ((A^(g',p)||p^=p^)^)i?.5. may be expressed in terms 
of the r] -coefficient of Ginzburg-Landau theory or diffusion operator D : 

^ n(,,c.,)=^(i + -;^ + -5-)-^(i + -;^)-r)g- (217) 



N{Q) ^' ^^2 AnT 2t:T' ^^2 AnT 

The definition of the critical temperature as the temperature Tc at which 
the pole of L(0, 0, Tc) takes place 

L'\0, 0, = g-' - n(0, 0, Te) = (218) 

permits us to express the fluctuation propagator in terms of the reduced 
temperature e = In(^) : 

""'^ ^ ~W)HB + ^(i + 4) - ^(i) + ^'''^ 

One has to remember that this expression was carried out in the assumption 
of the small momenta A^{q, p)\\p^=pj^ <S \ojn+n — <^-n| ~ max{T, r"-*^}, so 
the range of its applicability is restricted by the Ginzburg-Landau region of 
temperatures e = ln(;^) <^ 1, where the integrated functions in diagrammatic 
expressions have the singularities at small momenta. 

Nevertheless it is possible to generahze the last expression for arbitrary 
momenta. The first hint appears in the dirty case when rj ~ '0 (|) 
one can suppose that the term rj(f is nothing else as the first term of the 
expansion of the function V(^ + ^^^'i^)- One can find the confirmation of 
this hypothesis in the fact that the analytically continued expansion of the 
propagator 
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which is vahd for the region A^(g, -C max{T, r ^} and u <^ T, is 

nothing else as the Time Dependent Ginzburg-Landau (TDGL) equation's 



fundamental solution. As it is well known [p29|| the Maki-De Gennes equa- 




tion, which in the case discussed has the form 



M|)+v.(i + i;:!i±^^il^)-^*(i)|* = o (221) 

serves as the generalization of TDGL equation . So it is natural to suppose 
that the most general form of the propagator, valid for arbitrary momenta 
and frequencies, has the form 



L(g,tu„) = ^ ^- (222) 

(The Larmour frequency Ql appeared in ( p22| ) as electron eigen-energy of 



the Landau state in magnetic field). 

The second confirmation of the correctness of the hypothesis proposed 
is the direct calculation of the polarization operator from the formula ( pi4| ) 
for the case of a 2D spectrum, when it is possible to carry out exactly the 
angular average over the Fermi surface in V{q, ujn+^i, ^^-n) ■ Really, using the 
definition of \{q,u!n+^,uj^n) we can rewrite the polarization operator in the 
form: 

n(g,.;,)=T^-- — (223) 

uJr^ [y[q,UJn+f,,i^-n)\ - 27rN{0)T 

In Appendix 1 for 2D spectrum case it was found: 



^/^ N 27tN(0)Q(-UJiUJ2) , 

^(^' ^1' ^2) = . , (224) 
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what gives for the polarization operator 



n(g, uj,) = N{0) Y: I (225) 

j^/t -r 2 47rT ^ 4ttTt J ^ Ww^T^ A-kTt 

For max{ ^^Jtrp2 , ^ 1 the summation can be substituted by integration 
and one can find (for simphcity we put Tt ^ oo): 

L{q, u^) = ^ X i ^ (226) 

""^^^ + In (i + + v/(i + + j0fi) - m 

Taking into account that il){x ^ 1) — >■ Inx one can see that ( |226| ), in the 
hmit of maxj^jy^, ^ 1, coincides with 
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17 Glossary 



AL term ( paraconductivity, SaAi) is Aslamazov-Larkin contribution; 
ARPES is angular resolved photo-emission; 

Aa is a vector-potential; 

a is the interatomic distance; 

a — ^J^^ is the coefficient of the Ginzburg-Landau theory; 

ai — {{coskxa — coskyaY) ^ 1 is the coefficient of the Ginzburg-Landau 
theory for d-pairing; 

^2(0') = {{coskxa — coskyaYrf) ?s a^t^q^ is the kinetic energy part of the 
Ginzburg-Landau theory for d-pairing; 

^ _ 4r;2Dg dimeusiouless kinetic energy in the fluctuation propagator; 

Ba{q,ui^,ujiy) is the integrated three Green functions block of AL contri- 
bution; 

P = B/[2Tc \dBc2/dT\Tc] = 4:rieB = is reduced magnetic field; 

B is magnetic field; 

Bc2 is upper critical magnetic field at zero temperature; 
r is the gyromagnetic ratio; 
7 is the Euler constant; 

= ■;^;r^ ^ 8^ is the phase-breaking rate related with r^; 

7b = Ti^ + /3/2 is the phase-breaking rate related with in the presence 
of a magnetic field; 

D is diffusion coefficient; D ~ ^zl EeL jg j^g generalization from 

dirty to clean case; 

is the diffusion tensor; 

DOS term {Saoos) is the density of states fluctuation contribution; 

D is the space dimensionality; 

Di{z) = 2 In + y {z + r) is the function of complex variable; 

D2{z) = —\Jz{z + r) is the function of complex variable; 
ADi {z) = Di {z) — Di {e) is the function of complex variable; 
AD2 {z) = D2 {z) — D2 {e) is the function of complex variable; 
d is the fllms thickness; 
Ep is the Fermi energy; 

£'(p) = p^/(2m), is the kinetic energy of 2D free electron; 
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Eo{T) is the energy scale at which the DOS renormahzation occurs; 
e'^^^ ~ T - Te; 

e — ln(T/Tc) is reduced temperature; 

£b — s + (3/2 is reduced temperature in magnetic field; 

({x) is the Riemann zeta function; 

fas ~ T/Ep, the symmetry factor entering in the thermo-epf; 
f{e) is a universal function entering in 2D AL paraconductivity; 
Gi(^j:i) = Gi(£))(47rTT) is Ginzburg-Levanyuk parameter; 
G{V — 0, T) is the tunnel junction zero-bias conductance; 
G{V) is the differential tunnel conductance; 

G„(0) is the background value of the Ohmic conductance supposed to be 
bias independent; 

SG{V) = G{V) - Gn{0); 

G{p, ujn) = is the single quasiparticle normal state Green's func- 

tion; 

Iqpiy) is tunneling current; 
!{dc^) = !d^q/{2^)^- 

J d^q = J d'^qJ^^J'^gdQz is momentum space integral transformation for a 
layered superconductor; 

Iai3{q,uj^,u!i,) is the integrated four Green functions block entering in the 
MT contribution; 

I are the unit vectors along the axes; 

J is a hopping integral describing the Josephson interaction between lay- 
ers; 

and are heat-current and electric current operators in Heisenberg 
representation; 

Vd = (I + - ^ (I) - 4^^'(|) *he coefficient of the 

gradient term of /^-dimensional Ginzburg-Landau theory; rj = rj2] 

HTS - high temperature superconductors; 

Hc{0) is the zero-temperature thermodynamical critical field; 

h± are the components of the magnetic field at the nuclear site transverse 
with respect to the c axis both for NQR as well as for NMR; 

K{uJ,y) = J{dk)x+-(k,uj^); 

^ _ ii±i2 aj-g ii^Q functions of the impurity concentration entering in 
the DOS contribution in the case of an arbitrary impurity concentration; 
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k {uJ,T,T~^) is the function determinating the frequency and temperature 
dependence of the pseudogap in optical conductivity; 

Lx,y,z are the sample dimensions in appropriate directions; 
I is the intralayer mean-free path; 
Lt = is the diffusion length; 

L(g, uj^) is the fluctuation propagator; in the absence of the magnetic field 
and in the vicinity of Tc, it has the form : L~^{q, uj^) = —N{0) e + ip (^^ + + aqj — ip 

-^kk'(^) = I/(q)(cos fca^a — cos /c^a) (cos /c^a — cosk'yO) is the fluctuation 
propagator in the case of d-pairing; 

A(g, a;„, a;„') is the impurity vertex in the Cooper channel; 
/X is a chemical potential; 

MT term, anomalous and regular {SauT = 5(7^^(''^»)+(5(7^^("")), is the 
Maki-Thompson contribution; 

Nc.p. is the concentration of Cooper pairs; 

SMs.i. is the concentration of interfering Cooper pairs; 

•^i^^ ~ ^-^P one-electron concentration in 2D case; 

A^(2) = ^ is the density of states for 2D electron gas; 

A^l(O), A^r(O) are densities of states at the Fermi levels in each of elec- 
trodes of tunnel junction in the absence of interaction; 

SNjf{E, e) is the fluctuation correction to the density of states of the 2D 
electron gas; 

Hi is the impurity concentration; 

^(p) — E{p) + Jcos{pzs) — Ep is the spectrum of corrugated cylinder 
type; 



^0 is the superconductor coherence length at zero temperature: ^o,d 
iI^St^ clean case 

and ^o,d = It i^^ the dirty case; 

1 /2 

^(T) = ^0 i^ T^T ) temperature dependent coherence length of the 

Ginzburg-Landau 
theory; 

^ab is the in-plane BCS coherence length of layered superconductor; 
is the out of plane BCS coherence length of lavcrcd superconductor; 

V{q,LJi,U2) is the two one-electron Green's function loop integrated over 
internal electron momentum; 
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n(g, Ufj^) is the polarization operator defined as the two one-electron Green's 
function correlator averaged over impurities positions; 

p = (p,Pz) is a vector of momentum space; 

p = ipxiVy) is a two-dimensional intralayer wave-vector; 

Q(ch)Rj^^^ is the Fourier representation of the retarded correlation function 
of two current operators; 

g(eh)R(^_^/) ^ -Q{t-t'){{\j\X),J%X')\)) : X = (r,t) and ((■■■)) 
represents both thermodynamical averaging and averaging over random im- 
purity positions; 

Qapi^^u) is the electromagnetic response operator; 

Rn is the Ohmic resistance for unit 

Pab{T) and Pc{T) the components of the resistivity tensor in a6-plane and 
along c-axis; 

r = Ari2J'^ /v% — 4^1(0)/s^ is the Lawrence-Doniach anisotropy parame- 
ter; 

ri,r2 are the functions of impurity concentration accounting for the con- 
tributions of DOS diagrams; 

S — —-d/a is the Seebeck coefficient; 

crQ/3 is the conductivity tensor; 

(y^{uj) is Drudc conductivity; 

Tco is the BCS value of critical temperature; 

T* is the critical temperature reduced by the effect of fluctuations; 

5Tc is the fluctuation shift of the critical temperature; 

To and Ti are respectively the mean energy and the half width of the 

energy spread of the resonant defects referred to Fermi level;. 
T is the period of electron Bloch oscillations; 

1/Ti = 2Vr = ^ J{h+{t)h-{0))e~'''^^^dt is the NMR relaxation rate; 
VT is a temperature gradient; 

f^^ — D^~^ ~ TqI ~ T — Tc is the inverse of the time necessary for the 
electron to diffuse over a distance equal to the coherence length ^{T); 



t^^ ^ Vf^~^ ~ yTc{T — Tc) the same value for the ballistic motion; 
t is the transfer integral between the nearest-neighbor sites in the theory 
of d-pairing; 

T(i){e) is the one-electron phase-breaking time; 

Tgl = 7rh/8kB{T — Tc) is the characteristic time of the Timc-Dcpcndcnt 
Ginzburg-Landau Theory; plays the role of a fluctuation Cooper pair lifetime 
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in the vicinity of Tc] 

T is the quasiparticle scattering time; 

Thop is the characteristic time of anisotropic diffusion from one layer to 
the neighboring one; 

G(x) is the Heaviside step function; 
O/) is the Debye temperature; 

is the thermoelectric coefficient; 
$0 is the elementary magnetic flux; 
V is the volume of the sample; 

V^(k, k') = go{coskxa — cos ky a) {cos k'^a — cos kyO) is the interaction po- 
tential for the case of d-pairing; 

i^aip) = ^^-is the quasiparticle velocity; vp is the Fermi velocity in 
a6-plane; 

""zip) — = —Jssm{pzs) is the electron velocity along the c-axis 

direction; 

X+!^(k, u;) = x+-(k, icUi, uj + iO^) is the dynamic susceptibility; 
Xs = ^ 0,u; = 0) is the spin susceptibility; 

^'(r) - superconducting order parameter, ^'^ = J ^'(r) exp"*'^'' ciV^ is 
its Fourier transform; 

'il^{x) and t/j^'^^x) are the digamma function and its derivatives respec- 
tively; 

is the fluctuation part of the thermodynamical potential; 
f2x, is the Larmor frequency; 
Lo^ — {2u + l)7rT are the Matsubara frequencies; 

Un — ujn[i + l/(2|(^n|''')] the Matsubara frequencies renormalized by 
the impurity scattering; 

^n±u = <^n i ujiy, is the sum of Matsubara frequencies; 

TTOJ 

Cj — — — — -— is dimensionless electromagnetic field frequency; 
ujr is the resonance frequency. 
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FIGURE CAPTIONS 

Figure 1. The temperature dependence of the averaged order parameter 
after renormahsation by fluctuations. The dash-dotted hne represents the 
BCS curve for the unperturbed transition temperature, T^. The dashed hne 
is the BCS curve for the fluctuation- renormahsed transition temperature, T*. 
The solid line is the fluctuation-renormalised average of the squared order 
parameter, (T), as given in Eq. (p]) . 

Figure 2. The theoretical prediction of Eq. (^) for the normalised cor- 
rection 6N{E) to the single-particle density of states vs energy E (measured 
in units of Tc) for a clean two-dimensional superconductor above Tc. TqI 
assumes the values 0.02Tc, 0.04Tc and 0.06Tc. In the inset the dependence of 
the energy at which 6N{E) is a maximum, Eq, on TqI is shown |^ . 

Figure 3. The theoretical prediction of Eq. (|30|) for the fluctuation- 
induced zero-bias anomaly in tunnel- junction resistance as a function of volt- 
age for reduced temperatures t = 1.05 (top curve), t = 1.08 (middle curve) 
and t = 1.12 (bottom curve). The insert shows the experimentally observed 
differential resistance as a function of voltage in an Al-I-Sn junction just 
above the transition temperature []82| . 

Figure 4. Theoretical fit (solid line) of the experimentally observed tem- 
perature dependence zero-bias conductance of the YBaCuO/Pb junction of 
Ref. m. The theory used is that of Eqs. (|30|) and (|3|) with r = 0.07 and 
Tc = 90K. The inset shows the same results in a wider region of temperature. 

Figure 5. Normalised tunneling conductance data of the BSSCO-2212/Pb 



junctions of Ref. [13|. Pseudo-gap type non-linearities are seen in the tem- 
perature range from = 87 — 89K to llOK. 

Figure 6. Theoretical fit of the G{V, T) data from the BSSCO-2212 junc- 
tions of Ref. [0 . The solid lines are the experimental data at T = 90K and 
lOOi^ for a sample where Tc = 87K. The thin lines are the theory of Eq. 
(|30D with parameters r = and Gi(^2) = 0.008. 

Figure 7. The electronic spectrum for our model of HTSC materials as 
written in Eqn. (^4]). The Fermi surface takes the form of a corrugated 
cylinder. 

Figure 8. The origin of the Aslamazov-Larkin and Maki- Thompson con- 
tributions to fluctuation conductivity in terms of electron-electron pairing, 
(a) Correlations between electrons of opposite spin and momenta moving in 
a straight line in opposite directions lead to the Aslamazov-Larkin terms, (b) 
Correlations between electrons of opposite spin moving in opposite directions 
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along a self-intersecting trajectory in real space lead to the Maki-Tliompson 
terms. 

Figure 9. Feynman diagrams for the leading-order contributions to fluc- 
tuation conductivity. Wavy lines are fluctuation propagators, thin solid lines 
with arrows are impurity-averaged normal-state Green's functions, shaded 
semicircles are vertex corrections arising from impurities, dashed lines with 
central crosses are additional impurity renormalisations and shaded rectan- 
gles are impurity ladders. Diagram 1 is the Aslamazov-Larkin term; diagrams 
2-4 are the Maki- Thompson terms; diagrams 5-10 arise from corrections to 
the normal state density of states; diagram 11 is an example of a higher-order 
contribution. 

Figure 10. The minimal self-intersecting trajectory for the Maki-Thompson 
contribution to the c-axis conductivity. A quasiparticle hops from one layer 
to neighbouring layer, diffuses within that layer, hops back to the first layer, 
and finally diffuses back to its starting point. 

Figure 11. The graphical form of the Maki-Thompson contribution show- 
ing that it is the precursor phenomenon of the Josephson effect. Starting from 
the diagram for the Josephson current for T < Tc, we expand the Gorkov 
F-functions to linear order in A as T ^ Tc, and for T > Tc replace A, A* by 
the AA* propagator which is just the fluctuation propagator L. 

Figure 12. Theoretical predictions for the zero-field resistivities p^x/ Pxx 
(dashed curves) and Pzz/p^z (solid curves) vs reduced temperature T/T^o. In 
each plot tTcq = 1 and t^Tcq = 1 (top curves), T(j,Tco = 10 (middle curves) 
and T^Tco = 100 (bottom curves). All temperatures shown are in the region 
1 < T/Tco < 1.06. Plots (a)-(c) differ in the values of the parameters r and 
Ef: (a) r(Tco) = 0.1, Ef/T^o = 300, (b) r(T,o) = 0.01, Ep/T^o = 300, (c) 
r{Tco) = 0.001, Ef/Tco = 500. 

Figure 13. Theoretical predictions |Q for zero- field resistivities Pxx/Pxx 
(dashed curves) and Pzz/ Pzz (solid curves) vs reduced temperature T/Tcq 
for the parameters r(Tco = 0.01, Ep/Tco = 300, tTco = 0.1 and t^Tcq = 1 
(top curves), r^^Tco = 10 (middle curves) and r^^Tco = 100 (bottom curves). 
Changing the value of tTco into the dirty limit has suppressed the peak in 

Pzz/ Pzz- 

Figure 14. Theoretical predictions for zero-field resistivity Pxx/Pxx 
with (top curves) and without (bottom curves) the density of states contri- 
butions. The parameters used are Ep/Tco = 300, r^pTco = 10 and r(Tco) = 0.1 
(dashed curves) or r(Tco) = 0.01 (solid curves). We see that the DOS conti- 
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butions lead to an overall increase in resistivity. 

Figure 15. Resistance R{t) vs reduced temperature t = T/Tc for the 



three BSCCO-2212 samples utilised in ||136 | for fluctuation measurements. 
Resistances are normalised to their values at T = 1.33Tc. 

Figure 16. Temperature dependence of the excess conductivity of 4 YBCO 
samples as measured in [|131|| . The solid lines are the predictions of the AL 



theory for 3D and 2D cases. The dashed line is the modified theory of |91I . 

Figure 17. The normalised excess conductivity /(e) = (16/is/e^)A(T of the 
BSCCO samples in Fig. 15 plotted against e = In (T/Tc) in a In-ln plot as 
described in ||136|| . The sohd line is the extended theory of Reggiani, Vaglio 



and Varlamov ||117|| . The dashed line is the 2D AL theory. 

Figure 18. The normalised excess conductivity /(e) = {16hs/e^)Aa for 
samples of YBCO-123 (triangles), BSSCO-2212 (squares) and BSSCO-2223 
(circles) plotted against e = In (T/Tc) in a In-ln plot as described in [PT|]. 
The dotted and solid lines are the AL theory in 3D and 2D respectively. The 
dashed hne is the extended theory of [|117| . 

Figure 19. Fit of the temperature dependence of the transverse resistance 
of an underdoped BSCCO c-axis oriented film with the results of the fluctu- 



ation theory |32|. The inset shows the details of the fit in the temperature 
range between Tc and llOK. 

Figure 20. Comparison of the curvatures of experimental curves of pd Pah 
vs temperature for a BSSCO film with the predictions of the resonant 
tunneling model of Abrikosov (p4|). The solid lines are the theoretical fits; 
the points are the data for a reduced sample (open squares), as-grown sam- 
ple (closed circles) and oxidised sample (open diamonds). These fits show 
that the experimental data cannot be fit by a theory which predicts singular 
behaviour at T = instead of at T = Tc. The dashed line is the simulated 
behavior according to the fluctuation theory for the argon annealed sample. 



Figure 21. Plots of the possible behaviour |33] of the normal-state c- 
axis resistance of a BSSCO-2212 sample after the predicted fluctuation con- 
tributions have been subtracted from the experimental data (open circles). 
The solid lines are the subtracted data for three values of the Fermi energy: 
Ep = O.SeV (bottom curve), Ep = l.OeV (middle curve) and Ep = l.25eV 
(upper curve). We see that after subtraction of the fluctuation contribution, 
there may still be a weaker temperature dependence of the normal-state re- 
sistance to explain. 

Figure 22. Experimental data showing the effect of a magnetic field par- 
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allel to the c-axis on the in-plane resistivity Pab{T) and transverse resistivity 
Pc(T) of BSCCO single crystals The insets show the zero field behaviour. 
We see that the size of the peak increases and moves to lower temperature. 

Figure 23. Theoretical prediction for the magnetic field dependence 
of resistivities p^xiT)/ p^^ (dashed curves) and PzziT)/ p^^ (solid curves). The 
parameters used are tTcq = 1, t^Tcq = 10, r{Tco = 0.01, and Ep/Tco = 300, 
with field strengths corresponding to /3(Tco) = 0,0.05 and 0.1. 

Figure 24. Transverse magnetoconductivity of YBCO thin films (o) (from 
[ |1 56|] ) and single crystals (•) (calculated from [|180|] ) as a function of the re- 
duced temperature e (after ||187|| ). The broken line is a fit to the clean-limit 
theory [|165| , |179|| including all four contributions and the solid line represents 
the sum of the orbital AL contribution and the quasiparticle magnetocon- 
ductivity estimated from the Hall effect with A = 1.7. The latter two are 
also shown separately by dash-dotted and dotted curves, respectively. 

Figure 25. Negative c-axis magnetoresistance at various temperatures in 
2212 BSCCO single crystals after 



Figure 26. Fit of the experimental c-axis magnetoresistance in the 



BSCCO-2212 samples of Ref. to the fluctuation theory predictions of 
Eqs. (99)-(102). The fltting parameters used were vp = S.l x 10® cm/s, 
r = 1.0 X 10-^^s and = 8.7 x lO'^^s. 

Figure 27. Theoretical prediction of Eq. ( |136| ) for the temperature 
dependence of the c-axis magnetoresistivity of BSSCO-2212. The fltting 



parameters are the same as those used in Fig. 26 to flt the data of Ref. [p7 

Figure 28. Measured temperature dependence of the c-axis magneto- 
conductivity of two YBCO single crystals aX B = 12T 0. The solid lines 
are flts to the fluctuation theory prediction of Eq. ( |136| ) with fltting pa- 
rameters vf = 2x Wcm/s, t{100K) = r^(lOOi^) = (4 ± 1) x lO'^^s and 
J = (215 ± 10)i^r. The inset is an enlarged view of the temperature region 
close to Tc. 

Figure 29. Fit of c-axis resistance curves of BSCCO-2212 single crystals 
in a magnetic fleld with the Josephson coupling theory of Kim and Gray 



Figure 30. Fit of the measured temperature dependence of the c-axis 
magnetoresistance of a BSSCO-2212 fllm to the fluctuation theory ||3^. The 
points are experimental data at different magnetic flelds: B = OT (circles), 
B = 0.2T = /3 = 0.003 (squares) and B = OAT = /3 = 0.006 (triangles). The 
lines are the predictions of Eqs. (99)-(102) with fltting parameters r = (5.6± 
0.6) X 10-^\ = (8.6 ± 1.4) X W-^\ Ef = (1.07 ± 0.12)eV and J = (43 ± 
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4)K. For temperatures below the resistance maximum the fit underestimates 
the resistance, the discrepancy becoming larger with increasing field. 

Figure 31. (a). Fit of experimental Pab{T) curves of a 2212 BSCCO film 
in various external magnetic fields (see inset) with the fiuctuation theory in 
the Hartree approximation after (b). Fit of experimental Pc{T) curves of 
the same film as in (a) in various external magnetic fields with the fiuctuation 



theory in the Hartree approximation after Fitting parameters (see text) 
are the same as for (a). 

Figure 32. (a), c-axis resistance vs. magnetic field in 2212 BSCCO single 
crystals. Temperatures are 63.9 K, 68.1 K, 70.4 K, 74.0 K, 77.8 K, 80.1 K, 
84.4 K, 91.4 K, 95.7 K, 99.5 K from above. was 79 K. (b). Calculated c- 
axis magnetoresistance at several temperatures according to the fiuctuation 



theory including the DOS contribution. After |185 



Figure 33. The theoretical dependence of the real part of conduc- 
tivity, normalized on the Drude normal conductivity, on uj/T, '^[a'ioj)] = 
Re [o"(co')] /o"'^. The dashed line refers to the a6-plane component of the con- 
ductivity tensor whose Drude normal conductivity is = N{0)e'^TVp. The 
solid line refers to the c-axis component whose Drude normal conductivity 
is al = al.Ps'^/vl. In this plot we have put Tr = Q.?>,Ef/T = 50,r = 
0.01, £ = 0.04, Tr^ = 4. 

Figure 34. The theoretical behavior of the c-axis component of con- 
ductivity frequency dependence, for different values of temperature, is shown. 
The solid line refers to e = 0.04; the dashed line refers to e = 0.06; the dot- 
dashed line refers to e = 0.08. Tr = 0.2 for all the curves. The other 
parameters of this plot are the same used in Fig. 33 

Figure 35. The opening of the pseudo-gap in the c-axis conductivity 
measurements on Y Ba2Cu40% samples |§] 

Figure 36. The theoretical dependence of Re [o"±(t^)/cr5.] on tu/T for 
different values of Tr. The solid line refers to Tr = 0.4; the dot-dashed line 
refers to Tr = 0.3; the dashed line refers to Tr = 0.2. The other parameters 
of this plot are the same used in Fig. 33. 

Figure 37. The Feynman diagrams for the fluctuation correction to ther- 
moelectric coefficient are shown. Shaded partial circles are impurity vertex 
corrections (|153|) , dashed curves with central crosses are additional impu- 
rity renormalizations, and shaded thick lines are additional impurity vertex 
corrections. 

Figure 38. Diagrams for the fluctuation contribution to the dynamic spin 
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susceptibility x+-- 

Figure 39. The transformation of the MT type diagram for the NMR 
relaxation rate in self-intersecting trajectory. 

Figure 40. The ^^Cu relaxation rates in zero field 2iy(0) (from NQR 
relaxation) and 2W{H) in a field of 5.9 T (from NMR relaxation of the 
-1/2 1/2 line) in the oriented powders of YBCO, with Tc(0) = 90.5 
K and Tf.{H) = 87.5 K. In the inset the relaxation rates, normalized with 
respect to W{H) = W{0) for T ^ Tc, are reported as a function of T/Tc. 

Figure 41. The Dyson equation in the ladder approximation for the par- 
ticle particle channel vertex X{q,LL!i,uj2) renormalized by impurities (shaded 
partial circles). Here the dashed line, as usually, presents the "impurity 
propagator" l/27riV(0)r (which is nothing else as the square of the scatter- 
ing amplitude avaraged over the solid angles), solid lines are one-electron 
Green functions (^Sp already avareged over impurities configuration. 

Figure 42. The Dyson equation in the ladder approximation for the fluc- 
tuation propagator. 
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